HypoxiaDB: a database of hypoxia-regulated proteins

Surgical Tumor Resection Deregulates Hallmarks of Cancer in Resected Tissue and the Surrounding Microenvironment

Surgery exposes tumor tissue to severe hypoxia and mechanical stress leading to rapid gene expression changes in the tumor and its microenvironment, which remain poorly characterized. We biopsied tumor and adjacent normal tissues from patients with breast (n = 81) and head/neck squamous cancers (HNSC; n = 10) at the beginning (A), during (B), and end of surgery (C). Tumor/normal RNA from 46/81 patients with breast cancer was subjected to mRNA-Seq using Illumina short-read technology, and from nine patients with HNSC to whole-transcriptome microarray with Illumina BeadArray. Pathways and genes involved in 7 of 10 known cancer hallmarks, namely, tumor-promoting inflammation (TNF-A, NFK-B, IL18 pathways), activation of invasion and migration (various extracellular matrix-related pathways, cell migration), sustained proliferative signaling (K-Ras Signaling), evasion of growth suppressors (P53 signaling, regulation of cell death), deregulating cellular energetics (response to lipid, secreted factors, and adipogenesis), inducing angiogenesis (hypoxia signaling, myogenesis), and avoiding immune destruction (CTLA4 and PDL1) were significantly deregulated during surgical resection (time points A vs. B vs. C). These findings were validated using NanoString assays in independent pre/intra/post-operative breast cancer samples from 48 patients. In a comparison of gene expression data from biopsy (analogous to time point A) with surgical resection samples (analogous to time point C) from The Cancer Genome Atlas study, the top deregulated genes were the same as identified in our analysis, in five of the seven studied cancer types. This study suggests that surgical extirpation deregulates the hallmarks of cancer in primary tumors and adjacent normal tissue across different cancers.

IMPLICATIONS

Surgery deregulates hallmarks of cancer in human tissue.

Sleep Disorder Kleine-Levin Syndrome (KLS) Joins the List of Polygenic Brain Disorders Associated with Obstetric Complications

Kleine-Levin Syndrome is a rare neurological disorder with onset typically during adolescence that is characterized by recurrent episodes of hypersomnia, behavioral changes, and cognitive abnormalities, in the absence of structural changes in neuroimaging. As for many functional brain disorders, the exact disease mechanism in Kleine-Levin Syndrome is presently unknown, preventing the development of specific treatment approaches or protective measures. Here we review the pathophysiology and genetics of this functional brain disorder and then present a specific working hypothesis. A neurodevelopmental mechanism has been suspected based on associations with obstetric complications. Recent studies have focused on genetic factors whereby the first genome-wide association study (GWAS) in Kleine-Levin Syndrome has defined a linkage at the TRANK1 locus. A Gene x Environment interaction model involving obstetric complications was proposed based on concepts developed for other functional brain disorders. To stimulate future research, we here performed annotations of the genes under consideration for Kleine-Levin Syndrome in relation to factors expected to be associated with obstetric complications. Annotations used data-mining of gene/protein lists related to for hypoxia, ischemia, and vascular factors and targeted literature searches. Tentative links for TRANK1, four additional genes in the TRANK1 locus, and LMOD3-LMO2 are described. Protein interaction data for TRANK1 indicate links to CBX2, CBX4, and KDM3A, that in turn can be tied to hypoxia. Taken together, the neurological sleep disorder, Kleine-Levin Syndrome, shows genetic and mechanistic overlap with well analyzed brain disorders such as schizophrenia, autism spectrum disorder and ADHD in which polygenic predisposition interacts with external events during brain development, including obstetric complications.

Comparative genomics provides insights into the aquatic adaptations of mammals

Divergent lineages can respond to common environmental factors through convergent processes involving shared genomic components or pathways, but the molecular mechanisms are poorly understood. Here, we provide genomic resources and insights into the evolution of mammalian lineages adapting to aquatic life. Our data suggest convergent evolution, for example, in association with thermoregulation through genes associated with a surface heat barrier (NFIA) and internal heat exchange (SEMA3E). Combined with the support of previous reports showing that the UCP1 locus has been lost in many marine mammals independently, our results suggest that the thermostatic strategy of marine mammals shifted from enhancing heat production to limiting heat loss.

The ancestors of marine mammals once roamed the land and independently committed to an aquatic lifestyle. These macroevolutionary transitions have intrigued scientists for centuries. Here, we generated high-quality genome assemblies of 17 marine mammals (11 cetaceans and six pinnipeds), including eight assemblies at the chromosome level. Incorporating previously published data, we reconstructed the marine mammal phylogeny and population histories and identified numerous idiosyncratic and convergent genomic variations that possibly contributed to the transition from land to water in marine mammal lineages. Genes associated with the formation of blubber (NFIA), vascular development (SEMA3E), and heat production by brown adipose tissue (UCP1) had unique changes that may contribute to marine mammal thermoregulation. We also observed many lineage-specific changes in the marine mammals, including genes associated with deep diving and navigation. Our study advances understanding of the timing, pattern, and molecular changes associated with the evolution of mammalian lineages adapting to aquatic life.

Genomic and functional evidence reveals convergent evolution in fishes on the Tibetan Plateau

High-altitude environments are strong drivers of adaptive evolution in endemic organisms. However, little is known about the genetic mechanisms of convergent adaptation among different lineages, especially in fishes. There are three independent fish groups on the Tibetan Plateau: Tibetan Loaches, Schizothoracine fishes and Glyptosternoid fishes; all are well adapted to the harsh environmental conditions. They represent an excellent example of convergent evolution but with an unclear genetic basis. We used comparative genomic analyses between Tibetan fishes and fishes from low altitudes and detected genomic signatures of convergent evolution in fishes on the Tibetan Plateau. The Tibetan fishes exhibited genome-wide accelerated evolution in comparison with a control set of fishes from low altitudes. A total of 368 positively selected genes were identified in Tibetan fishes, which were enriched in functional categories related to energy metabolism and hypoxia response. Widespread parallel amino acid substitutions were detected among the Tibetan fishes and a subset of these substitutions occurred in positively selected genes associated with high-altitude adaptation. Functional assays suggested that von Hippel-Lindau (VHL) tumour suppressor genes from Tibetan fishes enhance hypoxia-inducible factor (HIF) activity convergently under hypoxia compared to low-altitude fishes. The results provide genomic and functional evidence supporting convergent genetic mechanisms for high-altitude adaptation in fishes on the Tibetan Plateau.

Gene network analysis to determine the effect of hypoxia-associated genes on brain damages and tumorigenesis using an avian model

BACKGROUND

Hypoxia refers to the condition of low oxygen pressure in the atmosphere and characterization of response to hypoxia as a biological complex puzzle, is challenging. Previously, we carried out a comparative genomic study by whole genome resequencing of highland and lowland Iranian native chickens to identify genomic variants associated with hypoxia conditions. Based on our previous findings, we used chicken as a model and the identified hypoxia-associated genes were converted to human's orthologs genes to construct the informative gene network. The main goal of this study was to visualize the features of diseases due to hypoxia-associated genes by gene network analysis.

RESULTS

It was found that hypoxia-associated genes contained several gene networks of disorders such as Parkinson, Alzheimer, cardiomyopathy, drug toxicity, and cancers. We found that biological pathways are involved in mitochondrion dysfunctions including peroxynitrous acid production denoted in brain injuries. Lewy body and neuromelanin were reported as key symptoms in Parkinson disease. Furthermore, calmodulin, and amyloid precursor protein were detected as leader proteins in Alzheimer's diseases. Dexamethasone was reported as the candidate toxic drug under the hypoxia condition that implicates diabetes, osteoporosis, and neurotoxicity. Our results suggested DNA damages caused by the high doses of UV radiation in high-altitude conditions, were associated with breast cancer, ovarian cancer, and colorectal cancer.

CONCLUSIONS

Our results showed that hypoxia-associated genes were enriched in several gene networks of disorders including Parkinson, Alzheimer, cardiomyopathy, drug toxicity, and different types of cancers. Furthermore, we suggested, UV radiation and low oxygen conditions in high-altitude regions may be responsible for the variety of human diseases.

A functional map of genomic HIF1α-DNA complexes in the eye lens revealed through multiomics analysis

BACKGROUND

During eye lens development the embryonic vasculature regresses leaving the lens without a direct oxygen source. Both embryonically and throughout adult life, the lens contains a decreasing oxygen gradient from the surface to the core that parallels the natural differentiation of immature surface epithelial cells into mature core transparent fiber cells. These properties of the lens suggest a potential role for hypoxia and the master regulator of the hypoxic response, hypoxia-inducible transcription factor 1 (HIF1), in the regulation of genes required for lens fiber cell differentiation, structure and transparency. Here, we employed a multiomics approach combining CUT&RUN, RNA-seq and ATACseq analysis to establish the genomic complement of lens HIF1α binding sites, genes activated or repressed by HIF1α and the chromatin states of HIF1α-regulated genes.

RESULTS

CUT&RUN analysis revealed 8375 HIF1α-DNA binding complexes in the chick lens genome. One thousand one hundred ninety HIF1α-DNA binding complexes were significantly clustered within chromatin accessible regions (χ2 test p < 1 × 10- 55) identified by ATACseq. Formation of the identified HIF1α-DNA complexes paralleled the activation or repression of 526 genes, 116 of which contained HIF1α binding sites within 10kB of the transcription start sites. Some of the identified HIF1α genes have previously established lens functions while others have novel functions never before examined in the lens. GO and pathway analysis of these genes implicate HIF1α in the control of a wide-variety of cellular pathways potentially critical for lens fiber cell formation, structure and function including glycolysis, cell cycle regulation, chromatin remodeling, Notch and Wnt signaling, differentiation, development, and transparency.

CONCLUSIONS

These data establish the first functional map of genomic HIF1α-DNA complexes in the eye lens. They identify HIF1α as an important regulator of a wide-variety of genes previously shown to be critical for lens formation and function and they reveal a requirement for HIF1α in the regulation of a wide-variety of genes not yet examined for lens function. They support a requirement for HIF1α in lens fiber cell formation, structure and function and they provide a basis for understanding the potential roles and requirements for HIF1α in the development, structure and function of more complex tissues.

Transcriptome Analysis of Hypoxic Lymphatic Endothelial Cells Indicates Their Potential to Contribute to Extracellular Matrix Rearrangement

Lymphedema (LE) affects millions of people worldwide. It is a chronic progressive disease with massive development of fibrosclerosis when untreated. There is no pharmacological treatment of lymphedema. The disease is associated with swelling of the interstitium of the affected organ, mostly arm or leg, impressive development of adipose tissue, fibrosis and sclerosis with accumulation of huge amounts of collagen, and Papillomatosis cutis. Malnutrition and reduced oxygenation of the affected tissues is a hallmark of lymphedema. Here, we investigated if the hypoxia of lymphatic endothelial cells (LECs) might contribute to fibrosis. We applied RNASeq and qPCR to study the concordant changes of the exome of three human foreskin-derived LEC isolates after 4 days of hypoxia (1% O₂) vs. normoxia (21% O₂). Of the approximately 16,000 genes expressed in LECs, 162 (1%) were up- or down-regulated by hypoxia. Of these, 21 genes have important functions in the production or modification of the extracellular matrix (ECM). In addition to the down-regulation of elastin, we found up-regulation of druggable enzymes and regulators such as the long non-coding RNA H19, inter-alpha-trypsin inhibitor heavy chain family member 5 (ITIH5), lysyl-oxidase (LOX), prolyl 4-hydroxylase subunit alpha 1 (P4HA1), procollagen-lysine 2-oxoglutarate 5-dioxygenase 2 (PLOD2), and others that are discussed in the paper. Initial lymphatics do not produce a continuous basement membrane; however, our study shows that hypoxic LECs have an unexpectedly high ability to alter the ECM.

Out of Tibet: Genomic Perspectives on the Evolutionary History of Extant Pikas

Pikas are widely distributed in the Northern Hemisphere and are highly adapted to cold and alpine environments. They are one of the most complex and problematic groups in mammalian systematics, and the origin and evolutionary history of extant pikas remain controversial. In this study, we sequenced the whole coding sequences of 105 pika samples (29 named species and 1 putative new species) and obtained DNA data for more than 10,000 genes. Our phylogenomic analyses recognized four subgenera of extant pikas: Alienauroa, Conothoa, Ochotona, and Pika. The interrelationships between the four subgenera were strongly resolved as (Conothoa, (Alienauroa, (Ochotona, Pika))), with the mountain group Conothoa being the sister group of all other pikas. Our divergence time and phylogeographic analyses indicated that the last common ancestor of extant pikas first occurred on in the middle Miocene, ∼14 Ma. The emergence of opportunities related to the climate, food supply, and spreading paths in concert promoted the dispersal of pikas from the Qinghai-Tibetan Plateau (QTP) to other parts of Eurasia and North America. We found that the genes that were positively selected in the early evolution of pikas were most concentrated in functional categories related to cold tolerance. These results suggest that the QTP may have served as a training ground for cold tolerance in early pikas, which gives pikas a great advantage when the climate continued to cool after the middle Miocene. Our study highlights the importance of the QTP as a center of origin for many cold-adapted animals.

Deferoxamine mesylate improves splicing and GAA activity of the common c.-32-13T>G allele in late-onset PD patient fibroblasts

Pompe disease (PD) is an autosomal recessive lysosomal storage disorder due to deficient activity of the acid alpha glucosidase enzyme (GAA). As a consequence of the enzymatic defect, undigested glycogen accumulates within lysosomes. Most patients affected by the late-onset (LO) phenotype carry in at least one allele the c.-32-13T>G variant, which leads to exon 2 exclusion from the pre-mRNA. These patients display a variable and suboptimal response to enzyme replacement therapy. To identify novel therapeutic approaches, we developed a fluorescent GAA exon 2 splicing assay and screened a library of US Food and Drug Administration (FDA)-approved compounds. This led to the identification of several drugs able to restore normal splicing. Among these, we further validated the effects of the iron chelator deferoxamine (Defe) in c.-32-13T>G fibroblasts. Defe treatment resulted in a 2-fold increase of GAA exon 2 inclusion and a 40% increase in enzymatic activity. Preliminary results suggest that this effect is mediated by the regulation of iron availability, at least partially. RNA-seq experiments also showed that Defe might shift the balance of splicing factor levels toward a profile promoting GAA exon 2 inclusion. This work provides the basis for drug repurposing and development of new chemically modified molecules aimed at improving the clinical outcome in LO-PD patients.

The role of long intergenic non-coding RNA for kinase activation (LINK-A) as an oncogene in non-small cell lung carcinoma

The oncogenic role of long intergenic non-coding RNA for kinase activation (LINK-A) has been appraised in triple-negative breast cancer. However, the molecular function of LINK-A is still unclear in most cancers including lung cancer. The present study aimed to evaluate the impact of down-regulation of LINK-A in A549 and Calu-3 cell lines as cellular models of non-small cell lung carcinoma (NSCLC). We used the RNA interference system to knock down LINK-A. LINK-A expression was significantly reduced by siRNA transfection in A549 and Calu-3 cell lines. LINK-A down-regulation significantly reduced cell viability, colony-forming ability and cell migration, as measured by MTT, colony formation and invasion assays. Finally, cell cycle analysis and Annexin-V/7AAD staining indicated that apoptosis was influenced by LINK-A silencing. Taken together, LINK-A can be proposed as an oncogene in NSCLC.

Comparative transcriptomic analysis of fireflies (Coleoptera: Lampyridae) to explore the molecular adaptations to fresh water

Aquatic insects are well adapted to freshwater environments, but the molecular basis of these adaptations remains largely unknown. Most firefly species (Coleoptera: Lampyridae) are terrestrial, but the larvae of several species are aquatic. Here, larval and adult transcriptomes from Aquatica leii (freshwater) and Lychnuris praetexta (terrestrial) were generated to test whether the genes associated with metabolic efficiency and morphology have undergone adaptive evolution to fresh water. The aquatic fireflies had a significantly lower ratio of nonsynonymous to synonymous substitutions than the terrestrial insects, indicating a genomewide evolutionary constraint in the aquatic fireflies. We identified 341 fast-evolving genes and 116 positively selected genes in the aquatic fireflies. Of these, 76 genes exhibiting both fast evolution and positive selection were primarily involved in ATP production, energy metabolism and the hypoxia response. We identified 7,271 differentially expressed genes (DEGs) in A. leii (adults versus larvae) and 8,309 DEGs in L. praetexta (adults versus larvae). DEGs specific to the aquatic firefly (n = 1,445) were screened via interspecific comparisons (A. leii versus L. praetexta) and were significantly enriched for genes involved in metabolic efficiency (e.g., ATP production, hypoxia, and immune responses) and certain aspects of morphology (e.g., cuticle chitin, tracheal and compound eye morphology). These results indicate that sequence and expression-level changes in genes associated with both metabolic efficiency and morphological attributes related to the freshwater lifestyle contributed to freshwater adaptation in fireflies. This study provides new insights into the molecular mechanisms of aquatic adaptation in insects.

Analysis of GWAS-Derived Schizophrenia Genes for Links to Ischemia-Hypoxia Response of the Brain

Obstetric complications (OCs) can induce major adverse conditions for early brain development and predispose to mental disorders, including schizophrenia (SCZ). We previously hypothesized that SCZ candidate genes respond to ischemia-hypoxia as part of OCs which impacts neurodevelopment. We here tested for an overlap between SCZ genes from genome-wide association study (GWAS) (n=458 genes from 145 loci of the most recent GWAS dataset in SCZ) and gene sets for ischemia-hypoxia response. Subsets of SCZ genes were related to (a) mutation-intolerant genes (LoF database), (b) role in monogenic disorders of the nervous system (OMIM, manual annotations), and (c) synaptic function (SynGO). Ischemia-hypoxia response genes of the brain (IHR genes, n=1,629), a gene set from RNAseq in focal brain ischemia (BH, n=2,449) and genes from HypoxiaDB (HDB, n=2,289) were overlapped with the subset of SCZ genes and tested for enrichment with Chi-square tests (p < 0.017). The SCZ GWAS dataset was enriched for LoF (n=112; p=0.0001), and the LoF subset was enriched for IHR genes (n=25; p=0.0002), BH genes (n=35; p=0.0001), and HDB genes (n=23; p=0.0005). N=96 genes of the SCZ GWAS dataset (21%) could be linked to a monogenic disorder of the nervous system whereby IHR genes (n=19, p=0.008) and BH genes (n=23; p=0.002) were found enriched. N=46 synaptic genes were found in the SCZ GWAS gene set (p=0.0095) whereby enrichments for IHR genes (n=20; p=0.0001) and BH genes (n=13; p=0.0064) were found. In parallel, detailed annotations of SCZ genes for a role of the hypoxia-inducible factors (HIFs) identified n=33 genes of high interest. Genes from SCZ GWAS were enriched for mutation-intolerant genes which in turn were strongly enriched for three sets of genes for the ischemia-hypoxia response that may be invoked by OCs. A subset of one fifth of SCZ genes has established roles in monogenic disorders of the nervous system which was enriched for two gene sets related to ischemia-hypoxia. SCZ genes related to synaptic functions were also related to ischemia-hypoxia. Variants of SCZ genes interacting with ischemia-hypoxia provide a specific starting point for functional and genomic studies related to OCs.

Targeting Cyclooxygenase-2 in Pheochromocytoma and Paraganglioma: Focus on Genetic Background

Cyclooxygenase 2 (COX-2) is a key enzyme of the tumorigenesis-inflammation interface and can be induced by hypoxia. A pseudohypoxic transcriptional signature characterizes pheochromocytomas and paragangliomas (PPGLs) of the cluster I, mainly represented by tumors with mutations in von Hippel–Lindau (VHL), endothelial PAS domain-containing protein 1 (EPAS1), or succinate dehydrogenase (SDH) subunit genes. The aim of this study was to investigate a possible association between underlying tumor driver mutations and COX-2 in PPGLs. COX-2 gene expression and immunoreactivity were examined in clinical specimens with documented mutations, as well as in spheroids and allografts derived from mouse pheochromocytoma (MPC) cells. COX-2 in vivo imaging was performed in allograft mice. We observed significantly higher COX-2 expression in cluster I, especially in VHL-mutant PPGLs, however, no specific association between COX-2 mRNA levels and a hypoxia-related transcriptional signature was found. COX-2 immunoreactivity was present in about 60% of clinical specimens as well as in MPC spheroids and allografts. A selective COX-2 tracer specifically accumulated in MPC allografts. This study demonstrates that, although pseudohypoxia is not the major determinant for high COX-2 levels in PPGLs, COX-2 is a relevant molecular target. This potentially allows for employing selective COX-2 inhibitors as targeted chemotherapeutic agents and radiosensitizers. Moreover, available models are suitable for preclinical testing of these treatments.

Genome-wide detection of copy number variations in polled yak using the Illumina BovineHD BeadChip

Background

Copy number variations (CNVs), which are genetic variations present throughout mammalian genomes, are a vital source of phenotypic variation that can lead to the development of unique traits. In this study we used the Illunima BovineHD BeadChip to conduct genome-wide detection of CNVs in 215 polled yaks.

Results

A total of 1066 CNV regions (CNVRs) were detected with a total length of 181.6 Mb, comprising ~ 7.2% of the bovine autosomal genome. The size of these CNVRs ranged from 5.53 kb to 1148.45 kb, with an average size of 170.31 kb. Eight out of nine randomly chosen CNVRs were successfully validated by qPCR. A functional enrichment analysis of the CNVR-associated genes indicated their relationship to a number of molecular adaptations that enable yaks to thrive at high altitudes. One third of the detected CNVRs were mapped to QTLs associated with six classes of economically important traits, indicating that these CNVRs may play an important role in variations of these traits.

Conclusions

Our genome-wide yak CNV map may thus provide valuable insights into both the molecular mechanisms of high altitude adaptation and the potential genomic basis of economically important traits in yak.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5759-1) contains supplementary material, which is available to authorized users.

Culture surfaces induce hypoxia-regulated genes in human mesenchymal stromal cells

Culturing human Mesenchymal stromal cells (hMSCs) in vitro in hypoxic conditions resulted in reduced senescence, enhanced pluripotency and altered proliferation rate. It has been known that in vitro hypoxia affects expression of cell surface proteins. However, the impact of culture surfaces on the hypoxia-regulated genes (HRG) have not yet been reported. This study utilized Next-Generation sequencing to analyse the changes in the gene expression levels of HRG for hMSCs cultured on different culture surfaces. The samples, which were cultured on four different synthesized surfaces (treatments) and tissue culture plate (control), resulted in a difference in growth rate. The sequencing results revealed that the transcription of a number of key genes involved in regulating hypoxic functions were significantly altered, including HIF2A, a marker for potency, differentiation, and various cellular functions. Significant alternations in the expression levels of previously reported oxygen-sensitive surface proteins were detected in this study, some of which closely correlate with the expression levels of HIF2A. Our analysis of the hMSCs transcriptome and HRG mapped out a list of genes encoding surface proteins which may directly regulate or be regulated by HIF2A. The findings from this study showed that culture surfaces have an impact on regulating the expression profile of HRG. Therefore, novel culture surfaces may be designed to selectively activate HIF2A and other HRG and pathways under in vitro normoxia. The understanding of the crosstalk between the regulating genes of hypoxia and culture surfaces may be utilized to strengthen desired hypoxic functions.

Expression-based analyses indicate a central role for hypoxia in driving tumor plasticity through microenvironment remodeling and chromosomal instability

Can transcriptomic alterations drive the evolution of tumors? We asked if changes in gene expression found in all patients arise earlier in tumor development and can be relevant to tumor progression. Our analyses of non-mutated genes from the non-amplified regions of the genome of 158 triple-negative breast cancer (TNBC) cases identified 219 exclusively expression-altered (EEA) genes that may play important role in TNBC. Phylogenetic analyses of these genes predict a "punctuated burst" of multiple gene upregulation events occurring at early stages of tumor development, followed by minimal subsequent changes later in tumor progression. Remarkably, this punctuated burst of expressional changes is instigated by hypoxia-related molecular events, predominantly in two groups of genes that control chromosomal instability (CIN) and those that remodel tumor microenvironment (TME). We conclude that alterations in the transcriptome are not stochastic and that early-stage hypoxia induces CIN and TME remodeling to permit further tumor evolution.

Species groups distributed across elevational gradients reveal convergent and continuous genetic adaptation to high elevations

Although many cases of genetic adaptations to high elevations have been reported, the processes driving these modifications and the pace of their evolution remain unclear. Many high-elevation adaptations (HEAs) are thought to have arisen in situ as populations rose with growing mountains. In contrast, most high-elevation lineages of the Qinghai-Tibetan Plateau appear to have colonized from low-elevation areas. These lineages provide an opportunity for studying recent HEAs and comparing them with ancestral low-elevation alternatives. Herein, we compare four frogs (three species of Nanorana and a close lowland relative) and four lizards (Phrynocephalus) that inhabit a range of elevations on or along the slopes of the Qinghai-Tibetan Plateau. The sequential cladogenesis of these species across an elevational gradient allows us to examine the gradual accumulation of HEA at increasing elevations. Many adaptations to high elevations appear to arise gradually and evolve continuously with increasing elevational distributions. Numerous related functions, especially DNA repair and energy metabolism pathways, exhibit rapid change and continuous positive selection with increasing elevations. Although the two studied genera are distantly related, they exhibit numerous convergent evolutionary changes, especially at the functional level. This functional convergence appears to be more extensive than convergence at the individual gene level, although we found 32 homologous genes undergoing positive selection for change in both high-elevation groups. We argue that species groups distributed along a broad elevational gradient provide a more powerful system for testing adaptations to high-elevation environments compared with studies that compare only pairs of high-elevation versus low-elevation species.

BHLHE40 confers a pro-survival and pro-metastatic phenotype to breast cancer cells by modulating HBEGF secretion

Background

Metastasis is responsible for a significant number of breast cancer-related deaths. Hypoxia, a primary driving force of cancer metastasis, induces the expression of BHLHE40, a transcription regulator. This study aimed to elucidate the function of BHLHE40 in the metastatic process of breast cancer cells.

Methods

To define the role of BHLHE40 in breast cancer, BHLHE40 expression was knocked down by a lentiviral construct expressing a short hairpin RNA against BHLHE40 or knocked out by the CRISPR/Cas9 editing system. Orthotopic xenograft and experimental metastasis (tail vein injection) mouse models were used to analyze the role of BHLHE40 in lung metastasis of breast cancer. Global gene expression analysis and public database mining were performed to identify signaling pathways regulated by BHLHE40 in breast cancer. The action mechanism of BHLHE40 was examined by chromatin immunoprecipitation (ChIP), co-immunoprecipitation (CoIP), exosome analysis, and cell-based assays for metastatic potential.

Results

BHLHE40 knockdown significantly reduced primary tumor growth and lung metastasis in orthotopic xenograft and experimental metastasis models of breast cancer. Gene expression analysis implicated a role of BHLHE40 in transcriptional activation of heparin-binding epidermal growth factor (HBEGF). ChIP and CoIP assays revealed that BHLHE40 induces HBEGF transcription by blocking DNA binding of histone deacetylases (HDAC)1 and HDAC2. Cell-based assays showed that HBEGF is secreted through exosomes and acts to promote cell survival and migration. Public databases provided evidence linking high expression of BHLHE40 and HBEGF to poor prognosis of triple-negative breast cancer.

Conclusion

This study reveals a novel role of BHLHE40 in promoting tumor cell survival and migration by regulating HBEGF secretion.

Role of miRNAs in hypoxia-related disorders

Hypoxia is a complex pathophysiological condition. The physiological and molecular responses to this stress have been extensively studied. However, the management of its ill effects still poses a challenge to clinicians. MicroRNAs (miRNAs) are short non-coding RNA molecules that control post-transcriptional gene expression. The regulatory role of miRNAs in hypoxic environments has been studied in many hypoxia-related disorders, however a comprehensive compilation and analysis of all data and the significance of miRNAs in hypoxia adaption is still lacking. This review summarizes the miRNAs related to various hypoxia-related disorders and highlights the computational approaches to study them. This would help in designing novel strategies toward efficient management of hypoxia-related disorders.

Programmed mitophagy is essential for the glycolytic switch during cell differentiation

Retinal ganglion cells (RGCs) are the sole projecting neurons of the retina and their axons form the optic nerve. Here, we show that embryogenesis-associated mouse RGC differentiation depends on mitophagy, the programmed autophagic clearance of mitochondria. The elimination of mitochondria during RGC differentiation was coupled to a metabolic shift with increased lactate production and elevated expression of glycolytic enzymes at the mRNA level. Pharmacological and genetic inhibition of either mitophagy or glycolysis consistently inhibited RGC differentiation. Local hypoxia triggered expression of the mitophagy regulator BCL2/adenovirus E1B 19-kDa-interacting protein 3-like (BNIP3L, best known as NIX) at peak RGC differentiation. Retinas from NIX-deficient mice displayed increased mitochondrial mass, reduced expression of glycolytic enzymes and decreased neuronal differentiation. Similarly, we provide evidence that NIX-dependent mitophagy contributes to mitochondrial elimination during macrophage polarization towards the proinflammatory and more glycolytic M1 phenotype, but not to M2 macrophage differentiation, which primarily relies on oxidative phosphorylation. In summary, developmentally controlled mitophagy promotes a metabolic switch towards glycolysis, which in turn contributes to cellular differentiation in several distinct developmental contexts.

Genetic Adaptation of Schizothoracine Fish to the Phased Uplifting of the Qinghai-Tibetan Plateau

Many species of Schizothoracine, a subfamily of Cyprinidae, are highly endemic to the Qinghai–Tibetan Plateau (QTP). To characterize the adaptive changes associated with the Schizothoracine expansion at high altitudes, we sequenced tissue transcriptomes of two highland and two subhighland Schizothoracines and analyzed gene evolution patterns by comparing with lowland cyprinids. Phylogenetic tree reconstruction and divergence time estimation indicated that the common ancestor of Schizothoracine fish lived ∼32.7 million years ago (MYA), coinciding with the timing of the first phase of QTP uplifting. Both high- and subhigh-Schizothoracines demonstrated elevated dN/dS ratios in the protein-coding genes compared to lowland cyprinids, from which some biological processes implicated in altitude adaptation were commonly identified. On the other hand, the highland and subhighland lineages presented drastically divergent landscapes of positively selected genes (PSGs), enriched with very different gene ontology (GO) profiles, including those in “sensory organ morphogenesis,” “regulation of protein ubiquitination,” “blood circulation,” and “blood vessel development.” These results indicated different selection pressures imposed on the highland and subhighland lineages of the Schizothoracine subfamily, with a higher number of genes in the high-altitude species involved in adaptations such as sensory perception, blood circulation, and protein metabolism. Our study indicated divergent genetic adaptations in the aquatic species facing the phased uplifting of QTP.

Transcription factor HIF1A: downstream targets, associated pathways, polymorphic hypoxia response element (HRE) sites, and initiative for standardization of reporting in scientific literature

Hypoxia-inducible factor-1α (HIF-1α) has crucial role in adapting cells to hypoxia through expression regulation of many genes. Identification of HIF-1α target genes (HIF-1α-TGs) is important for understanding the adapting mechanism. The aim of the present study was to collect known HIF-1α-TGs and identify their associated pathways. Targets and associated genomics data were retrieved using PubMed, WoS ( http://apps.webofknowledge.com/ ), HGNC ( http://www.genenames.org/ ), NCBI ( http://www.ncbi.nlm.nih.gov/ ), Ensemblv.84 ( http://www.ensembl.org/index.html ), DAVID Bioinformatics Resources ( https://david.ncifcrf.gov /), and Disease Ontology database ( http://disease-ontology.org/ ). From 51 papers, we collected 98 HIF-1α TGs found to be associated with 20 pathways, including metabolism of carbohydrates and pathways in cancer. Reanalysis of genomic coordinates of published HREs (hypoxia response elements) revealed six polymorphisms within HRE sites (HRE-SNPs): ABCG2, ACE, CA9, and CP. Due to large heterogeneity of results presentation in scientific literature, we also propose a first step towards reporting standardization of HIF-1α-target interactions consisting of ten relevant data types. Suggested minimal checklist for reporting will enable faster development of a complete catalog of HIF-1α-TGs, data sharing, bioinformatics analyses, and setting novel more targeted hypotheses. The proposed format for data standardization is not yet complete but presents a baseline for further optimization of the protocol with additional details, for example, regarding the experimental validation.

Evidence for Adaptation to the Tibetan Plateau Inferred from Tibetan Loach Transcriptomes

Triplophysa fishes are the primary component of the fish fauna on the Tibetan Plateau and are well adapted to the high-altitude environment. Despite the importance of Triplophysa fishes on the plateau, the genetic mechanisms of the adaptations of these fishes to this high-altitude environment remain poorly understood. In this study, we generated the transcriptome sequences for three Triplophysa fishes, that is, Triplophysa siluroides, Triplophysa scleroptera, and Triplophysa dalaica, and used these and the previously available transcriptome and genome sequences from fishes living at low altitudes to identify potential genetic mechanisms for the high-altitude adaptations in Triplophysa fishes. An analysis of 2,269 orthologous genes among cave fish (Astyanax mexicanus), zebrafish (Danio rerio), large-scale loach (Paramisgurnus dabryanus), and Triplophysa fishes revealed that each of the terminal branches of the Triplophysa fishes had a significantly higher ratio of nonsynonymous to synonymous substitutions than that of the branches of the fishes from low altitudes, which provided consistent evidence for genome-wide rapid evolution in the Triplophysa genus. Many of the GO (Gene Ontology) categories associated with energy metabolism and hypoxia response exhibited accelerated evolution in the Triplophysa fishes compared with the large-scale loach. The genes that exhibited signs of positive selection and rapid evolution in the Triplophysa fishes were also significantly enriched in energy metabolism and hypoxia response categories. Our analysis identified widespread Triplophysa-specific nonsynonymous mutations in the fast evolving genes and positively selected genes. Moreover, we detected significant evidence of positive selection in the HIF (hypoxia-inducible factor)-1A and HIF-2B genes in Triplophysa fishes and found that the Triplophysa-specific nonsynonymous mutations in the HIF-1A and HIF-2B genes were associated with functional changes. Overall, our study provides new insights into the adaptations and evolution of fishes in the high-altitude environment of the Tibetan Plateau and complements previous findings on the adaptations of mammals and birds to high altitudes.

Comprehensive transcriptome analysis reveals accelerated genic evolution in a Tibet fish, Gymnodiptychus pachycheilus

Elucidating the genetic mechanisms of organismal adaptation to the Tibetan Plateau at a genomic scale can provide insights into the process of adaptive evolution. Many highland species have been investigated and various candidate genes that may be responsible for highland adaptation have been identified. However, we know little about the genomic basis of adaptation to Tibet in fishes. Here, we performed transcriptome sequencing of a schizothoracine fish (Gymnodiptychus pachycheilus) and used it to identify potential genetic mechanisms of highland adaptation. We obtained totally 66,105 assembled unigenes, of which 7,232 were assigned as putative one-to-one orthologs in zebrafish. Comparative gene annotations from several species indicated that at least 350 genes lost and 41 gained since the divergence between G. pachycheilus and zebrafish. An analysis of 6,324 orthologs among zebrafish, fugu, medaka, and spotted gar identified consistent evidence for genome-wide accelerated evolution in G. pachycheilus and only the terminal branch of G. pachycheilus had an elevated Ka/Ks ratio than the ancestral branch. Many functional categories related to hypoxia and energy metabolism exhibited rapid evolution in G. pachycheilus relative to zebrafish. Genes showing signature of rapid evolution and positive selection in the G. pachycheilus lineage were also enriched in functions associated with energy metabolism and hypoxia. The first genomic resources for fish in the Tibetan Plateau and evolutionary analyses provided some novel insights into highland adaptation in fishes and served as a foundation for future studies aiming to identify candidate genes underlying the genetic bases of adaptation to Tibet in fishes.