Physiological and genetic convergence supports hypoxia resistance in high-altitude songbirds
Authors:
Ying Xiong aff001; Liqing Fan aff003; Yan Hao aff001; Yalin Cheng aff001; Yongbin Chang aff001; Jing Wang aff001; Haiyan Lin aff001; Gang Song aff001; Yanhua Qu aff001; Fumin Lei aff001
Authors place of work:
Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
aff001; University of Chinese Academy of Sciences, Beijing, China
aff002; National Forest Ecosystem Observation & Research Station of Nyingchi Tibet, Institute of Plateau Ecology, Tibet Agriculture & Animal Husbandry University, Linzhi City, China
aff003; Key Laboratory of Forest Ecology in Tibet Plateau (Tibet Agriculture & Animal Husbandry University), Ministry of Education, Linzhi City, China
aff004; Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, China
aff005
Published in the journal:
Physiological and genetic convergence supports hypoxia resistance in high-altitude songbirds. PLoS Genet 16(12): e1009270. doi:10.1371/journal.pgen.1009270
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pgen.1009270
Summary
Skeletal muscle plays a central role in regulating glucose uptake and body metabolism; however, highland hypoxia is a severe challenge to aerobic metabolism in small endotherms. Therefore, understanding the physiological and genetic convergence of muscle hypoxia tolerance has a potential broad range of medical implications. Here we report and experimentally validate a common physiological mechanism across multiple high-altitude songbirds that improvement in insulin sensitivity contributes to glucose homeostasis, low oxygen consumption, and relative activity, and thus increases body weight. By contrast, low-altitude songbirds exhibit muscle loss, glucose intolerance, and increase energy expenditures under hypoxia. This adaptive mechanism is attributable to convergent missense mutations in the BNIP3L gene, and METTL8 gene that activates MEF2C expression in highlanders, which in turn increases hypoxia tolerance. Together, our findings from wild high-altitude songbirds suggest convergent physiological and genetic mechanisms of skeletal muscle in hypoxia resistance, which highlights the potentially medical implications of hypoxia-related metabolic diseases.
Keywords:
Bird genetics – Bird physiology – Birds – Glucose – Hypoxia – Skeletal muscles
Introduction
Hypoxia, including ambient hypoxia (high-altitude exposure) [1] and pathological hypoxia (diseases) [2], induces metabolic alterations, muscle loss and insulin resistance [3–5]. At high altitude, hypoxia resulting from the low barometric pressure is a significantly severe challenge to aerobic exercise and thermogenesis in small endotherms, because a high rate of O2 flux should be concurrently sustained to thermogenesis in cold temperature [6]. Skeletal muscle plays a central role in regulating glucose uptake and body metabolism [7]. Therefore, muscle physiological modifications, including typically greater capillary density, more oxidative fiber, and higher proportion of subsarcolemmal mitochondria have been found in highland animals [6, 8, 9]. However, the physiological and genetic changes in skeletal muscles supporting hypoxia resistance are poorly understood, particularly for high-altitude songbirds.
In evolutionary biology, a growing body of literature has shown convergent physiological and genetic adaptations to high elevations across variant species. For example, high-altitude mice have evolved a physiological strategy of increasing carbohydrate utilization to economize oxygen use [10]. Convergent gene expression and genetic mechanisms for hypoxia adaptation have been elucidated in a number of mammals [11], birds [12], reptiles [13], and amphibians [14]. Generally, similar phenotypes occurring in different species under the same environmental stressors are proposed to depend on common physiological and genetic mechanisms. However, most physiological convergence lacks genetic evidence or genetic convergence is confirmed merely at the protein or cellular levels in vitro [10, 13]. To date, whether convergent gene expression shifts or mutants for hypoxia adaptation indeed drive a similar physiology change and subsequently a morphological phenotype in whole organisms remains unclear.
Here we mainly focused on snowfinches (Onychostruthus taczanowskii, and Pyrgilauda ruficollis) and tree sparrow (Passer montanus) belonging to the Old World sparrows (Passeridae). Snowfinches are endemic natives with a strict elevation distribution ranging from 3,500 to 5,100m on the Qinghai-Tibet Plateau (QTP), whereas tree sparrows are an introduced colonizer on the QTP despite its altitudinal distribution from sea level to 4,400m [15, 16]. We integrated behavioral, physiological and transcriptomic data across altitudinal populations and species to chiefly reveal common physiological and genetic mechanisms of muscle hypoxia tolerance. We report and experimentally validate using RNA interference in vivo a convergent mechanism that improves insulin sensitivity based on natural selection of standing genetic variations are essential to muscle hypoxia resistance in high-altitude songbirds, which provides novel and valuable insights into understanding vertebrate hypoxia adaptation.
Results
Muscle mass and phenotype variation in high-altitude songbirds
We found a tendency that body mass increases with the elevation in snow finches and tree sparrow (p value <0.05) (S1 Fig and S1 Data). And high-altitude songbirds have greater pectoralis muscle mass compared to lowlanders (p values <0.05) (Fig 1A). To exclude the influence of body size, we compared their muscle mass index (MMI, calculated as muscle mass/ body length2), which is also greater in highlanders (p values <0.05) (Fig 1B). Therefore, for highland songbirds, skeletal muscle is an important contributor to the increase in body weight.
Considering muscle fiber-type transition when ascending to high altitudes in mammals [17] and a greater numerical density of oxidative fibres in high-altitude deer mice [18], we further compared muscle phenotypic variations among altitudinal species and populations. Histological examination revealed that the pectoralis muscles of all birds contained only fast oxidative glycolytic fibers. Despite the homogeneous fiber-type composition, highland tree sparrows and snow finches respectively increased fiber size and myofibril size compared to lowlanders (p values <0.05) (Fig 1C). Previous studies found two opposite phenotypes that lower or higher mitochondrial volume densities occurred in human or high-altitude deer mice [19–21]. We found that highland birds have higher total mitochondrial volume density and proportion of subsarcolemmal mitochondria in the pectoralis (p values <0.05) (Figs 1D and S2A). Overall, highland songbirds evolved substantial muscle phenotypes differing from these of lowland birds.
High-altitude songbirds have improved insulin sensitivity
Skeletal muscle is the major tissue involved in the control of glucose homeostasis and a major determinant of resting energy expenditure [7, 22]. We observed that fasting plasma glucose was lower in highlanders compared to lowlanders, whereas plasma insulin was the same (Figs S2B and 1E). Furthermore, glucose and glycogen concentrations significantly increased in highland bird muscle (p values <0.05) (S2C Fig), highlighting a potential improvement in glucose utilization and insulin sensitivity. To further test this hypothesis, we compared the control of glucose homeostasis using glucose and insulin tolerance tests (GTT and ITT). Highlanders exhibited a more rapid normalization of blood glucose during GTT and also had a greater decrease in blood glucose concentration at 30 min and continued to 60 min after insulin injection (p values <0.05) (Fig 1F), indicating increased insulin sensitivity is an adaptive strategy in highlanders. To determine whether improved insulin sensitivity is associated with changes in energy expenditure, we conducted respirometry to measure metabolic features at rest. Highland songbirds showed a significant increase in respiratory exchange ratio (RER) (p values <0.05), but lower O2 consumption compared to lowlanders (ANCOVA with body mass as a covariate; F2,17 = 5.055, p value = 0.019) (Fig 1G and S1 and S2 Tables).
Additionally, discordant to the report that exercise improves glucose uptake through inducing expression of Glut4 and its translocation [23], lower relative activity was observed in highlanders compared to lowlanders (p values <0.05) (S2D Fig), indicating that songbirds keep quiet to reduce hypoxia stress and insulin sensitivity does not result from physical exercise at high altitudes. Taken together, improvement in insulin sensitivity is probably related with glucose utilization and muscle mass increase, as well as the decrease of oxygen consumption at the resting state.
Molecular regulation of muscle physiology under hypoxia
To extract relevant information on the muscle phenotype, we then performed weighted gene correlation network analysis (WGCNA) to generate modules of genes correlated with muscle phenotypes using all expressed genes. Five, six and seven modules were identified to have a significant correlation with the first principal components of the PCA on the basis of 12 phenotypic traits at the Pa.mo (L)—Pa.mo (H), Pa.mo (L)—Py.ru and Pa.mo (L)—On.ta levels, respectively (correlation coefficient> 0.65, p values< 0.05) (S3A, S3B and S3C Fig). MEF2C was an up-regulated gene in all highland birds and identified as a shared hub gene in both the Pa.mo(L)—Pa.mo (H) and Pa.mo (L)—Py.ru networks of muscle phenotype-related modules (correlation coefficient> 0.65, p values< 0.05) (S3A and S3B Fig and S2 and S3 Datasets). We also found that MEF2C was a hub gene in yellow module of Pa.mo (L)—On.ta network, although correlation coefficient of this module was 0.50 (p value = 0.009) (S3C Fig and S3 Data). Genes of these modules containing MEF2C were related with muscle structure development (p value< 0.05), blood vessel development or remodeling (p values< 0.05) and insulin signaling pathway or negative regulation of insulin secretion (p values< 0.05) (S3A, S3B and S3C Fig). In addition, EPAS1 played an important role in the Pa.mo (L)—Py.ru and Pa.mo (L)—On.ta turquoise modules and was also an up-regulated gene involving in hypoxic response (S3B and S3C Fig and S3 Data). The expression of MEF2C and EPAS1 was significantly correlated with PC1 of muscle phenotypes (S3D Fig). Protein levels of these two genes in pectoral major muscle increased in snow finches; however, in highland tree sparrow, only MEF2C was up-regulated (Fig 2A). The expression of MEF2C showed a significantly positive association with muscle mass across species (R = 0.75, p value <0.0001) (Fig 2B). These results indicate that MEF2C and EPAS1 could be crucial to maintain muscle mass and insulin response under hypoxia.
We next tested the function of MEF2C and EPAS1 using RNA interference in the muscle when songbirds were exposed to low oxygen tension (9% O2). To avoid hypobaropathy resulting from a rapid drop in oxygen tension, we gradually decreased oxygen pressure by 2% per day from 20.9% to 9% and exposed siRNA-treated birds to 9% O2 for 15 days (Fig 2C). We found that all experimental birds underwent a reduction in body weight, but MEF2C and EPAS1 double knockdown caused greater weight loss (S4A and S4B Fig), which resulted from a significant decrease in muscle mass (p values <0.05) (S4C Fig). This indicates that MEF2C rather than EPAS1 maintains body weight through repression of muscle loss under chronic hypoxia. Protein expression levels of MEF2C and EPAS1 markedly decreased after 15 days of hypoxia in the target siRNA groups (Fig 2D). Histological examination showed EPAS1 knockdown exhibited significantly large muscle fibers (p values <0.05) (Fig 2E). A larger myofibril diameter occurred in control birds with both high expression of MEF2C and EPAS1(p values <0.05) (Fig 2F), which probably suppresses protein breakdown under hypoxia. In addition, EPAS1 knockdown showed a lower capillary density (S4D Fig). We also observed severe mitophagy but no decrease of mitochondrial volume density when MEF2C and EPAS1 were simultaneously interfered (p values <0.05) (Figs 2F and S4E). Overall, MEF2C and EPAS1 can control muscle fiber or myofibril disintegration, capillarity, and mitophagy in hypoxia.
We further investigated whether MEF2C and EPAS1 contributed to whole body glucose homeostasis. Birds with low MEF2C expression showed an insulin resistance with significantly higher blood glucose levels after 3 h of fasting and had similar insulin concentration. Additionally, these birds cleared glycemia more slowly in 45 min (also 60 min in hypoxia control group) during GTT and had a slow glucose clearance after 15 min during ITT than the other groups (p values <0.05) (Figs 2G and S4F). Furthermore, birds also exhibited reduced glycogen deposits and showed more relative activity but less RER at both resting and movement states (p values <0.05), which were observed with increasing oxygen consumption during the resting state (ANCOVA with body mass as a covariate; F2,10 = 7.158, p value = 0.012) (Figs 2H and S4G and S3 and S4 Tables). These results suggest that MEF2C might sustain insulin sensitivity, glucose homeostasis and low physical activity.
Oxygen response in high- and low-altitude songbirds
To exclude the effects of plasticity and temperature on the high-altitude songbirds, we tried to perform hypoxia (9%) and normoxia (20.9%) treatments at 25°C for highland tree sparrows and rufous-necked snow finch (natives from 4200m and 4400m, respectively). The deaths usually occurred during the first week for all wild-captured songbirds but snow finches also had a second outbreak of death during the last week because of the persistent losses of body weight (S5A and S5B Fig). High-altitude tree sparrows were slightly decreased body weights and muscle masses in normoxia but no significant difference compared to hypoxia treatment (Figs 3A and S5C). Flight muscle mass, muscle fiber size and capillary per fiber of high-altitude tree sparrows were increased in normoxia and hypoxia compared to low-altitude counterparts (p values < 0.05; Figs 3B, S5C and S5D) despite body weight in normoxia closed to these of low-altitude tree sparrows (S5A Fig), showing that these were probably not plastic in adult high-altitude tree sparrows. Mitochondrial volume density of high-altitude tree sparrows was similar to this of lowlanders in normoxia, suggesting a plasticity (Fig 3C), but high-altitude tree sparrows in hypoxia did not have an obvious mitophagy compared to in normoxia and lowland tree sparrow (S5E Fig). Additionally, hypoxia slightly induced but did not decrease the protein level of MEF2C in high-altitude tree sparrows which was still increased MEF2C abundance compared to lowlanders both in hypoxia and normoxia(p values <0.05, Fig 3D). High-altitude birds showed lower glycemia after three hours fasting (p values <0.05; S5F Fig), a slight lower resting oxygen consumption (Fig 3E) and a significant higher RER and exercise RER in hypoxia compared to in normoxia (p values <0.05; Fig 3F). Highland tree sparrow in hypoxia but not normoxia decreased oxygen consumption compared to lowland tree sparrows both in normoxia and hypoxia (ANCOVA with body mass as a covariate; F2,11 = 9.958, p value = 0.003) (Fig 3E and S3 and S5 Tables). We also found that insulin sensitivity was increased in 30 min and 45 min (Fig 3G) but similar insulin level was in hypoxia compared to normoxia (S5G Fig), which likely drove glucose utilization and glycogen storage in high-altitude tree sparrows (p value <0.01; S5H Fig). Overall, hypoxia failed to induce muscle loss in highland adult tree sparrows but improved insulin sensitivity to sustain glucose metabolism.
To examine the effect of hypoxia on the lowland tree sparrows, we compared the changes in muscle physiology of lowland tree sparrows under hypoxic and normoxic conditions. Body weight significantly decreased on the third day and this low body weight was sustained until the 15th day under 9% oxygen conditions (p values <0.05) (Fig 3A). Hypoxia induced a significant decrease in muscle mass (S5C Fig). Hypoxic exposure causes mitochondrial autophagy at the cellular level [24]. Surprisingly, we observed an increase in mitochondrial volume density although there was an obvious mitophagy in hypoxic birds (p values <0.05) (Fig 3C), which was probably caused by an increase in metabolic rate (ANCOVA with body mass as a covariate; F1,7 = 5.740, p value = 0.048) (Fig 3E and S3 and S6 Tables). Hypoxic birds had lower glycemia after three hours fasting (S5F Fig) and higher RER of exercise, suggesting the increase of glucose utilization (p values <0.05) (Fig 3F). However, we found glucose intolerance during GTT but failed to find the change on insulin sensitivity during ITT in hypoxia birds (Fig 3G). High insulin content in hypoxia probably contributed to the increase of glucose utilization and glycogen storage (p values <0.05) (S5G and S5H Fig). From our experiment, hypoxia slightly decreased the protein abundance of MEF2C (Fig 3D). These results suggest that lowland native exposure to chronic hypoxia undergo muscle loss and increases energy expenditure rather than the similar muscle physiology of highlander, indicating a different genetic background between high- and low-altitude populations.
Genetic basis of muscle physiology in high-altitude songbirds
Next, we identified genes with high amounts of genetic differentiation between altitudinal tree sparrows. We performed a principal component (PC) analysis that suggested that individuals from the two populations could not be distinguished on PC1 of CDS region SNPs, with obvious overlap between high- and low-altitude populations (Fig 4A). Moreover, the weighted Fst between the two populations was 0.014 (mean Fst = 0.008), indicating that the overall genetic differentiation was extremely low. In the RNA-seq data, we identified 112 SNPs in 97 genes with significantly different allele frequencies (p value <0.001), and 23 genes with high fixation index (CDS region Fst >0.20) between high and low altitudinal populations (Fig 4B and 4C and Table 1). PC1 of a principle component analysis of SNPs in these candidate genes explained over 20% of the observed variations in the data and could distinguish the two altitudinal populations (Fig 4E).
Remarkably, among these 23 candidates were four genes with highly differentiated SNPs and with nonsynonymous changes which of three genes were convergent with high-altitude endemic species and were good candidates for influencing hypoxia response (Table 1). One of these genes was hemoglobin subunit beta (HBB), of which two substitutions (Ser80Gly and Val112Ile) were convergent with the Py. ru. The Ser80Gly mutant in Andean high-altitude house wren (Troglodytes aedon) and Val112Ile mutant in band-winged nightjar (Hydropsalis longirostris) are thought to increase Hb–O2 affinity (S6A and S6B Fig) [25, 26]. One candidate gene was methyltransferase-like 8 (METTL8), which has been associated with metabolic disease in humans [27], body weight in chicken [28] and myogenic differentiation by influencing myogenic gene expression (SRF and MEF2) in mouse [29]. The other candidate gene, BCL2 interacting protein 3 like (BNIP3L) which was crucial for hypoxia-induced mitophagy and tumor progression [30–32], had two convergent missense mutations with the Py. ru. but only one with On. ta.
We used an association test to investigate whether any of the SNPs with the most significant allele frequency differences between the populations are associated with body and muscle mass. Notably, two SNPs in the METTL8 gene encoding two missense mutations occurred in all high-altitude birds and were significantly related with muscle mass and body weight (Figs 4D, S6C–S6F and Table 1). The METTL8 gene showed a G-to-A transition at site 326 and 395, causing cysteine-to- tyrosine and serine-to- asparagine nonsynonymous substitutions, respectively, with a higher A frequency in high-altitude songbirds (Fig 4G and 4F). Homozygous individuals for the highland genotype had greater body weight and muscle mass than both heterozygous and homozygous individuals for lowland genotype (p values <0.05) (Figs 4H and S6G). MEF2C expression was also markedly upregulated in the homozygous highland genotype (p values <0.05) (Fig 4G). These results suggest that the convergence in METTL8 may contribute to insulin sensitivity and further muscle mass increase via activating MEF2C in high-altitude birds (Fig 4I).
Discussion
In this study, we revealed convergent physiological and genetic mechanisms of the skeletal muscle resistance to hypoxia in high-altitude songbirds. High-altitude birds suffer from severely limited oxygen supplement due to decrease in the partial pressure of oxygen. Therefore, the physiological responses to highland hypoxia are generally found by altering oxygen utilization in pectoralis muscle [9]. We found an increase of muscle mass resulting from improvement in insulin sensitivity in high-altitude songbirds. Muscle mass and phenotype variations in highland songbirds are probably associated with the flight performance in hypobaric environment and thermogenesis in cold temperature [17, 33]. Animals usually increase their resistance to environmental challenges by altering their behavior responses, such as torpor and hibernation in birds and rodents in response to cold stress [34]. Therefore, reduced physical activity can decrease energy expenditure in high-altitude hypoxia. Additionally, it is confirmed an increase in carbohydrate consumption and a lower plasma glucose in high-altitude animals [10, 35], whilst how birds regulate blood glucose remains unclear. Here we present the first compelling evidence that improving insulin sensitivity may play an important role on carbohydrate utilization and low physical activity in high-altitude songbirds.
Comparative genome and transcriptome analyses across altitudinal species have revealed that convergent evolutionary changes in expression and sequence levels are mainly related to energy metabolism [12, 14]. However, only a few studies have confirmed that these convergent signatures do indeed improve individual fitness to high-altitude environments and physiological evidence are currently unavailable. In this study, combined transcriptomic and physiological analyses demonstrated a convergence on muscle physiology and energy metabolism in high-altitude songbirds. In addition, we verified using RNA interference that convergent MEF2C signaling was pivotal to maintain muscle mass and glucose homeostasis and thus decreased physical activity under hypoxia. However, EPAS1 controls mitophagy, protein breakdown and capillarity, which suggests the evolutionary conservation of function in birds and mammals [36, 37]. MEF2C is a key transcriptional regulator that drives skeletal muscle development during embryogenesis, and hippocampal-dependent memory [38, 39]. Our findings are discordant with these observations, as we highlight a novel function for MEF2C in sustaining postnatal muscle mass and glucose metabolism under hypoxia in songbirds.
Our previous study revealed that high-altitude tree sparrows colonized the QTP for millennia and underwent weak genomic differentiation with the lowland population, whereas polygenic adaptation rather than phenotypic plasticity promoted muscle phenotypic adaptation to high-elevations [26]. Consistent with this result, low-altitude tree sparrows suffered from muscle loss and increased energy expenditures under hypoxia, which was conserved in many mammals and even human [3–5], but converse with high-altitude songbirds. In this study, we identified multiple candidate genes for muscle physiological adaptation to hypoxia in highland tree sparrows including HBB which was the first most highly differentiated gene between the populations and might increase Hb-O2 affinity at two nonsynonymous substitutions [25, 26]. In addition, lowland birds had an increase in capacity of mitochondrial autophagy and muscle loss under chronic hypoxia, suggesting that convergent sites in BNIP3L of highland songbirds probably prevented mitophagy and apoptosis under hypoxia as reported by previous studies [30–32]. We further investigated one of the candidate genes, METTL8, in detail. Two missense mutations in METTL8 with high genetic differentiation between altitudinal populations was significantly associated with body weight, muscle mass and MEF2C expression. These two mutations in high-altitude tree sparrows also occurred in other two endemic species of the QTP, suggesting that convergent substitutions of this gene may come from standing variations in high-altitude songbirds. METTL8 is identified as a new mRNA m3C writer enzyme, and plays a crucial role in myogenic differentiation and metabolic disease [27, 29]. Knockdown of METTL8 decreased MEF2 expression in embryonic mesenchymal cells [29]. These results highlight that high-altitude genotype of METTL8 could cause unique muscle physiological resistance to hypoxia by regulating MEF2C transcript abundance, and rapid evolution of muscle hypoxia tolerance can be explained by large-effect alleles rather than several loci of small effect.
Methods
Ethics statement
All procedures performed on animal were approved by the Animal Care Committee of the Institute of Zoology, Chinese Academy of Sciences. The approval number was IOZ20160049.
Data collection and sampling
We measured 467 museum specimens belonging to 3 Passeridae species from the National Zoological Museum, Institute of Zoology, Chinese Academy of Sciences (see S1 Data). Using mist-net in summer of 2016, highland birds were caught from Qinghai Province at about 3,200 m (10 tree sparrows) and 3,900 m (10 rufous-necked snow finch and 10 white-rumped snowfinch), and lowland birds (n = 8) of the tree sparrow were from Yanqi Lake of Beijing at 80 m. Blood samples were collected and centrifuged at 2,500 rpm for 20 min after fasting for two hours. Pectoral major muscle was dissected immediately following euthanasia. Plasma and tissues were flash-frozen in liquid nitrogen and stored at -80°C until analysis.
Behavioral measures
In the field, we mainly conducted direct observation on birds and combined video record on captive birds. Briefly, subjects were observed continuously for a period of 2 h and behaviors were recorded. In the laboratory, general activity patterns were measured using video recordings, and each treatment was recoded for at least four hours in daylight. The behaviors were divided into three parts, namely, eating, moving and sedentary. Then we calculated the time when birds had been actively eating, moving and inactive (sedentary).
Histology and transmission electron microscopy
Oxidative muscle type and capillarity were analyzed as previously described [9]. Briefly, the pectoral major muscle was dissected and samples were taken third way along the sternum, covered with OCT, and frozen in isopentane (cooled in liquid N2). Muscle sections (10 microns) were obtained in a Cryostat Microtome (Leica CM900, Germany) at -20°C. Succinate dehydrogenase and alkaline phosphatase staining were performed to identify oxidative muscle fibers and muscle capillaries, respectively. The sections were imaged using 10X objective of light microscopy (Leicas DM750) and at least eight images were analyzed for each sample using image J software.
Flight muscles were removed from an intermediate depth and then fixed at 4°C for 24–48 h in 2% glutaraldehyde in 0.1M PBS buffer at pH 7.4. Small muscle blocks (2 mm × 2 mm) were prepared and post-fixed in 1% osmium tetroxide in 0.1 M PBS for 1 h, dehydrated in a graded ethanol series (50%, 70%, 70%, 95%, 95%, 100%, 100%), and embedded in epoxy resin. Ultra-thin sections were cut on a Leica UC7 ultramicrotome and placed on copper grids. The sections were post-stained with uranyl acetate and lead citrate. Images were collected using a transmission electron microscope (Tecnai G2 F20 TWIN TMP, USA). We measured mitochondrial volume density using stereological methods as previously described [40, 41]. Grid size of 90 nm was used and mitochondrial volume was estimated through number of intersection in a mitochondrion dividing total number of intersection at a square size of 4,460 × 4,460 nm. Fracted mitochondria were considered to be autophagy and the ratio of number of fracted mitochondria to normal mitochondria were calculated. We identified mitochondria nearby myolemma and not surrounded by myofibril as subsarcolemmal mitochondria.
Plasma insulin and tissue glycogen measurement
Fasting plasma was obtained from each species. Insulin concentrations were detected using insulin ELISA kit (Cloud-Clone), and muscle glycogen were determined with a relative assay kit (Solarbio) according to the manufacturer’s instructions.
Glucose quantification
Flight muscle (n = 10 each) and plasma (n = 7 each) metabolites were extracted by methanol/chloroform protocol, as previously described [42]. Approximately 1 mL of a chloroform/methanol/distilled water mixture (1:4:1) was added to Eppendorf tubes containing ∼50 mg flight muscle powdered in lipid nitrogen. Then, 450 μL of methanol/distilled water mixture (8:1) was added to Eppendorf tubes containing ∼50 μL plasma. 20 μg, and 10 μg of heptadecanoic acid as well as decanoic acid as double internal standard were added into the muscle and plasma tube, separately. After vortexing for 10 s, the samples were placed on ice for 15 min and then in a sonication bath for 15 min. After centrifugation at 12,000 rpm for 15 min at 4°C, 200 μL of the supernatant were transferred to a 2 mL auto-sampler vial and dried in a vacuum oven. Dried samples were derivatized using methoxyamine hydrochloride solution (50 μL of 15 mg·mL−1 methoxyamine in pyridine). The mixture was kept for 16 h at room temperature for methoxymation, and then 50 μL of TMSFA containing 1% TMCS were added for trimethylsilylated. After 1 h trimethysilylation, 30 μL of hexane was added and then transferred to an insert in a 2-mL autosampler vial for GC-MS analysis. A 2 μL of the derivatized sample was injected by an Agilent 4683B GC auto-sampler (Agilent Technologies, Atlanta, GA, USA) into an Agilent 6890N gas chromatograph with 5973 mass spectrometry at 250°C without splitting. The gas chromatography-mass spectrometry analysis program was as follows: the GC oven temperature was set at 90°C for 1 min, and then increased to 175°C at 5°C min−1, held for 3 min, increased to 270°C at 3°C min−1, 310°C at 20°C min−1, and finally held for 15 min. Identification of glucose from GC-MS analysis was performed by using glucose (Sigma-Aldrich) as standard.
Respirometry
We measured the rates of oxygen consumption VO2 by open flow-through respirometry (Foxbox). For resting metabolic rate (RMR), birds were maintained in a 1.0L transparent metabolic chamber with an air flow rate of 500 mL/minute and O2 and CO2 concentration within the chamber were recorded for over 15 min. The lowest 1-min average of oxygen consumption was taken as the RMR. For exercise oxygen consumption, we used a metabolic flight wheel as previously described [43]. An incurrent flow rate of 5 L/min and a subsampling rate of 100 mL/min were used. We exercised the birds using an initial speed of wheel rotation (approximately 0.3 m/s) and increased the rotation speed by approximately 0.1 m/s every 2 min until birds could no longer keep position and VO2 no longer increased with increasing speed. We took the mean of two minutes of peak O2 consumption and CO2 production to calculate exercise RER. The tests lasted for 20 min. We corrected for the effects of body mass on RMR and exercise oxygen consumption using an analysis of covariance (ANCOVA) with body mass as covariate. Statistical tests were performed using IBM SPSS 22 (SPSS Inc.) with a significance level set to p value = 0.05.
RNA sequencing and differential expression analysis
Total RNA was extracted from each of the 38 sample (8 lowland tree sparrows, 10 highland tree sparrows, 10 rufous-necked snow finches, and 10 white-rumped snow finches) using Trizol RNA isolation reagents (Invitrogen Corp., Carlsbad, CA). RNA integrity was assessed using the RNA Nano 6000 Assay Kit of the Agilent Bioanalyzer 2100 system (Agilent Technologies, CA, USA). A total amount of 3 μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using NEBNext Ultra RNA Library Prep Kit for Illumina(NEB, USA) following manufacturer’s recommendations and index codes were added to attribute sequences to each sample.
Trimmomatic [44] was used to filter reads containing adapter, reads containing ploy-N and low quality reads based on read quality checked with FASTQC [45]. The parameters used were as follows: sliding window = 4-bp; Phred33 quality scores = 20; min read length = 50. Adapter sequences, when detected, were removed. Clean data with high quality were mapped from each species to respective genomes [16] using STAR with default parameters [46]. We used the reciprocal best-hit method to generate tree sparrow-rufous-necked snowfinch and tree sparrow- white-rumped snowfinch orthologs, respectively. The orthologs shared by three species were obtained by intersecting the lists of the above two orthologs. After the reads were mapped to the reference genomes, expression quantification of the genes and transcripts were performed using RSEM [47]. Expression levels for genes with one-to-one orthologs in all three bird species (n = 12,951) were normalized with a RLE (relative log expression) method across muscle samples [48, 49].
To identify genes related with flight muscle variation across altitudinal songbirds, we performed differential expression analysis. Differentially expressed genes (DEGs) were calculated based on the negative binomial distribution and independent filtering was enabled in the program Deseq2 with a false discovery rate (FDR) < 0.05 based on Benjamini-Hochberg method to account for type I error in identifying significantly DEGs [50]. The cut off values for log2-fold change were set at 0.59 and -0.59. Only genes with count number > 10 in at least four samples were included in differential gene expression analysis (see S2 Data).
Weighted gene co-expression network analysis (WGCNA)
We performed WGCNA to identify gene modules associated with muscle phenotypic variance and its potentially genetic basis. The first three principal components of the PCA on 12 phenotypes were used to analysis involved in body weight, lipid droplet volume density, muscle fiber (fiber density, fiber area and myofibril size), capillarity (capillary density, capillary size and capillary per fiber) and mitochondrial phenotypes (distribution of subsarcolemmal mitochondria, volume densities of subsarcolemmal and intermyofibrillar mitochondria, and total mitochondria volume density). All expressed genes were analyzed. Briefly, we constructed the weighted gene co-expression network using the normalized, log2-transformed counts to analyze all expressed genes with the blockwiseModules function in WGCNA [51] for the Pa.mo (L)—Pa.mo (H), Pa.mo (L)—Py.ru and Pa.mo (L)—On.ta levels, respectively. For our analysis, the parameters used were as follows: for Pa.mo (L)—Pa.mo (H) network, maximum block size = 5000 genes, power (β) = 4; for Pa.mo (L)—Py.ru network, maximum block size = 5000 genes, power (β) = 12; for Pa.mo (L)—On.ta network, maximum block size = 5000 genes, power (β) = 18; minimum module size = 25; minimum height for merging modules = 0.25; maximum height for cutting the tree = 0.90. The remaining parameters were kept at the default settings.
Co-expression modules associated with phenotypes were identified using a principal component analysis (PCA) of gene expression with the blockwiseModules function in WGCNA. Each module was summarized by an eigengene, which is the first principal component of the scaled module expression. Trait-related modules were decided by correlation coefficient> 0.65 and p values< 0.05. To specify genes potentially explaining muscle phenotypic variation, we identified hub genes which may be central to the architecture of the regulatory networks represented by each co-expression module. The hub genes were calculated by their first principal component (PC1), the module eigengene (a summary of overall module expression) [52]. Genes which had higher degree (top 25% of genes in module) were identified as hub genes in trait-related module network. We then implemented gene ontology categories within modules positively and negatively correlated with muscle traits using G:PROFILER (FDR < 0.05) [53]. According to the results of functional GO enrichment analysis, intramodular hub genes were identified as candidate genes potentially driving muscle variation.
Genetic variation analysis
Clean data of RNA sequence from high and low altitude populations (12 vs. 10 individuals) was aligned to tree sparrow CDS reference using BWA [54] with default parameters. Only reads that were mapped uniquely to the reference were retained. After high-quality filtering with custom scripts, SNPs were called on all 22 individuals with Samtools. Local realignment over indel positions were performed using GATK [55]. SNPs with missing data in two or more individuals in altitudinal populations were excluded. Allele frequencies and fixed index (Fst) were estimated using Vcftools [56]. Fisher’s exact test was used to determine the statistical significance of the allele frequency differences between highland and lowland populations at exonic SNP positions. Genotype–phenotype correlations between about 120000 SNPs and muscle mass as well as body weight was performed using Plink though the default settings. We further verified the association from results of Plink using a linear mixed model as implemented in the R package lme4 [57]. Candidate genes and SNPs were identified according to Fst >0.3, Fisher’s exact test P-value <0.001 and P-value transformation of genotype–phenotype association >4.
Western blot
Western blots were performed on samples of pectoral major muscle, as previously described [58]. 50 mg of muscle were grinded with liquid nitrogen and then the powder was lysed in 1 ml lysis buffer, 50 mM Tris HCl (pH 7.5), 150 mM NaCl, 5 mM EDTA, 0.5% Nonidet P-40, and completeMini protease inhibitor cocktail (BC3640; Solarbio). The samples were placed on ice for 4 min and then centrifugated at 12,000 rpm for 15 min at 4°C. 350 μL of the supernatant were transferred to a 2 mL auto-sampler vial and were diluted by 3 fold of loading buffer. Protein concentration was quantified using the BCA Protein Assay Kit (PC0020; Solarbio). Supernatant were then separated using SDS-PAGE for 1 h at 140V and transferred to PVDF for 2 h at 90V. Nonspecific antibody bind in was blocked in 5% nonfat dry milk powder in 1 Tris-buffered saline for 1 h. Blots were incubated in 1% nonfat dry milk powder in Tris-buffered saline with 0.05% Tween 20 overnight at 4 C with the following antibodies: rabbit anti-MEF2C (bs-4130R; 1:1,000; Bioss), rabbit anti-EPAS1 (bs-1447R; 1:1,000; Bioss), and mouse anti-β-actin (bsm-33036M; 1:1,000; Bioss). Blots were incubated with goat anti-rabbit (bs-0295G; Bioss) or goat anti-mouse antibodies (bs-0296G; 1:3,000; Bioss) for 1 h at room temperature. Blots were detected using protein detection reagents (AY0371; Solarbio). Membranes were quantified using the ImageJ software and expression was quantified.
In vivo experiments and RNA interference
The experimental birds were transported to the laboratory at the Institute of Zoology, Chinese Academy of Sciences after capture. Birds were weighed and housed individually in screen cages to avoid bird injury (35 cm × 35 cm × 35 cm) at a constant temperature (25±1°C) and 12 h:12 h L:D photoperiod. After acclimation to laboratory for about 14 days, 15 individuals with similar weights were randomly separated into three groups for hypoxia experiment. For hypoxia trails, a hypoxic chamber (FLYDWC-50; Fenglei Co., Ltd, China) was used to control oxygen concentration in normobaric state (20.9kPa). After gradual decrease of oxygen concentration from normoxia (20.9%) to hypoxia (9%), two groups of birds (n = 5 each) were injected with 1 nmol of a mixed double-stranded EPAS1 siRNA or MEF2C+EPAS1 siRNAs (siMEF2C: CCAACAAACTGTTCCAGTA; GCCGAACAAATTCAGATAT; GCACTGATATGGACAAAGT; siEPAS1: TCATCCATGTGACCATGAA; GCATGAGTACCATCTTCCA; TGATCTACAATACCCGAAA; Ribobio, Guangzhou, China), respectively. And then these three groups were maintained for 15 days in 9% oxygen tension. Normoxia and hypoxia control groups were injected with equal normal saline. siRNA was vertically injected to 2 millimetres deep of the pectoralis major using microsyringe. The site of injection was closed to bird carina about 0.5 cm. Birds were injected with target siRNA and weighed every three days during cage cleaning for a brief period (<1 h) and were provided with water and mixed seed ad libitum.
High-altitude tree sparrows and rufous-necked snow finch were respectively captured from 4200m and 4400m, and then transplanted to transported to the laboratory at Tibet Agriculture & Animal Husbandry University in Linzhi (elevation 3000m). Birds were kept at a constant temperature (25±1°C) and 12 h:12 h L:D photoperiod for 12 days. After acclimation, 12 tree sparrows and 7 rufous-necked snow finches with similar weights were randomly separated into two groups for hypoxia and normoxia experiment. Oxygen concentration were gradually increased and decreased from 15% to 20.9% and from 15% to 9% (2% per day). And then birds were maintained for 15 days in 20.9% and 9% oxygen tensions.
Glucose and insulin tolerance tests
For the GTT, an oral gavage of 50% D-glucose (3 mg/g body wt; Sigma) was performed after 16 h fasting. Blood glucose was measured at 0, 15, 30, 45, 60, and 120 min via the medial metatarsal vein. For the ITT, human insulin (1.5 mU/g body wt; Novonordisk) was injected to 2 millimetres deep of the enterocoelia using microsyringe after the birds underwent three hours of fasting. Blood glucose was measured at 0, 15, 30, 45, and 60 min via the medial metatarsal vein. Glucose and insulin were measured using an Accu-Check blood glucose meter (Roche Diagnostics) and an insulin ELISA kit (Cloud-Clone), respectively.
Statistical analyses
One-way ANOVA to analyse was used to determine statistical significance across species and treatments. But we also performed two-tailed Student’s t-test if no statistical significance of comparative analyses using above method. For all figures, p value < 0.05 were considered significant. The black and green asterisks respectively present significant differences for One-way ANOVA and two-tailed Student’s t-test. *P < 0.05, **P < 0.01, ***P < 0.001.
Supporting information
S1 Fig [g]
Altitudinal distribution and body variation.
S2 Fig [a]
Muscle phenotypes and metabolic variation.
S3 Fig [a]
Weighted gene correlation network analysis (WGCNA) based on all expressed genes.
S4 Fig [a]
and effect glucose homeostasis and muscle phenotypes.
S5 Fig [a]
Oxygen response in high- and low-altitude tree sparrows.
S6 Fig [a]
Differentiation and regions under selection across altitudinal populations.
S7 Fig [tif]
The full gels for western blots.
S1 Table [doc]
Body mass and RMR for . (L), . (H) and ..
S2 Table [doc]
Influences of body mass (g) and species (. (L), . (H) and .; ANCOVA with mass as a covariate) on RMR (mLO/h).
S3 Table [doc]
Body mass, RMR and EMR and for different treatments.
S4 Table [doc]
Influences of body mass (g) and RNAi treatment (siMEF2C-EPAS1, siEPAS1 and hypoxia control; ANCOVA with mass as a covariate) on RMR (mLO/h).
S5 Table [doc]
Influences of body mass (g) and treatment (normoxia . (L), hypoxia . (L) and hypoxia . (H); ANCOVA with mass as a covariate) on RMR (mLO/h).
S6 Table [doc]
Influences of body mass (g) and treatment (hypoxia . (L) and normoxia . (L); ANCOVA with mass as a covariate) on RMR (mLO/h).
S1 Data [xlsx]
Body mass of three Passeridae species from the National Zoological Museum, Institute of Zoology, Chinese Academy of Sciences.
S2 Data [xlsx]
Differentially expressed genes (DEGs) across altitudinal songbirds.
S3 Data [xlsx]
The hub genes in some trait-related modules at the . (L)—. (H), . (L)—. and . (L)—. levels.
Zdroje
1. West JB. Alexander M. Kellas and the physiological challenge of Mt. Everest. J Appl Physiol. 1987;63(1):3–11. doi: 10.1152/jappl.1987.63.1.3 3305469
2. Semenza GL. Oxygen sensing, hypoxia-inducible factors, and disease pathophysiology. Annu Rev Pathol-Mech. 2014;9:47–71.
3. Tsoli M, Robertson G. Cancer cachexia: malignant inflammation, tumorkines, and metabolic mayhem. Trends Endocrin Met. 2013;24(4):174–83.
4. Fearon KC, Glass DJ, Guttridge DC. Cancer cachexia: mediators, signaling, and metabolic pathways. Cell Metab. 2012;16(2):153–66. doi: 10.1016/j.cmet.2012.06.011 22795476
5. Baracos VE, Martin L, Korc M, Guttridge DC, Fearon KC. Cancer-associated cachexia. Nat Rev Dis Primers. 2018;4:17105. doi: 10.1038/nrdp.2017.105 29345251
6. Scott GR, Elogio TS, Lui MA, Storz JF, Cheviron ZA: Adaptive modifications of muscle phenotype in high-altitude deer mice are associated with evolved changes in gene regulation. Mol Biol Evol 2015; 32(8):1962–1976. doi: 10.1093/molbev/msv076 25851956
7. Amoasii L, Sanchez-Ortiz E, Fujikawa T, Elmquist JK, Bassel-Duby R, Olson EN. NURR1 activation in skeletal muscle controls systemic energy homeostasis. Proc Natl Acad Sci USA. 2019;116(23):11299–308. doi: 10.1073/pnas.1902490116 31110021
8. Leonvelarde F, Sanchez J, Bigard AX, Brunet A, Lesty C, Monge C. High-Altitude Tissue Adaptation in Andean Coots—Capillarity, Fiber Area, Fiber Type and Enzymatic-Activities of Skeletal-Muscle. J Comp Physiol B. 1993; 163:52–58. doi: 10.1007/BF00309665 8459054
9. Scott GR, Egginton S, Richards JG, Milsom WK: Evolution of muscle phenotype for extreme high altitude flight in the bar-headed goose. P R Soc B. 2009; 276(1673):3645–3653. doi: 10.1098/rspb.2009.0947 19640884
10. Schippers M-P, Ramirez O, Arana M, Pinedo-Bernal P, McClelland GB. Increase in carbohydrate utilization in high-altitude Andean mice. Curr Biol. 2012;22(24):2350–4. doi: 10.1016/j.cub.2012.10.043 23219722
11. Wu D-D, Yang C-P, Wang M-S, Dong K-Z, Yan D-W, Hao Z-Q, et al. Convergent genomic signatures of high altitude adaptation among domestic mammals. bioRxiv. 2019:743955.
12. Hao Y, Xiong Y, Cheng Y, Song G, Jia C, Qu Y, et al. Comparative transcriptomics of 3 high-altitude passerine birds and their low-altitude relatives. Proc Natl Acad Sci USA. 2019;116(24):11851–6. doi: 10.1073/pnas.1819657116 31127049
13. Li J-T, Gao Y-D, Xie L, Deng C, Shi P, Guan M-L, et al. Comparative genomic investigation of high-elevation adaptation in ectothermic snakes. Proc Natl Acad Sci USA. 2018;115(33):8406–11. doi: 10.1073/pnas.1805348115 30065117
14. Sun Y-B, Fu T-T, Jin J-Q, Murphy RW, Hillis DM, Zhang Y-P, et al. Species groups distributed across elevational gradients reveal convergent and continuous genetic adaptation to high elevations. Proc Natl Acad Sci USA. 2018;115(45):E10634–E41. doi: 10.1073/pnas.1813593115 30348757
15. del Hoyo J, Elliot A, Sartagal J, Christie D, De Juana E. Handbook of the birds of the world alive. Barcelona: Lynx Edicions. 2014. 2018.
16. Qu Y, Chen C, Xiong Y, She H, Zhang YE, Cheng Y, et al. Rapid phenotypic evolution with shallow genomic differentiation during early stages of high elevation adaptation in Eurasian Tree Sparrows. Natl Sci Rev. 2019.
17. Chaillou T. Skeletal muscle fiber type in hypoxia: adaptation to high-altitude exposure and under conditions of pathological hypoxia. Front Physiol. 2018;9:1450. doi: 10.3389/fphys.2018.01450 30369887
18. Lui MA, Mahalingam S, Patel P, Connaty AD, Ivy CM, Cheviron ZA, et al. High-altitude ancestry and hypoxia acclimation have distinct effects on exercise capacity and muscle phenotype in deer mice. Am J Physiol-Reg I. 2015;308(9):R779–R91. doi: 10.1152/ajpregu.00362.2014 25695288
19. Hoppeler H, Kleinert E, Schlegel C, Claassen H, Howald H, Kayar S, et al. II. Morphological adaptations of human skeletal muscle to chronic hypoxia. Int J Sports Med. 1990;11(S 1):S3–S9.
20. Mahalingam S, McClelland GB, Scott GR. Evolved changes in the intracellular distribution and physiology of muscle mitochondria in high-altitude native deer mice. J Physiol-London. 2017;595(14):4785–801. doi: 10.1113/JP274130 28418073
21. Hoppeler H, Vogt M. Muscle tissue adaptations to hypoxia. J Exp Biol. 2001;204(18):3133–9. 11581327
22. Zurlo F, Larson K, Bogardus C, Ravussin E. Skeletal muscle metabolism is a major determinant of resting energy expenditure. J Clin Invest. 1990;86(5):1423–7. doi: 10.1172/JCI114857 2243122
23. Kennedy JW, Hirshman MF, Gervino EV, Ocel JV, Forse RA, Hoenig SJ, et al. Acute exercise induces GLUT4 translocation in skeletal muscle of normal human subjects and subjects with type 2 diabetes. Diabetes. 1999;48(5):1192–7. doi: 10.2337/diabetes.48.5.1192 10331428
24. Zhang H, Bosch-Marce M, Shimoda LA, Tan YS, Baek JH, Wesley JB, et al. Mitochondrial autophagy is an HIF-1-dependent adaptive metabolic response to hypoxia. J Biol Chem. 2008;283(16):10892–903. doi: 10.1074/jbc.M800102200 18281291
25. Galen SC, Natarajan C, Moriyama H, Weber RE, Fago A, Benham PM, et al. Contribution of a mutational hot spot to hemoglobin adaptation in high-altitude Andean house wrens. Proc Natl Acad Sci USA. 2015;112(45):13958–63. doi: 10.1073/pnas.1507300112 26460028
26. Kumar A, Natarajan C, Moriyama H, Witt CC, Weber RE, Fago A, et al. Stability-mediated epistasis restricts accessible mutational pathways in the functional evolution of avian hemoglobin. Mol Biol Evol. 2017;34(5):1240–51. doi: 10.1093/molbev/msx085 28201714
27. Tobi EW, Slieker RC, Luijk R, Dekkers KF, Stein AD, Xu KM, et al. DNA methylation as a mediator of the association between prenatal adversity and risk factors for metabolic disease in adulthood. Sci Adv. 2018;4(1):eaao4364. doi: 10.1126/sciadv.aao4364 29399631
28. Wang W, Wang J, Zhang T, Wang Y, Zhang Y, Han K. Genome-wide association study of growth traits in Jinghai Yellow chicken hens using SLAF-seq technology. Anim Genet. 2019;50(2):175–6. doi: 10.1111/age.12346 26365057
29. Jakkaraju S, Zhe X, Pan D, Choudhury R, Schuger L. TIPs are tension-responsive proteins involved in myogenic versus adipogenic differentiation. Dev Cell. 2005;9(1):39–49. doi: 10.1016/j.devcel.2005.04.015 15992539
30. Fei P, Wang W, Kim S-h, Wang S, Burns TF, Sax JK, et al. Bnip3L is induced by p53 under hypoxia, and its knockdown promotes tumor growth. Cancer cell. 2004;6(6):597–609. doi: 10.1016/j.ccr.2004.10.012 15607964
31. Gan Z, Fu T, Kelly DP, Vega RB. Skeletal muscle mitochondrial remodeling in exercise and diseases. Cell Res. 2018;28(10):969–80. doi: 10.1038/s41422-018-0078-7 30108290
32. Liu L, Sakakibara K, Chen Q, Okamoto K. Receptor-mediated mitophagy in yeast and mammalian systems. Cell Res. 2014;24(7):787. doi: 10.1038/cr.2014.75 24903109
33. Sun YF, Ren ZP, Wu YF, Lei FM, Dudley R, Li DM: Flying high: limits to flight performance by sparrows on the Qinghai-Tibet Plateau. J Exp Biol. 2016; 219(22):3642–3648. doi: 10.1242/jeb.142216 27609759
34. Ruf T, Geiser F. Daily torpor and hibernation in birds and mammals. Biol Rev. 2015;90(3):891–926. doi: 10.1111/brv.12137 25123049
35. Sweazea KL, Tsosie K, Beckman EJ, Benham PM, Witt CC. Seasonal and elevational variation in glucose and glycogen in two songbird species. Comp Biochem Phys A. 2020:110703. doi: 10.1016/j.cbpa.2020.110703 32283178
36. Scortegagna M, Ding K, Oktay Y, Gaur A, Thurmond F, Yan L-J, et al. Multiple organ pathology, metabolic abnormalities and impaired homeostasis of reactive oxygen species in Epas1−/− mice. Nat Genet. 2003;35(4):331. doi: 10.1038/ng1266 14608355
37. Befani C, Liakos P. The role of hypoxia-inducible factor-2 alpha in angiogenesis. J Cellular Physiol. 2018;233(12):9087–98. doi: 10.1002/jcp.26805 29968905
38. Lin Q, Schwarz J, Bucana C, Olson EN. Control of mouse cardiac morphogenesis and myogenesis by transcription factor MEF2C. Science. 1997;276(5317):1404–7. doi: 10.1126/science.276.5317.1404 9162005
39. Barbosa AC, Kim M-S, Ertunc M, Adachi M, Nelson ED, McAnally J, et al. MEF2C, a transcription factor that facilitates learning and memory by negative regulation of synapse numbers and function. Proc Natl Acad Sci USA. 2008;105(27):9391–6. doi: 10.1073/pnas.0802679105 18599438
40. Weibel E. Profile size and particle size Weibel ER eds. Stereological Methods: Practical Methods for Biological Morphometry, Vol. 1: 51–60. Academic Press, Inc. New York; 1979.
41. Nielsen J, Gejl KD, Hey-Mogensen M, Holmberg HC, Suetta C, Krustrup P, et al. Plasticity in mitochondrial cristae density allows metabolic capacity modulation in human skeletal muscle. J Physiol-London. 2017;595(9):2839–47. doi: 10.1113/JP273040 WOS:000400357800012. 27696420
42. Shi YL, Chi QS, Liu W, Fu HP, Wang DH. Environmental metabolomics reveal geographic variation in aerobic metabolism and metabolic substrates in Mongolian gerbils (Meriones unguiculatus). Comp Biochem Phys D. 2015;14:42–52. doi: 10.1016/j.cbd.2015.03.001 WOS:000355027900005. 25817427
43. Chappell MA, Bech C, Buttemer WA. The relationship of central and peripheral organ masses to aerobic performance variation in house sparrows. J Exp Biol. 1999;202(17):2269–79. 10441080
44. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20. doi: 10.1093/bioinformatics/btu170 WOS:000340049100004. 24695404
45. Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. 2016.
46. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21. doi: 10.1093/bioinformatics/bts635 23104886
47. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC bioinformatics. 2011;12(1):323.
48. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10). ARTN R106 doi: 10.1186/gb-2010-11-10-r106 WOS:000287378900008. 20979621
49. Maza E. In papyro comparison of TMM (edgeR), RLE (DESeq2), and MRN normalization methods for a simple two-conditions-without-replicates RNA-seq experimental design. Front Genet. 2016;7:164. doi: 10.3389/fgene.2016.00164 27695478
50. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12). ARTN 550 doi: 10.1186/s13059-014-0550-8 WOS:000346609500022. 25516281
51. Langfelder P, Zhang B, Horvath S. Defining clusters from a hierarchical cluster tree: the Dynamic Tree Cut package for R. Bioinformatics. 2008;24(5):719–20. doi: 10.1093/bioinformatics/btm563 WOS:000253746400020. 18024473
52. Han JDJ, Bertin N, Hao T, Goldberg DS, Berriz GF, Zhang LV, et al. Evidence for dynamically organized modularity in the yeast protein-protein interaction network. Nature. 2004;430(6995):88–93. doi: 10.1038/nature02555 WOS:000222356800049. 15190252
53. Reimand J, Kull M, Peterson H, Hansen J, Vilo J. g: Profiler—a web-based toolset for functional profiling of gene lists from large-scale experiments. Nucleic Acids Res. 2007;35:W193–W200. doi: 10.1093/nar/gkm226 WOS:000255311500037. 17478515
54. Li H, Durbin R. Fast and accurate long-read alignment with Burrows–Wheeler transform. Bioinformatics. 2010;26(5):589–95. doi: 10.1093/bioinformatics/btp698 20080505
55. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303. doi: 10.1101/gr.107524.110 20644199
56. Danecek P, Auton A, Abecasis G, Albers CA, Banks E, DePristo MA, et al. The variant call format and VCFtools. Bioinformatics. 2011;27(15):2156–8. doi: 10.1093/bioinformatics/btr330 21653522
57. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. arXiv preprint arXiv:14065823. 2014. doi: 10.1016/j.addbeh.2014.05.011 24935795
58. Frolova A, Flessner L, Chi M, Kim ST, Foyouzi-Yousefi N, Moley KH. Facilitative glucose transporter type 1 is differentially regulated by progesterone and estrogen in murine and human endometrial stromal cells. Endocrinology. 2008;150(3):1512–20. doi: 10.1210/en.2008-1081 18948400
Článek vyšel v časopise
PLOS Genetics
2020 Číslo 12
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
- Co obsahují nové pokyny WHO pro dozor nad antibiotickou rezistencí?
- AI místo lékaře? Výzkum zkoumal ochotu Čechů konzultovat své zdraví s umělou inteligencí
- Prof. Jan Škrha: Metformin je bezpečný, ale je třeba jej bezpečně užívat a léčbu kontrolovat
- Ibuprofen jako alternativa antibiotik při léčbě infekcí močových cest
Nejčtenější v tomto čísle
- Exploiting codon usage identifies intensity-specific modifiers of Ras/MAPK signaling in vivo
- Common maternal and fetal genetic variants show expected polygenic effects on risk of small- or large-for-gestational-age (SGA or LGA), except in the smallest 3% of babies
- PEA15 loss of function and defective cerebral development in the domestic cat
- Precision medicine in cats—The right biomedical model may not be the mouse!