Genomic insights into adaptation to high-altitude environments

Chromosome-level genome assembly of Hippophae tibetana provides insights into high-altitude adaptation and flavonoid biosynthesis

BACKGROUND

As an endemic shrub of the Qinghai-Tibetan Plateau (QTP), the distribution of Hippophae tibetana Schlecht. ranges between 2800 and 5200 m above sea level. As the most basal branch of the Hippophae genus, H. tibetana has an extensive evolutionary history. The H. tibetana is a valuable tree for studying the ecological evolution of species under extreme conditions.

RESULTS

Here, we generated a high-quality chromosome-level genome of H. tibetana. The total size of the assembly genome is 917 Mb. The phylogenomic analysis of 1064 single-copy genes showed a divergence between 3.4 and 12.8 Mya for H. tibetana. Multiple gene families associated with DNA repair and disease resistance were significantly expanded in H. tibetana. We also identified many genes related to DNA repair with signs of positive selection. These results showed expansion and positive selection likely play important roles in H. tibetana's adaptation to comprehensive extreme environments in the QTP. A comprehensive genomic and transcriptomic analysis identified 49 genes involved in the flavonoid biosynthesis pathway in H. tibetana. We generated transgenic sea buckthorn hairy root producing high levels of flavonoid.

CONCLUSIONS

Taken together, this H. tibetana high-quality genome provides insights into the plant adaptation mechanisms of plant under extreme environments and lay foundation for the functional genomic research and molecular breeding of H. tibetana.

The adaptation of bumblebees to extremely high elevation associated with their gut microbiota

Bumblebees are among the most abundant and important pollinators for sub-alpine and alpine flowering plant species in the Northern Hemisphere, but little is known about their adaptations to high elevations. In this article, we focused on two bumblebee species, Bombus friseanus and Bombus prshewalskyi, and their respective gut microbiota. The two species, distributed through the Hengduan Mountains of southwestern China, show species replacement at different elevations. We performed genome sequencing based on 20 worker bee samples of each species. Applying evolutionary population genetics and metagenomic approaches, we detected genes under selection and analyzed functional pathways between bumblebees and their gut microbes. We found clear genetic differentiation between the two host species and significant differences in their microbiota. Species replacement occurred in both hosts and their bacteria (Snodgrassella) with an increase in elevation. These extremely high-elevation bumblebees show evidence of positive selection related to diverse biological processes. Positively selected genes involved in host immune systems probably contributed to gut microbiota changes, while the butyrate generated by gut microbiota may influence both host energy metabolism and immune systems. This suggests a close association between the genomes of the host species and their microbiomes based on some degree of natural selection.IMPORTANCETwo closely related and dominant bumblebee species, distributed at different elevations through the Hengduan Mountains of southwestern China, showed a clear genomic signature of adaptation to elevation at the molecular level and significant differences in their respective microbiota. Species replacement occurred in both hosts and their bacteria (Snodgrassella) with an increase in elevation. Bumblebees' adaptations to higher elevations are closely associated with their gut microbiota through two biological processes: energy metabolism and immune response. Information allowing us to understand the adaptive mechanisms of species to extreme conditions is implicit if we are to conserve them as their environments change.

Stronger evidence for relaxed selection than adaptive evolution in high-elevation animal mtDNA

Mitochondrial (mt) genes are the subject of many adaptive hypotheses due to the key role of mitochondria in energy production and metabolism. One widespread adaptive hypothesis is that selection imposed by life at high elevation leads to the rapid fixation of beneficial alleles in mtDNA, reflected in the increased rates of mtDNA evolution documented in many high-elevation species. However, the assumption that fast mtDNA evolution is caused by positive, rather than relaxed purifying selection has rarely been tested. Here, we calculated the dN/dS ratio, a metric of nonsynonymous substitution bias, and explicitly tested for relaxed selection in the mtDNA of over 700 species of terrestrial vertebrates, freshwater fishes, and arthropods, with information on elevation and latitudinal range limits, range sizes, and body sizes. We confirmed that mitochondrial genomes of high-elevation taxa have slightly higher dN/dS ratios compared to low-elevation relatives. High-elevation species tend to have smaller ranges, which predict higher dN/dS ratios and more relaxed selection across species and clades, while absolute elevation and latitude do not predict higher dN/dS. We also find a positive relationship between body mass and dN/dS, supporting a role for small effective population size leading to relaxed selection. We conclude that higher mt dN/dS among high-elevation species is more likely to reflect relaxed selection due to smaller ranges and reduced effective population size than adaptation to the environment. Our results highlight the importance of rigorously testing adaptive stories against non-adaptive alternative hypotheses, especially in mt genomes.

Whole genome analyses reveal weak signatures of population structure and environmentally associated local adaptation in an important North American pollinator, the bumble bee Bombus vosnesenskii

Studies of species that experience environmental heterogeneity across their distributions have become an important tool for understanding mechanisms of adaptation and predicting responses to climate change. We examine population structure, demographic history and environmentally associated genomic variation in Bombus vosnesenskii, a common bumble bee in the western USA, using whole genome resequencing of populations distributed across a broad range of latitudes and elevations. We find that B. vosnesenskii exhibits minimal population structure and weak isolation by distance, confirming results from previous studies using other molecular marker types. Similarly, demographic analyses with Sequentially Markovian Coalescent models suggest that minimal population structure may have persisted since the last interglacial period, with genomes from different parts of the species range showing similar historical effective population size trajectories and relatively small fluctuations through time. Redundancy analysis revealed a small amount of genomic variation explained by bioclimatic variables. Environmental association analysis with latent factor mixed modelling (LFMM2) identified few outlier loci that were sparsely distributed throughout the genome and although a few putative signatures of selective sweeps were identified, none encompassed particularly large numbers of loci. Some outlier loci were in genes with known regulatory relationships, suggesting the possibility of weak selection, although compared with other species examined with similar approaches, evidence for extensive local adaptation signatures in the genome was relatively weak. Overall, results indicate B. vosnesenskii is an example of a generalist with a high degree of flexibility in its environmental requirements that may ultimately benefit the species under periods of climate change.

Insights into the Adaptation to High Altitudes from Transcriptome Profiling: A Case Study of an Endangered Species, Kingdonia uniflora

Kingdonia uniflora is an endangered alpine herb that is distributed along an altitudinal gradient. The unique traits and important phylogenetic position make K. uniflora an ideal model for exploring how endangered plants react to altitude variation. In this study, we sampled nine individuals from three representative locations and adopted RNA-seq technology to sequence 18 tissues, aiming to uncover how K. uniflora responded to different altitudes at the gene expression level. We revealed that genes that responded to light stimuli and circadian rhythm genes were significantly enriched in DEGs in the leaf tissue group, while genes that were related to root development and peroxidase activity or involved in the pathways of cutin, suberin, wax biosynthesis, and monoterpenoid biosynthesis were significantly enriched in DEGs in the flower bud tissue group. All of the above genes may play an important role in the response of K. uniflora to various stresses, such as low temperatures and hypoxia in high-altitude environments. Furthermore, we proved that the discrepancy in gene expression patterns between leaf and flower bud tissues varied along the altitudinal gradient. Overall, our findings provide new insights into the adaptation of endangered species to high-altitude environments and further encourage parallel research to focus on the molecular mechanisms of alpine plant evolution.

Oxygen stress on age-stage, two-sex life tables and transcriptomic response of diamondback moth (Plutella xylostella)

Elucidating the genetic basis of local adaption is one of the important tasks in evolutionary biology. The Qinghai-Tibet Plateau has the highest biodiversity for an extreme environment worldwide, and provides an ideal natural laboratory to study adaptive evolution. The diamondback moth (DBM), Plutella xylostella, is one of the most devastating pests of the global Brassica industry. A highly heterozygous genome of this pest has facilitated its adaptation to a variety of complex environments, and so provides an ideal model to study fast adaptation. We conducted a pilot study combining RNA-seq with an age-stage, two-sex life table to study the effects of oxygen deprivation on DBM. The developmental periods of all instars were significantly shorter in the hypoxic environment. We compared the transcriptomes of DBM from Fuzhou, Fujian (low-altitude) and Lhasa, Tibet (high-altitude) under hypoxia treatment in a hypoxic chamber. Some DEGs are enriched in pathways associated with DNA replication, such as DNA repair, nucleotide excision repair, base excision repair, mismatch repair and homologous recombination. The pathways with significant changes were associated with metabolism process and cell development. Thus, we assumed that insects could adapt to different environments by regulating their metabolism. Our findings indicated that although adaptive mechanisms to hypoxia in different DBM strains could be similar, DBM individuals from Tibet had superior tolerance to hypoxia compared with those of Fuzhou. Local adaptation of the Tibetan colony was assumed to be responsible for this difference. Our research suggests novel mechanisms of insect responses to hypoxia stress.

Deciphering local adaptation of native Indian cattle (Bos indicus) breeds using landscape genomics and in-silico prediction of deleterious SNP effects on protein structure and function

India has 50 registered breeds of native cattle (Bos indicus) which are locally adapted to diverse environmental conditions. This study aimed to investigate the genomic basis of adaptation of native Indian cattle and to predict the impact of key SNPs on the amino acid changes that affect protein function. The Illumina 777 K BovineHD BeadChip was used to genotype 178 native cattle belonging to contrasting landscapes and agro-climatic conditions. The genotype-environment association was investigated with R. SamBada, using 5,74,382 QC passed SNPs and 11 predictor variables (10 multi-collinearity controlled environmental variables and 1 variable as "score of PCA" on ancestry coefficients of individuals). In total, 1,12,780 models were selected as significant (q < 0.05) based on G score. The pathway ontology of the annotated genes revealed many important pathways and genes having a direct and indirect role in cold and hot adaptation. Only ten SNP variants had a SIFT score of < 0.05 (deleterious), and only two of them, each lying in the genes CRYBA1 and USP18, were predicted to be deleterious with high confidence. RaptorX predicted the tertiary structures of proteins encoded by wild and mutant variants of these genes. The quality of the models was determined using Ramachandran plots and RaptorX parameters, indicating that they are accurate. RaptorX and I-Mutant 2.0 softwares revealed significant differences among wild and mutant proteins. Adaptive alleles identified in the present investigation might be responsible for the local adaptation of these cattle breeds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03493-3.

Analysis of serum antioxidant capacity and gut microbiota in calves at different growth stages in Tibet

Introduction

The hypoxic environment at high altitudes poses a major physiological challenge to animals, especially young animals, as it disturbs the redox state and induces intestinal dysbiosis. Information about its effects on Holstein calves is limited.

Methods

Here, serum biochemical indices and next-generation sequencing were used to explore serum antioxidant capacity, fecal fermentation performance, and fecal microbiota in Holstein calves aged 1, 2, 3, 4, 5, and 6 months in Tibet.

Results and Discussion

Serum antioxidant capacity changed with age, with the catalase and malondialdehyde levels significantly decreasing (p < 0.05), and superoxide dismutase levels significantly increasing (p < 0.05) with age. No significant differences (p > 0.05) in total volatile fatty acid levels were noted between the groups. In all groups, Firmicutes, Bacteroidetes, and Actinobacteria were the three most dominant phyla in the gut. Gut microbial alpha diversity significantly increased (p < 0.05) with age. Principal coordinate analysis plot based on Bray-Curtis dissimilarity revealed significant differences (p = 0.001) among the groups. Furthermore, the relative abundance of various genera changed dynamically with age, and the serum antioxidant capacity was associated with certain gut bacteria. The study provides novel insights for feeding Holstein calves in high-altitude regions.

Mimicking Gene-Environment Interaction of Higher Altitude Dwellers by Intermittent Hypoxia Training: COVID-19 Preventive Strategies

Cyclooxygenase 2 (COX2) inhibitors have been demonstrated to protect against hypoxia pathogenesis in several investigations. It has also been utilized as an adjuvant therapy in the treatment of COVID-19. COX inhibitors, which have previously been shown to be effective in treating previous viral and malarial infections are strong candidates for improving the COVID-19 therapeutic doctrine. However, another COX inhibitor, ibuprofen, is linked to an increase in the angiotensin-converting enzyme 2 (ACE2), which could increase virus susceptibility. Hence, inhibiting COX2 via therapeutics might not always be protective and we need to investigate the downstream molecules that may be involved in hypoxia environment adaptation. Research has discovered that people who are accustomed to reduced oxygen levels at altitude may be protected against the harmful effects of COVID-19. It is important to highlight that the study's conclusions only applied to those who regularly lived at high altitudes; they did not apply to those who occasionally moved to higher altitudes but still lived at lower altitudes. COVID-19 appears to be more dangerous to individuals residing at lower altitudes. The downstream molecules in the (COX2) pathway have been shown to adapt in high-altitude dwellers, which may partially explain why these individuals have a lower prevalence of COVID-19 infection. More research is needed, however, to directly address COX2 expression in people living at higher altitudes. It is possible to mimic the gene-environment interaction of higher altitude people by intermittent hypoxia training. COX-2 adaptation resulting from hypoxic exposure at altitude or intermittent hypoxia exercise training (IHT) seems to have an important therapeutic function. Swimming, a type of IHT, was found to lower COX-2 protein production, a pro-inflammatory milieu transcription factor, while increasing the anti-inflammatory microenvironment. Furthermore, Intermittent Hypoxia Preconditioning (IHP) has been demonstrated in numerous clinical investigations to enhance patients' cardiopulmonary function, raise cardiorespiratory fitness, and increase tissues' and organs' tolerance to ischemia. Biochemical activities of IHP have also been reported as a feasible application strategy for IHP for the rehabilitation of COVID-19 patients. In this paper, we aim to highlight some of the most relevant shared genes implicated with COVID-19 pathogenesis and hypoxia. We hypothesize that COVID-19 pathogenesis and hypoxia share a similar mechanism that affects apoptosis, proliferation, the immune system, and metabolism. We also highlight the necessity of studying individuals who live at higher altitudes to emulate their gene-environment interactions and compare the findings with IHT. Finally, we propose COX2 as an upstream target for testing the effectiveness of IHT in preventing or minimizing the effects of COVID-19 and other oxygen-related pathological conditions in the future.

Effects of altitude on meat quality difference and its relationship with HIF-1α during postmortem maturation of beef

This study investigated the differences in meat quality during postmortem aging of yak meat from different altitudes as well as the relationship between the release of hypoxic factor HIF-1α and meat quality. The results showed that the HIF-1α increased with altitude but during aging process, there was an initial increase before a subsequent decrease (p < .05). Moreover, significant increases were showed in glycolytic potential, a* value, pH, HIF-1α mRNA expression, HIF-1α protein expression and shear force with altitude (p < .05). Additionally, the b* value, L* value, water holding power and MFI decreased significantly (p < .05). HIF-1α was shown, by PLS-DA method analysis, to be the main protein marker for differences in the quality during aging time of meat from three altitude groups. HIF-1α protein expression was high correlated with glycolytic potential, pH value, meat color, tenderness and water holding capacity during postmortem aging. The results demonstrated that HIF-1α is a novel marker protein that influences meat quality in yak from different altitudes and that HIF-1α-mediated glycolytic pathway was key to the meat quality during postmortem aging. PRACTICAL APPLICATIONS: Yak meat has the advantages of high protein, low fat, good amino acid and fatty acid composition, so the nutritional value of yak meat is in line with the current best-selling beef with less fat in domestic and foreign markets. But consumers often think that the meat tenderness of yak meat is worse than that of beef and improving the quality of yak meat was worthy of attention specifically. This study investigated the differences in meat quality during postmortem aging of yak meat at different altitudes and the relationship between hypoxic factor HIF-1α release and meat quality.

Dynamic alterations in yak (Bos grunniens) rumen microbiome in response to seasonal variations in diet

Rumen microorganisms play important roles in the healthy growth of yaks. This study investigated changes in yak rumen microbiome during natural grazing at the warm seasons and supplementary feeding at cold seasons. High-throughput sequencing of 16S rRNA and metagenome analysis were conducted to investigate the structures and functions of yak rumen microbial communities. The results indicated that Bacteroidetes and Firmicutes were the most abundant phyla. In addition, Bacteroidetes might play a more important role than Firmicutes during the supplementary feeding stage (spring and winter), but less during natural grazing stage (summer and autumn). KEGG analysis showed that the amino sugar and nucleotide sugar metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, starch and sucrose metabolism, and fructose and mannose metabolism were the main pathways in the microbial community, which were significantly different between seasons. The carbohydrate-active enzymes (CAZyme) annotation revealed that cellulose was an important carbon source for microorganisms in yak rumen. Glycoside hydrolases (GHs) were the most abundant class of CAZymes, followed by glycosyl transferases (GTs), which were important to digestion of oil, cellulose, and hemicellulose in food. These results contribute to the understanding of microbial components and functions in yak rumen.

Glucose and oxygen in the early intrauterine environment and their role in developmental abnormalities

The early life environment can have profound impacts on the developing conceptus in terms of both growth and morphogenesis. These impacts can manifest in a variety of ways, including congenital fetal anomalies, placental dysfunction with subsequent effects on fetal growth, and adverse perinatal outcomes, or via effects on long-term health outcomes that may not be detected until later childhood or adulthood. Two key examples of environmental influences on early development are explored: maternal hyperglycaemia and gestational hypoxia. These are increasingly common pregnancy exposures worldwide, with potentially profound impacts on population health. We explore what is known regarding the mechanisms by which these environmental exposures can impact early intrauterine development and thus result in adverse outcomes in the immediate, short, and long term.

Gene expression underlying parenting and being parented shows limited plasticity in response to different ambient temperatures

Flexible interactions between parents and offspring are essential for buffering families against variable, unpredictable, and challenging environmental conditions. In the subsocial carrion beetle, Nicrophorus orbicollis, mid-summer temperatures impose steep fitness costs on parents and offspring but do not elicit behavioural plasticity in parents. Here, we ask if plasticity of gene expression underpins this behavioural stability or facilitates independent compensation by larvae. To test this, we characterized gene expression of parents and offspring before and during active parenting under benign (20°C) and stressful (24°C) temperatures to identify genes of parents and offspring associated with thermal response, parenting/being parented, and gene expression plasticity associated with behavioural stability of parental care. The main effects of thermal and social condition each shaped patterns of gene expression in females, males, and larvae. In addition, we implicated 79 genes in females as "buffering" parental behaviour across environments. The majority of these underwent significant changes in expression in actively parenting mothers at the benign temperature, but not at the stressful temperature. Our results suggest that neither genetic programmes for parenting nor their effects on offspring gene expression are fundamentally different under stressful conditions, and that behavioural stability is associated primarily with the maintenance of existing genetic programmes rather than replacement or supplementation. Thus, while selection for compensatory gene expression could expand the range of thermal conditions parents will tolerate, without expanding the toolkit of genes involved selection is unlikely to lead to adaptive changes of function.

Comparative genomics reveals putative evidence for high-elevation adaptation in the American pika (Ochotona princeps)

High-elevation environments have lower atmospheric oxygen content, reduced temperatures, and higher levels of UV radiation than found at lower elevations. As such, species living at high elevations must overcome these challenges to survive, grow, and reproduce. American pikas (Ochotona princeps) are alpine lagomorphs that are habitat specialists typically found at elevations >2,000 m. Previous research has shown putative evidence for high-elevation adaptation; however, investigations to date have been limited to a fraction of the genome. Here, we took a comparative genomics approach to identify putative regions under selection using a chromosomal reference genome assembly for the American pika relative to 8 other mammalian species targeted based on phylogenetic relatedness and (dis)similarity in ecology. We first identified orthologous gene groups across species and then extracted groups containing only American pika genes as well as unclustered pika genes to inform functional enrichment analyses; among these, we found 141 enriched terms with many related to hypoxia, metabolism, mitochondrial function/development, and DNA repair. We identified 15 significantly expanded gene families within the American pika across all orthologous gene groups that displayed functionally enriched terms associated with hypoxia adaptation. We further detected 196 positively selected genes, 41 of which have been associated with putative adaptation to hypoxia, cold tolerance, and response to UV following a literature review. In particular, OXNAD1, NRDC, and those genes critical in DNA repair represent important targets for future research to examine their functional implications in the American pika, especially as they may relate to adaptation to rapidly changing environments.

Identifying climate refugia for high-elevation Alpine birds under current climate warming predictions

Identifying climate refugia is key to effective biodiversity conservation under a changing climate, especially for mountain-specialist species adapted to cold conditions and highly threatened by climate warming. We combined species distribution models (SDMs) with climate forecasts to identify climate refugia for high-elevation bird species (Lagopus muta, Anthus spinoletta, Prunella collaris, Montifringilla nivalis) in the European Alps, where the ecological effects of climate changes are particularly evident and predicted to intensify. We considered future (2041-2070) conditions (SSP585 scenario, four climate models) and identified three types of refugia: (1) in-situ refugia potentially suitable under both current and future climate conditions, ex-situ refugia suitable (2) only in the future according to all future conditions, or (3) under at least three out of four future conditions. SDMs were based on a very large, high-resolution occurrence dataset (2901-12,601 independent records for each species) collected by citizen scientists. SDMs were fitted using different algorithms, balancing statistical accuracy, ecological realism and predictive/extrapolation ability. We selected the most reliable ones based on consistency between training and testing data and extrapolation over distant areas. Future predictions revealed that all species (with the partial exception of A. spinoletta) will undergo a range contraction towards higher elevations, losing 17%-59% of their current range (larger losses in L. muta). We identified ~15,000 km² of the Alpine region as in-situ refugia for at least three species, of which 44% are currently designated as protected areas (PAs; 18%-66% among countries). Our findings highlight the usefulness of spatially accurate data collected by citizen scientists, and the importance of model testing by extrapolating over independent areas. Climate refugia, which are only partly included within the current PAs system, should be priority sites for the conservation of Alpine high-elevation species and habitats, where habitat degradation/alteration by human activities should be prevented to ensure future suitability for alpine species.

Cyclophilin A represses reactive oxygen species generation and death of hypoxic non-small-cell lung cancer cells by degrading thioredoxin-interacting protein

Cyclophilin A (cypA) is overexpressed in many types of carcinomas, including non-small-cell lung cancer (NSCLC). However, the effect of anoxia, a critical feature of the carcinoma cell microenvironment, on cypA expression in NSCLC is unknown. Here, formaldehyde-fixed and paraffin-embedded samples were collected from 60 subjects with NSCLC. The protein expression levels of cypA and hypoxia-inducible factor-1α (HIF-1α) were evaluated using immunohistochemistry. Kaplan-Meier analysis showed that subjects with high cypA expression had remarkably shorter progression-free survival than those with low cypA expression. Furthermore, cypA expression levels were significantly related to HIF-1α expression levels (Spearman's correlation=0.34, P<0.0001). To further assess the effect of cypA, an anoxic carcinoma cell model was established. CypA expression was remarkably upregulated in H1299 and A549 cell lines under hypoxic conditions. Overexpression of cypA restored hypoxia-impaired cell growth and prevented reactive oxygen species (ROS) production and cell death in hypoxic A549 and H1299 cells. However, these phenotypes were not altered by the inactive R55A mutant of cypA. Mechanistic studies demonstrated that cypA can bind to and degrade the tumor suppressor protein TXNIP in H1299 and A549 cells. Restored TXNIP expression in cypA-overexpressed and hypoxic NSCLC cells led to increased ROS levels and apoptotic cell numbers and decreased cell growth compared with cypA-overexpressed and hypoxic NSCLC cells. These findings indicate that anoxia results in an increase in cypA expression in NSCLC. Additionally, cypA served as an oncogene during hypoxia by interacting with TXNIP.

Effects of Altitude on the Digestion Performance, Serum Antioxidative Characteristics, Rumen Fermentation Parameters, and Rumen Bacteria of Sanhe Heifers

The production efficiency of dairy cows is affected by altitude, with lower efficiency reported at higher altitudes. However, only a few studies have investigated the digestion performance, serum antioxidative characteristics, rumen fermentation performance, and rumen bacteria of Sanhe heifers at different altitudes. Therefore, in this study, we explored the effects of altitude on these aspects of Sanhe heifers. We evaluated the effects of altitude on the apparent digestibility of nutrients, serum antioxidative characteristics, rumen fermentation parameters, and rumen bacteria in Sanhe heifers. Twenty Sanhe heifers from the same herd and managed with the same practice were used. However, the heifers were from two regions in China: 10 were fed in Hulunbuir City, Inner Mongolia Autonomous Region (119°57'E, 47°17'N; approximately 700 m altitude, named LA) and 10 were fed in Lhasa City, Tibet Autonomous Region (91°06'E, 29°36'N; approximately 3,750 m altitude, named HA). The dry matter intake (DMI), average daily gain (ADG), and DMI/ADG ratio were higher (p < 0.05) in LA than in HA heifers, whereas the apparent total tract digestibility of dry matter, ether extract, and crude proteins were higher (p < 0.05) in the HA group. Compared with LA heifers, the HA heifers showed decreased (p < 0.05) serum concentrations of superoxide dismutase and glutathione peroxidase and increased serum concentration of hydrogen peroxide (p < 0.05). Altitude did not significantly affect the volatile fatty acid concentration in the rumen, but HA presented a lower acetate-to-propionate ratio than LA. The 16S rRNA gene sequencing data showed that altitude significantly affected the rumen microbial composition. At the phylum level, the HA heifers presented a lower relative abundance of Actinobacteria (p < 0.05) and higher relative abundance of Spirochaetae (p < 0.05) than the LA heifers. The correlation analysis revealed that the operational taxonomic units belonging to the genus Prevotella_1 were correlated (p < 0.05) with altitude and DMI. The results indicate that altitude can influence the apparent digestibility of nutrients, serum antioxidant capacity, rumen fermentation, and rumen bacteria composition of Sanhe heifers. The study provides insights into the adaptation mechanism of Sanhe heifers to high-altitude areas.

Comparative analysis of long noncoding RNA and mRNA expression provides insights into adaptation to hypoxia in Tibetan sheep

Tibetan sheep have lived on the Qinghai-Tibetan Plateau for thousands of years and have good adaptability to the hypoxic environment and strong disease resistance. However, the molecular mechanism by which Tibetan sheep adapt to this extreme environment, especially the role of genetic regulation, is still unknown. Emerging evidence suggests that long noncoding RNAs (lncRNAs) participate in the regulation of a diverse range of biological processes. To explore the potential lncRNAs involved in the adaptation to high-altitude hypoxia of Tibetan sheep, we analysed the expression profile of lncRNAs and mRNAs in the liver and lung tissues of sheep using comparative transcriptome analysis between four Tibetan sheep populations (high altitude) and one Hu sheep population (low altitude). The results showed a total of 7848 differentially expressed (DE) lncRNA transcripts, and 22,971 DE mRNA transcripts were detected by pairwise comparison. The expression patterns of selected mRNAs and lncRNAs were validated by qRT-PCR, and the results correlated well with the transcriptome data. Moreover, the functional annotation analysis based on the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases showed that DE mRNAs and the target genes of the lncRNAs were significantly enriched in organ morphogenesis, response to stimulus, haem binding, the immune system, arginine and proline metabolism, and fatty acid biosynthesis. The prediction of mRNA–mRNA and lncRNA–mRNA interaction networks further revealed transcripts potentially involved in adaptation to high-altitude hypoxia, and the hub genes DDX24, PDCD11, EIF4A3, NDUFA11, SART1, PRPF8 and TCONS_00306477, TCONS_00306029, TCONS_00139593, TCONS_00293272, and TCONS_00313398 were selected. Additionally, a set of target genes, PIK3R1, IGF1R, FZD6, IFNB2, ATF3, MB, CYP2B4, PSMD13, and TGFB1, were also identified as candidate genes associated with high-altitude hypoxia adaptation. In conclusion, a collection of novel expressed lncRNAs, a set of target genes and biological pathways known to be relevant for altitude adaptation were identified by comparative transcriptome analysis between Tibetan sheep and Hu sheep. Our results are the first to identify the characterization and expression profile of lncRNAs between Tibetan sheep and Hu sheep and provide insights into the genetic regulation mechanisms by which Tibetan sheep adapt to high-altitude hypoxic environments.

Do habitat and elevation promote hybridization during secondary contact between three genetically distinct groups of warbling vireo (Vireo gilvus)?

Following postglacial expansion, secondary contact can occur between genetically distinct lineages. These genetic lineages may be associated with specific habitat or environmental variables and therefore, their distributions in secondary contact could reflect such conditions within these areas. Here we used mtDNA, microsatellite, and morphological data to study three genetically distinct groups of warbling vireo (Vireo gilvus) and investigate the role that elevation and habitat play in their distributions. We studied two main contact zones and within each contact zone, we examined two separate transects. Across the Great Plains contact zone, we found that hybridization between eastern and western groups occurs along a habitat and elevational gradient, whereas hybridization across the Rocky Mountain contact zone was not as closely associated with habitat or elevation. Hybrids in the Great Plains contact zone were more common in transitional areas between deciduous and mixed-wood forests, and at lower elevations (<1000 m). Hybridization patterns were similar along both Great Plains transects indicating that habitat and elevation play a role in hybridization between distinct eastern and western genetic groups. The observed patterns suggest adaptation to different habitats, perhaps originating during isolation in multiple Pleistocene refugia, is facilitating hybridization in areas where habitat types overlap.

Insights into the plateau adaptation of Salvia castanea by comparative genomic and WGCNA analyses

INTRODUCTION

Salvia castanea, a wild plant species is adapted to extreme Qinghai-Tibetan plateau (QTP) environments. It is also used for medicinal purposes due to high ingredient of tanshinone IIA (T-IIA). Despite its importance to Chinese medicinal industry, the mechanisms associated with secondary metabolites accumulation (i.e. T-IIA and rosmarinic acid (RA)) in this species have not been characterized. Also, the role of special underground tissues in QTP adaptation of S. castanea is still unknown.

OBJECTIVES

We explored the phenomenon of periderm-like structure in underground stem center of S. castanea with an aim to unravel the molecular evolutionary mechanisms of QTP adaptation in this species.

METHODS

Morphologic observation and full-length transcriptome of S. castanea plants were conducted. Comparative genomic analyses of S. castanea with other 14 representative species were used to reveal its phylogenetic position and molecular evolutionary mechanisms. RNA-seq and WGCNA analyses were applied to understand the mechanisms of high accumulations of T-IIA and RA in S. castanea tissues.

RESULTS

Based on anatomical observations, we proposed a "trunk-branches" developmental model to explain periderm-like structure in the center of underground stem of S. castanea. Our study suggested that S. castanea branched off from cultivated Danshen around 16 million years ago. During the evolutionary process, significantly expanded orthologous gene groups, 24 species-specific and 64 positively selected genes contributed to morphogenesis and QTP adaptation in S. castanea. RNA-seq and WGCNA analyses unraveled underlying mechanisms of high accumulations of T-IIA and RA in S. castanea and identified NAC29 and TGA22 as key transcription factors.

CONCLUSION

We proposed a "trunk-branches" developmental model for the underground stem in S. castanea. Adaptations to extreme QTP environment in S. castanea are associated with accumulations of high secondary metabolites in this species.

The relevance of genetic structure in ecotype designation and conservation management

The concept of ecotypes is complex, partly because of its interdisciplinary nature, but the idea is intrinsically valuable for evolutionary biology and applied conservation. The complex nature of ecotypes has spurred some confusion and inconsistencies in the literature, thereby limiting broader theoretical development and practical application. We provide suggestions for how incorporating genetic analyses can ease confusion and help define ecotypes. We approach this by systematically reviewing 112 publications across taxa that simultaneously mention the terms ecotype, conservation and management, to examine the current use of the term in the context of conservation and management. We found that most ecotype studies involve fish, mammals and plants with a focus on habitat use, which at 60% was the most common criterion used for categorization of ecotypes. Only 53% of the studies incorporated genetic analyses, and major discrepancies in available genomic resources among taxa could have contributed to confusion about the role of genetic structure in delineating ecotypes. Our results show that the rapid advances in genetic methods, also for nonmodel organisms, can help clarify the spatiotemporal distribution of adaptive and neutral genetic variation and their relevance to ecotype designations. Genetic analyses can offer empirical support for the ecotype concept and provide a timely measure of evolutionary potential, especially in changing environmental conditions. Genetic variation that is often difficult to detect, including polygenic traits influenced by small contributions from several genes, can be vital for adaptation to rapidly changing environments. Emerging ecotypes may signal speciation in progress, and findings from genome‐enabled organisms can help clarify important selective factors driving ecotype development and persistence, and thereby improve preservation of interspecific genetic diversity. Incorporation of genetic analyses in ecotype studies will help connect evolutionary biology and applied conservation, including that of problematic groups such as natural hybrid organisms and urban or anthropogenic ecotypes.

Disentangling direct and indirect effects of local temperature on abundance of mountain birds and implications for understanding global change impacts

Unravelling the environmental factors driving species distribution and abundance is crucial in ecology and conservation. Both climatic and land cover factors are often used to describe species distribution/abundance, but their interrelations have been scarcely investigated. Climatic factors may indeed affect species both directly and indirectly, e.g., by influencing vegetation structure and composition. We aimed to disentangle the direct and indirect effects (via vegetation) of local temperature on bird abundance across a wide elevational gradient in the European Alps, ranging from montane forests to high-elevation open areas. In 2018, we surveyed birds by using point counts and collected fine-scale land cover and temperature data from 109 sampling points. We used structural equation modelling to estimate direct and indirect effects of local climate on bird abundance. We obtained a sufficient sample for 15 species, characterized by a broad variety of ecological requirements. For all species we found a significant indirect effect of local temperatures via vegetation on bird abundance. Direct effects of temperature were less common and were observed in seven woodland/shrubland species, including only mountain generalists; in these cases, local temperatures showed a positive effect, suggesting that on average our study area is likely colder than the thermal optimum of those species. The generalized occurrence of indirect temperature effects within our species set demonstrates the importance of considering both climate and land cover changes to obtain more reliable predictions of future species distribution/abundance. In fact, many species may be largely tracking suitable habitat rather than thermal niches, especially among homeotherm organisms like birds.

Parallel genomic responses to historical climate change and high elevation in East Asian songbirds

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

RETSAT Mutation Selected for Hypoxia Adaptation Inhibits Tumor Growth

Hypoxia occurs not only in natural environments including high altitude, underground burrows and deep sea, but also in human pathological conditions, such as hypoxic solid tumors. It has been well documented that hypoxia related signaling pathway is associated with a poor clinical outcome. Our group has recently identified multiple novel genes critical for solid tumor growth comparing the genome-wide convergent/parallel sequence evolution of highland mammals. Among them, a single mutation on the retinol saturase gene (RETSAT) containing amino acid switch from glutamine (Q) to arginine (R) at the position 247 was identified. Here, we demonstrate that RETSAT is mostly downregulated in multiple types of human cancers, whose lower expression correlates with worse clinical outcome. We show that higher expression of RETSAT is positively associated with immune infiltration in different human cancers. Furthermore, we identify that the promoter region of RETSAT is highly methylated, which leads to its decreased expressions in tumor tissues comparing to normal tissues. Furthermore, we show that RETSAT knockdown promotes, while its overexpression inhibits, the cell proliferation ability of mouse embryonic fibroblasts (MEFs) and B16 in vitro. In addition, the mice carrying homozygous Q247R mutation (RETSATR/R) is more resistant to xenograft tumor formation, as well as DMBA/TPA induced cutaneous keratinocyte carcinoma formation, compared to littermate wild-type (RETSATQ/Q) mice. Mechanistic study uncovers that the oncogenic factor, the prolyl isomerase (PPIase) Pin1 and its related downstream signaling pathway, were both markedly repressed in the mutant mice compared to the wild-type mice. In summary, these results suggest that interdisciplinary study between evolution and tumor biology can facilitate identification of novel molecular events essential for hypoxic solid tumor growth in the future.

Large- and small-scale geographic structures affecting genetic patterns across populations of an Alpine butterfly

Understanding factors influencing patterns of genetic diversity and the population genetic structure of species is of particular importance in the current era of global climate change and habitat loss. These factors include the evolutionary history of a species as well as heterogeneity in the environment it occupies, which in turn can change across time. Most studies investigating spatio-temporal genetic patterns have focused on patterns across wide geographic areas rather than local variation, but the latter can nevertheless be important particularly in topographically complex areas. Here, we consider these issues in the Sooty Copper butterfly (Lycaena tityrus) from the European Alps, using genome-wide SNPs identified through RADseq. We found strong genetic differentiation within the Alps with four genetic clusters, indicating western, central, and eastern refuges, and a strong reduction of genetic diversity from west to east. This reduction in diversity may suggest that the southwestern refuge was the largest one in comparison to other refuges. Also, the high genetic diversity in the west may result from (a) admixture of different western refuges, (b) more recent demographic changes, or (c) introgression of lowland L. tityrus populations. At small spatial scales, populations were structured by several landscape features and especially by high mountain ridges and large river valleys. We detected 36 outlier loci likely under altitudinal selection, including several loci related to membranes and cellular processes. We suggest that efforts to preserve alpine L. tityrus should focus on the genetically diverse populations in the western Alps, and that the dolomite populations should be treated as genetically distinct management units, since they appear to be currently more threatened than others. This study demonstrates the usefulness of SNP-based approaches for understanding patterns of genetic diversity, gene flow, and selection in a region that is expected to be particularly vulnerable to climate change.

High-altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high-altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high-altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein-coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high-altitude vertebrates than in low-altitude vertebrates. The seven rapidly evolving genes were shared by the high-altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high-altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high-altitude vertebrates was higher than in low-altitude vertebrates as their evolution requires more energy in a high-altitude environment. Our study demonstrates the high-altitude environment (low atmospheric O2 levels) drives vertebrate evolution in mtDNA PCGs.

Chromosome-Level Reference Genome Assembly for the American Pika (Ochotona princeps)

The American pika (Ochotona princeps) is an alpine lagomorph found throughout western North America. Primarily inhabiting talus slopes at higher elevations (>2000 m), American pikas are well adapted to cold, montane environments. Warming climates on both historical and contemporary scales have contributed to population declines in American pikas, positioning them as a focal mammalian species for investigating the ecological effects of climate change. To support and expand ongoing research efforts, here, we present a highly contiguous and annotated reference genome assembly for the American pika (OchPri4.0). This assembly was produced using Dovetail de novo proximity ligation methods and annotated through the NCBI Eukaryotic Genome Annotation pipeline. The resulting assembly was chromosome- scale, with a total length of 2.23 Gb across 9350 scaffolds and a scaffold N50 of 75.8 Mb. The vast majority (>97%) of the total assembly length was found within 36 large scaffolds; 33 of these scaffolds correlated to whole autosomes, while the X chromosome was covered by 3 large scaffolds. Additionally, we identified 17 enriched gene ontology terms among American pika-specific genes putatively related to adaptation to high-elevation environments. This high-quality genome assembly will serve as a springboard for exploring the evolutionary underpinnings of behavioral, ecological, and taxonomic diversification in pikas as well as broader-scale eco-evolutionary questions pertaining to cold-adapted species in general.

Genomic signatures of parallel alpine adaptation in recently evolved flightless insects

Natural selection along elevational gradients has potential to drive predictable adaptations across distinct lineages, but the extent of such repeated evolution remains poorly studied for many widespread alpine taxa. We present parallel genomic analyses of two recently evolved flightless alpine insect lineages to test for molecular signatures of repeated alpine adaptation. Specifically, we compare low-elevation vs. alpine stonefly ecotypes from parallel stream populations in which flightless upland ecotypes have been independently derived. We map 67,922 polymorphic genetic markers, generated across 176 Zelandoperla fenestrata specimens from two independent alpine stream populations in New Zealand's Rock and Pillar Range, to a newly developed plecopteran reference genome. Genome-wide scans revealed 31 regions with outlier single nucleotide polymorphisms (SNPs) differentiating lowland vs. alpine ecotypes in Lug Creek, and 37 regions with outliers differentiating ecotypes in Six Mile Creek. Of these regions, 13% (8/60) yielded outlier SNPs across both within-stream ecotype comparisons, implying comparable genomic shifts contribute to this repeated alpine adaptation. Candidate genes closely linked to repeated outlier regions include several with documented roles in insect wing-development (e.g., dishevelled), suggesting that they may contribute to repeated alpine wing reduction. Additional candidate genes have been shown to influence insect fecundity (e.g., ovo) and lifespan (e.g., Mrp4), implying that they might contribute to life history differentiation between upland and lowland ecotypes. Additional outlier genes have potential roles in the evolution of reproductive isolation among ecotypes (hedgehog and Desaturase 1). These results demonstrate how replicated outlier tests across independent lineages can potentially contribute to the discovery of genes underpinning repeated adaptation.

Transcriptomes of Saussurea (Asteraceae) Provide Insights into High-Altitude Adaptation

Understanding how species adapt to extreme environments is an extension of the main goals of evolutionary biology. While alpine plants are an ideal system for investigating the genetic basis of high-altitude adaptation, genomic resources in these species are still limited. In the present study, we generated reference-level transcriptomic data of five Saussurea species through high-throughput sequencing and de novo assembly. Three of them are located in the highland of the Qinghai-Tibet Plateau (QTP), and the other two are close relatives distributed in the lowland. A series of comparative and evolutionary genomics analyses were conducted to explore the genetic signatures of adaptive evolution to high-altitude environments. Estimation of divergence time using single-copy orthologs revealed that Saussurea species diversified during the Miocene, a period with extensive tectonic movement and climatic fluctuation on the QTP. We characterized gene families specific to the alpine species, including genes involved in oxidoreductase activity, pectin catabolic process, lipid transport, and polysaccharide metabolic process, which may play important roles in defense of hypoxia and freezing temperatures of the QTP. Furthermore, in a phylogenetic context with the branch model, we identified hundreds of genes with signatures of positive selection. These genes are involved in DNA repair, membrane transport, response to UV-B and hypoxia, and reproductive processes, as well as some metabolic processes associated with nutrient intake, potentially responsible for Saussurea adaptation to the harsh environments of high altitude. Overall, our study provides valuable genomic resources for alpine species and gained helpful insights into the genomic basis of plants adapting to extreme environments.

Exposomes to Exosomes: Exosomes as Tools to Study Epigenetic Adaptive Mechanisms in High-Altitude Humans

Humans on earth inhabit a wide range of environmental conditions and some environments are more challenging for human survival than others. However, many living beings, including humans, have developed adaptive mechanisms to live in such inhospitable, harsh environments. Among different difficult environments, high-altitude living is especially demanding because of diminished partial pressure of oxygen and resulting chronic hypobaric hypoxia. This results in poor blood oxygenation and reduces aerobic oxidative respiration in the mitochondria, leading to increased reactive oxygen species generation and activation of hypoxia-inducible gene expression. Genetic mechanisms in the adaptation to high altitude is well-studied, but there are only limited studies regarding the role of epigenetic mechanisms. The purpose of this review is to understand the epigenetic mechanisms behind high-altitude adaptive and maladaptive phenotypes. Hypobaric hypoxia is a form of cellular hypoxia, which is similar to the one suffered by critically-ill hypoxemia patients. Thus, understanding the adaptive epigenetic signals operating in in high-altitude adjusted indigenous populations may help in therapeutically modulating signaling pathways in hypoxemia patients by copying the most successful epigenotype. In addition, we have summarized the current information about exosomes in hypoxia research and prospects to use them as diagnostic tools to study the epigenome of high-altitude adapted healthy or maladapted individuals.

Application of a novel haplotype-based scan for local adaptation to study high-altitude adaptation in rhesus macaques

When natural populations split and migrate to different environments, they may experience different selection pressures that can lead to local adaptation. To capture the genomic patterns of a local selective sweep, we develop XP-nSL, a genomic scan for local adaptation that compares haplotype patterns between two populations. We show that XP-nSL has power to detect ongoing and recently completed hard and soft sweeps, and we then apply this statistic to search for evidence of adaptation to high altitude in rhesus macaques. We analyze the whole genomes of 23 wild rhesus macaques captured at high altitude (mean altitude > 4000 m above sea level) to 22 wild rhesus macaques captured at low altitude (mean altitude < 500 m above sea level) and find evidence of local adaptation in the high-altitude population at or near 303 known genes and several unannotated regions. We find the strongest signal for adaptation at EGLN1, a classic target for convergent evolution in several species living in low oxygen environments. Furthermore, many of the 303 genes are involved in processes related to hypoxia, regulation of ROS, DNA damage repair, synaptic signaling, and metabolism. These results suggest that, beyond adapting via a beneficial mutation in one single gene, adaptation to high altitude in rhesus macaques is polygenic and spread across numerous important biological systems.

Evolution of the "world's only alpine parrot": Genomic adaptation or phenotypic plasticity, behaviour and ecology?

Climate warming, in particular in island environments, where opportunities for species to disperse are limited, may become a serious threat to cold adapted alpine species. In order to understand how alpine species may respond to a warming world, we need to understand the drivers that have shaped their habitat specialisation and the evolutionary adaptations that allow them to utilize alpine habitats. The endemic, endangered New Zealand kea (Nestor notabilis) is considered the only alpine parrot in the world. As a species commonly found in the alpine zone it may be highly susceptible to climate warming. But is it a true alpine specialist? Is its evolution driven by adaptation to the alpine zone, or is the kea an open habitat generalist that simply uses the alpine zone to, for example, avoid lower lying anthropogenic landscapes? We use whole genome data of the kea and its close, forest adapted sister species, the kākā (Nestor meridionalis) to reconstruct the evolutionary history of both species and identify the functional genomic differences that underlie their habitat specialisations. Our analyses do not identify major functional genomic differences between kea and kākā in pathways associated with high-altitude. Rather, we found evidence that selective pressures on adaptations commonly found in alpine species are present in both Nestor species, suggesting that selection for alpine adaptations has not driven their divergence. Strongly divergent demographic responses to past climate warming between the species nevertheless highlight potential future threats to kea survival in a warming world.

Mitogenomes Reveal Multiple Colonization of Mountains by Rattus in Sundaland

Tropical mountains are cradles of biodiversity and endemism. Sundaland, tropical Southeast Asia, hosts 3 species of Rattus endemic to elevations above 2000 m with an apparent convergence in external morphology: Rattus korinchi and R. hoogerwerfi from Sumatra, and R. baluensis from Borneo. A fourth one, R. tiomanicus, is restricted to lowland elevations across the whole region. The origins of these endemics are little known due to the absence of a robust phylogenetic framework. We use complete mitochondrial genomes from the 3 high altitude Rattus, and several related species to determine their relationships, date divergences, reconstruct their history of colonization, and test for selection on the mitochondrial DNA. We show that mountain colonization happened independently in Borneo (<390 Kya) and Sumatra (~1.38 Mya), likely from lowland lineages. The origin of the Bornean endemic R. baluensis is very recent and its genetic diversity is nested within the diversity of R. tiomanicus. We found weak evidence of positive selection in the high-elevation lineages and attributed the greater nonsynonymous mutations on these branches (specially R. baluensis) to lesser purifying selection having acted on the terminal branches in the phylogeny.

Shifts in the relative fitness contributions of fecundity and survival in variable and changing environments

Organisms respond to shifts in climate means and variability via distinct mechanisms. Accounting for these differential responses and appropriately aggregating them is central to understanding and predicting responses to climate variability and change. Separately considering fitness components can clarify organismal responses: fecundity is primarily an integrated, additive response to chronic environmental conditions over time via mechanisms such as energy use and acquisition, whereas survival can be strongly influenced by short-term, extreme environmental conditions. In many systems, the relative importance of fecundity and survival constraints changes systematically along climate gradients, with fecundity constraints dominating at high latitudes or altitudes (i.e. leading range edges as climate warms), and survival constraints dominating at trailing range edges. Incorporating these systematic differences in models may improve predictions of responses to recent climate change over models that assume similar processes along environmental gradients. We explore how detecting and predicting shifts in fitness constraints can improve our ability to forecast responses to climate gradients and change.

A Re-Assessment of Positive Selection on Mitochondrial Genomes of High-Elevation Phrynocephalus Lizards

Due to their integral roles in oxidative phosphorylation, mitochondrially encoded proteins represent common targets of selection in response to altitudinal hypoxia across high-altitude taxa. While previous studies revealed evidence of positive selection on mitochondrial genomes of high-altitude Phrynocephalus lizards, their conclusions were restricted by out-of-date phylogenies and limited taxonomic sampling. Using topologies derived from both nuclear and mitochondrial DNA phylogenies, we re-assessed the evidence of positive selection on the mitochondrial genomes of high-altitude Phrynocephalus. We sampled representative species from all four main lineages and sequenced the mitochondrial genome of P. maculatus, a putative sister taxon to the high-altitude group. Positive selection was assessed through two widely used branch-site tests: the branch-site model in PAML and BUSTED in HyPhy. No evidence of positive selection on mitochondrial genes was detected on branches leading to two most recent common ancestors of high-altitude species; however, we recovered evidence of positive selection on COX1 on the P. forsythii branch, which represents a reversal from high- to low-elevation environments. A positively selected site therein marked a threonine to valine substitution at position 419. We suggest this bout of selection occurred as the ancestors of P. forsythii re-colonized lower altitude environments north of the Tibetan Plateau. Despite their role in oxidative phosphorylation, we posit that mitochondrial genes are unlikely to have represented historical targets of selection for high-altitude adaptation in Phrynocephalus. Consequently, future studies should address the roles of nuclear genes and differential gene expression.

Molecular convergent and parallel evolution among four high-elevation anuran species from the Tibetan region

BACKGROUND

To date, evidence for the relative prevalence or rarity of molecular convergent and parallel evolution is conflicting, and understanding of how these processes contribute to adaptation is limited. We compared four high-elevation anuran species (Bufo tibetanus, Nanorana parkeri, Rana kukunoris and Scutiger boulengeri) from the Tibetan region, and examined convergent and parallel amino acid substitutions between them and how they may have contributed to high-elevation adaptation.

RESULTS

Genomic data of the four high-elevation species and eight of their low-elevation close relatives were gathered. A total of 1098 orthologs shared by all species were identified. We first conducted pairwise comparisons using Zhang and Kumar's test. Then, the R_conv index was calculated and convergence/divergence correlation plotting was conducted. Furthermore, genes under positive selection and with elevated evolutionary rate were examined. We detected a large number of amino acid sites with convergent or parallel substitutions. Several pairs of high-elevation species, in particular, R. kukunoris vs N. parkeri and B. tibetanus vs S. boulengeri, had excessive amounts of convergent substitutions compared to neutral expectation. Nevertheless, these sites were mostly concentrated in a small number of genes (3-32), and no genome-wide convergence was detected. Furthermore, the majority of these convergent genes were neither under detectable positive selection nor had elevated evolutionary rates, although functional prediction analysis suggested some of the convergent genes could potentially contribute to high-elevation adaptation.

CONCLUSIONS

There is a substantial amount of convergent evolution at the amino-acid level among high-elevation amphibians, although these sites are concentrated in a few genes, not widespread across the genomes. This may attribute to the fact that all the target species are from the same environment. The relative prevalence of convergent substitutions among high-elevation amphibians provides an excellent opportunity for further study of molecular convergent evolution.

Anthropogenic and natural barriers affect genetic connectivity in an Alpine butterfly

Dispersal is a key biological process serving several functions including connectivity among populations. Habitat fragmentation caused by natural or anthropogenic structures may hamper dispersal, thereby disrupting genetic connectivity. Investigating factors affecting dispersal and gene flow is important in the current era of anthropogenic global change, as dispersal comprises a vital part of a species' resilience to environmental change. Using finescale landscape genomics, we investigated gene flow and genetic structure of the Sooty Copper butterfly (Lycaena tityrus) in the Alpine Ötz valley system in Austria. We found surprisingly high levels of gene flow in L. tityrus across the region. Nevertheless, ravines, forests, and roads had effects on genetic structure, while rivers did not. The latter is surprising as roads and rivers have a similar width and run largely in parallel in our study area, pointing towards a higher impact of anthropogenic compared with natural linear structures. Additionally, we detected eleven loci potentially under thermal selection, including ones related to membranes, metabolism, and immune function. This study demonstrates the usefulness of molecular approaches in obtaining estimates of dispersal and population processes in the wild. Our results suggest that, despite high gene flow in the Alpine valley system investigated, L. tityrus nevertheless seems to be vulnerable to anthropogenically-driven habitat fragmentation. With anthropogenic rather than natural linear structures affecting gene flow, this may have important consequences for the persistence of species such as the butterfly studied here in altered landscapes.

The Genome Sequence of Alpine Megacarpaea delavayi Identifies Species-Specific Whole-Genome Duplication

Megacarpaea delavayi (Brassicaceae), a plant found the high mountains of southwest China at high altitudes (3000–4800 m), is used as a vegetable or medicine. Here, we report a draft genome for this species. The assembly genome of M. delavayi is 883 Mb, and 61.59% of the genome is composed of repeat sequences. Annotation of the genome identified a total of 41,114 protein-coding genes. We found that M. delavayi experienced an independent whole-genome duplication (WGD), paralleling those independent WGDs in Iberis, Biscutella, and Anastatica in the early Miocene. Phylogenetic analyses based on the single-copy genes confirmed the position of the genus Megacarpaea within the expanded lineage II of the family and resolved its basal divergence to a subclade consisting of Anastatica, Iberis, and Biscutella. Species-specific and fast-evolving genes in M. delavayi are mainly involved in “DNA repair” and “response to UV-B radiation.” These genetic changes may together help this species survive in high-altitude environments. The reference genome reported here provides a valuable resource for studying adaptation of this and other alpine plants to the high-altitude habitats.

Comparative microRNA Transcriptomes in Domestic Goats Reveal Acclimatization to High Altitude

High-altitude acclimatization is a representative example of vertebrates' acclimatization to harsh and extreme environments. Previous studies reported sufficient evidence for a molecular genetic basis of high-altitude acclimatization, and genomic patterns of genetic variation among populations and species have been widely elucidated in recent years. However, understanding of the miRNA role in high-altitude acclimatization have lagged behind, especially in non-model species. To investigate miRNA expression alterations of goats that were induced by high-altitude stress, we performed comparative miRNA transcriptome analysis on six hypoxia-sensitive tissues (heart, kidney, liver, lung, skeletal muscle, and spleen) in two goat populations from distinct altitudes (600 and 3000 m). We obtained the expression value of 1391 mature miRNAs and identified 138 differentially expressed (DE) miRNAs between high and low altitudes. Combined with tissue specificity analysis, we illustrated alterations of expression levels among altitudes and tissues, and found that there were coexisting tissue-specific and -conserved mechanisms for hypoxia acclimatization. Notably, the interplay between DE miRNA and DE target genes strongly indicated post-transcriptional regulation in the hypoxia inducible factor 1, insulin, and p53 signaling pathways, which might play significant roles in high-altitude acclimatization in domestic goats. It's also worth noting that we experimentally confirmed miR-106a-5p to have a negative regulation effect on angiogenesis by directly targeting FLT-1. These results provide insight into the complicated miRNA expression patterns and regulatory mechanisms of high-altitude acclimatization in domestic goats.

Is there low maternal genetic variation in West Asian populations of leopard?

Persian leopards Panthera pardus saxicolor have been extirpated from over 84% of their historic range and are now limited to rugged landscapes of West Asia and the Caucasus. Understanding and maintaining genetic diversity and population connectivity is important for preventing inbreeding and genetic drift, both of which can threaten population viability. All previous analyses of intraspecific genetic variation of West Asian leopards based on the NADH dehydrogenase subunit 5 gene have reported low mitogenomic diversity. In the current study, we sequenced 959 bp of the mtDNA cytochrome b gene to describe the spatial genetic structure of 22 wild Persian leopards across Iran, which hosts most of the subspecies extant range. The findings based on phylogenetic trees and median-joining network indicated that leopards from Iran formed a distinct subclade, i.e., P. p. saxicolor. The AMOVA analysis showed significant differentiation (88.55%) between the subclades of Persian leopards and other Asian leopards. The lowest levels of haplotype (0.247) and nucleotide (0.00078) diversity were estimated in Persian leopards from Iran. Mitochondrial genome sequencing revealed only two closely related haplotypes. There was no evidence for recent sudden demographic expansion scenario in Persian leopards. The low diversity in cytochrome b gene could potentially be brought about by selective pressure on mitochondria to adapt to oxidative stress and higher metabolic rates in cold environments.

Analysis of Multiplicity of Hypoxia-Inducible Factors in the Evolution of Triplophysa Fish (Osteichthyes: Nemacheilinae) Reveals Hypoxic Environments Adaptation to Tibetan Plateau

HIF (Hypoxia-inducible factor) gene family members function as master regulators of cellular and systemic oxygen homeostasis during changes in oxygen availability. Qinghai-Tibet Plateau is a natural laboratory for for long-term hypoxia and cold adaptation. In this context, T. scleroptera that is restricted to >3500 m high-altitude freshwater rivers was selected as the model to compare with a representative species from the plain, P. dabryanus. We cloned different HIF-α and carried out a phylogenetic analysis from invertebrates to vertebrates for identifying HIF-α genes and analyzing their evolutionary history. Intriguingly, the HIF-α has undergone gene duplications might be due to whole-genome duplication (WGD) events during evolution. PAML analysis indicated that HIF-1αA was subjected to positive selection acted on specific sites in Triplophysa lineages. To investigate the relationship between hypoxia adaptation and the regulation of HIF-α stability by pVHL in plateau and plain fish, a series of experiments were carried out. Comparison the luciferase transcriptional activity and protein levels of HIF-αs and the differing interactions of HIF-αs with pVHL, show clear differences between plateau and plain fish. T. scleroptera pVHL could enhance HIF-α transcriptional activity under hypoxia, and functional validation through pVHL protein mutagenesis showed that these mutations increased the stability of HIF-α and its hetero dimerization affinity to ARNT. Our research shows that missense mutations of pVHL induced evolutionary molecular adaptation in Triplophysa fishes living in high altitude hypoxic environments.

Potential sex-dependent effects of weather on apparent survival of a high-elevation specialist

Mountain ecosystems are inhabited by highly specialised and endemic species which are particularly susceptible to climatic changes. However, the mechanisms by which climate change affects species population dynamics are still largely unknown, particularly for mountain birds. We investigated how weather variables correlate with survival or movement of the white-winged snowfinch Montifringilla nivalis, a specialist of high-elevation habitat. We analysed a 15-year (2003-2017) mark-recapture data set of 671 individuals from the Apennines (Italy), using mark-recapture models. Mark-recapture data allow estimating, forgiven time intervals, the probability that individuals stay in the study area and survive, the so called apparent survival. We estimated annual apparent survival to be around 0.44-0.54 for males and around 0.51-0.64 for females. Variance among years was high (range: 0.2-0.8), particularly for females. Apparent survival was lower in winter compared to summer. Female annual apparent survival was negatively correlated with warm and dry summers, whereas in males these weather variables only weakly correlated with apparent survival. Remarkably, the average apparent survival measured in this study was lower than expected. We suggest that the low apparent survival may be due to recent changes in the environment caused by global warming. Possible, non-exclusive mechanisms that potentially also could explain sexual differential apparent survival act via differential breeding dispersal, hyperthermia, weather-dependent food availability, and weather-dependent trade-off between reproduction and self-maintenance. These results improve our current understanding of the mechanisms driving population dynamics in high-elevation specialist birds, which are particularly at risk due to climate change.

Hypoxia and cancer related pathology

Hypoxic environments occur normally at high altitude, or in underground burrows and in deep sea habitats. They also occur pathologically in human ischemia and in hypoxic solid tumors. Hypoxia in various cancer types and its related molecular mechanisms are associated with a poor clinical outcome. This review will discuss how hypoxia can influence two aspects of tumorigenesis, namely the direct, cell-intrinsic oncogenic effects, as well as the indirect effects on tumor progression mediated by an altered tumor microenvironment. We will also discuss recent progress in identifying the functional roles of hypoxia-related factors (HIFs), along with their regulators and downstream target genes, in cancer stem cells and therapy. Importantly, we propose, using convergent evolution schemes to identify novel biomarkers for both hypoxia adaptation and hypoxic solid tumors as an important strategy in the future.

A linked-read approach to museomics: Higher quality de novo genome assemblies from degraded tissues

High-throughput sequencing technologies are a proposed solution for accessing the molecular data in historical specimens. However, degraded DNA combined with the computational demands of short-read assemblies has posed significant laboratory and bioinformatics challenges for de novo genome assembly. Linked-read or "synthetic long-read" sequencing technologies, such as 10× Genomics, may provide a cost-effective alternative solution to assemble higher quality de novo genomes from degraded tissue samples. Here, we compare assembly quality (e.g., genome contiguity and completeness, presence of orthogroups) between four new deer mouse (Peromyscus spp.) genomes assembled using linked-read technology and four published genomes assembled from a single shotgun library. At a similar price-point, these approaches produce vastly different assemblies, with linked-read assemblies having overall higher contiguity and completeness, measured by larger N50 values and greater number of genes assembled, respectively. As a proof-of-concept, we used annotated genes from the four Peromyscus linked-read assemblies and eight additional rodent taxa to generate a phylogeny, which reconstructed the expected relationships among species with 100% support. Although not without caveats, our results suggest that linked-read sequencing approaches are a viable option to build de novo genomes from degraded tissues, which may prove particularly valuable for taxa that are extinct, rare or difficult to collect.

Effect of hemoglobin polymorphism on performance traits in indigenous chicken genotypes in Nigeria

The effect of hemoglobin polymorphism on performance traits in Nigerian indigenous chicken types was investigated in this research. The chickens were obtained through the pure mating of the sire and dam of each population of frizzle feathered, normal feathered, and naked neck chickens to produce F₁ offsprings. One hundred fifty-five chicks (37 frizzle, 79 normal, and 39 naked neck) were measured for body weight (g), breast girth (cm), and tibia length (cm). At 20 weeks, 5 mL of blood was collected from the wing vein of each chicken into heparinized tubes and labeled according to its tag number for electrophoresis. The electrophoresis procedures outlined by RIKEN BRC (2006) was used. 0.6 μl of undiluted blood constituent was taken. The prepared buffer (Tris 10.91 g, EDTA 0.60 g, boric acid 3.10 g) at pH 8.5 was used for the gel preparation. The cellulose acetate membrane was used as a supporting medium. The electrophoresis was carried out at a voltage of 150 V for about 50 min at a temperature of 4 °C. The migration of the genotype was from cathode (-) to anode (+). Each bird was scored as either fast (AA), midway (AB), or slow (BB) according to the mobility on the cellulose acetate paper for hemoglobin. Data were analyzed using statistical analysis system (SAS 2002), and significant means were separated using Tukey's Honestly Significant Difference. Hardy Weinberg's equation was used to calculate genotypic and allelic frequencies and tested using chi-square (χ²). Hemoglobin had three polymorphic forms viz AA, AB, and BB. The genotypic frequencies of hemoglobin (Hb^AA, Hb^AB, and Hb^BB) were 49, 56, and 50, respectively, while the allelic frequencies were 0.50 for both Hb^A and Hb^B. The effect of the polymorphic forms on body weight (g), breast girth (cm), and tibia length (cm) showed that the AA had significantly higher (P < 0.05) body weight (g) than AB and BB (1296.43 g, 1029.59 g, and 884.46 g, respectively). The AA was also higher (P < 0.05) than the AB and BB for breast girth (cm) and tibia length (cm). Heterozygotes adapted and survived better than the homozygotes. The effect of the polymorphic forms of hemoglobin on body weight (g), breast girth (cm), and tibia length (cm) showed that it could be used for body weight selection.

Phylogeny and highland adaptation of Chinese species in Allium section Daghestanica (Amaryllidaceae) revealed by transcriptome sequencing

Allium L. is one of the largest monocotyledonous genera with extensive distribution in the Northern Hemisphere. The fundamental phylogenies of Allium have been investigated using many morphological and molecular characters. However, the morphological characters may not agree with the molecular results in some Allium groups or sections (such as the Chinese Allium section Daghestanica), which may result in ambiguous species relationships and hinder further evolutionary and adaptive researches. Here, transcriptome sequences of the six Chinese endemics from Allium section Daghestanica were collected, with their single-copy genes (SCGs) were extracted. The interspecies relationships were analyzed using concatenation and coalescent methods. The branch-site model (BSM) was conducted to detect the positively selected genes (PSGs) in five highland species of this section. Based on 1644, 1281 and 1580 SCGs in flowers, leaves, and flowers-leaves combination respectively, a robust consistent and well-resolved phylogeny was generated from the concatenation method. Strong conflicts among individual gene trees were detected in the coalescent method, and morphological characters were incongruent with molecular relationships to some degree. Many PSGs were involved in responses of various stresses and stimuli (e.g. hypoxia, low temperature, aridity), DNA repair, metabolism, nutrient or energy intake, photosynthesis, and signal transduction. Our study revealed a clear interspecies relationship of Chinese endemics in Allium section Daghestanica and suggested that the discordance between morphological characters and molecular relationships might result from that the former are more susceptible to convergence compared with the latter. PSGs detected in our study may provide some insights into highland adaptation in Allium species.

Convergent genomic signatures of high-altitude adaptation among domestic mammals

Abundant and diverse domestic mammals living on the Tibetan Plateau provide useful materials for investigating adaptive evolution and genetic convergence. Here, we used 327 genomes from horses, sheep, goats, cattle, pigs and dogs living at both high and low altitudes, including 73 genomes generated for this study, to disentangle the genetic mechanisms underlying local adaptation of domestic mammals. Although molecular convergence is comparatively rare at the DNA sequence level, we found convergent signature of positive selection at the gene level, particularly the EPAS1 gene in these Tibetan domestic mammals. We also reported a potential function in response to hypoxia for the gene C10orf67, which underwent positive selection in three of the domestic mammals. Our data provide an insight into adaptive evolution of high-altitude domestic mammals, and should facilitate the search for additional novel genes involved in the hypoxia response pathway.

Parallel Molecular Evolution in Pathways, Genes, and Sites in High-Elevation Hummingbirds Revealed by Comparative Transcriptomics

High-elevation organisms experience shared environmental challenges that include low oxygen availability, cold temperatures, and intense ultraviolet radiation. Consequently, repeated evolution of the same genetic mechanisms may occur across high-elevation taxa. To test this prediction, we investigated the extent to which the same biochemical pathways, genes, or sites were subject to parallel molecular evolution for 12 Andean hummingbird species (family: Trochilidae) representing several independent transitions to high elevation across the phylogeny. Across high-elevation species, we discovered parallel evolution for several pathways and genes with evidence of positive selection. In particular, positively selected genes were frequently part of cellular respiration, metabolism, or cell death pathways. To further examine the role of elevation in our analyses, we compared results for low- and high-elevation species and tested different thresholds for defining elevation categories. In analyses with different elevation thresholds, positively selected genes reflected similar functions and pathways, even though there were almost no specific genes in common. For example, EPAS1 (HIF2α), which has been implicated in high-elevation adaptation in other vertebrates, shows a signature of positive selection when high-elevation is defined broadly (>1,500 m), but not when defined narrowly (>2,500 m). Although a few biochemical pathways and genes change predictably as part of hummingbird adaptation to high-elevation conditions, independent lineages have rarely adapted via the same substitutions.