Divergent adaptation to Qinghai-Tibetan Plateau implicated from transciptome study of Gymnocypris dobula and Schizothorax nukiangensis

The Schizothoracine fishes are widely distributed in the Qinghai-Tibetan Plateau (QTP) area and its peripheral regions, which provide a prime example of adaptation in highland aquatic environments. Recent progresses have revealed various genetic adaptations of these fishes by comparing to distantly related lowerland species, however, comparative studies on closely-related species of different altitudes are still lacking. In this study, we sequenced and annotated a primitive Schizothoracine fish Schizothorax nukiangensis Tsao and a highly specialized one Gymnocypris dobula. We performed evolutionary analyses to investigate the candidate genes and signaling pathways involved QTP highland adaptation in both Schizothoracine fishes. Analysis of the 11,007 one-copy orthologs to the primitive cyprinid species, Danio rerio, revealed that both G. dobula and S. nukiangensis showed elevated evolutionary rates. A large number of genes related to hypoxia, including genes involved metabolic processes and cardiovascular system development, exhibited signatures of positive selection in both Schizothoracine fishes, but very few positively selected genes were found overlapping among these Schizothoracines. Our results indicated divergent genetic adaptation to highland environment for aquatic species living in QTP.     

Whole genome sequencing of Ethiopian highlanders reveals conserved hypoxia tolerance genes

Background

Although it has long been proposed that genetic factors contribute to adaptation to high altitude, such factors remain largely unverified. Recent advances in high-throughput sequencing have made it feasible to analyze genome-wide patterns of genetic variation in human populations. Since traditionally such studies surveyed only a small fraction of the genome, interpretation of the results was limited.

Results

We report here the results of the first whole genome resequencing-based analysis identifying genes that likely modulate high altitude adaptation in native Ethiopians residing at 3,500 m above sea level on Bale Plateau or Chennek field in Ethiopia. Using cross-population tests of selection, we identify regions with a significant loss of diversity, indicative of a selective sweep. We focus on a 208 kbp gene-rich region on chromosome 19, which is significant in both of the Ethiopian subpopulations sampled. This region contains eight protein-coding genes and spans 135 SNPs. To elucidate its potential role in hypoxia tolerance, we experimentally tested whether individual genes from the region affect hypoxia tolerance in Drosophila. Three genes significantly impact survival rates in low oxygen: cic, an ortholog of human CIC, Hsl, an ortholog of human LIPE, and Paf-AHα, an ortholog of human PAFAH1B3.

Conclusions

Our study reveals evolutionarily conserved genes that modulate hypoxia tolerance. In addition, we show that many of our results would likely be unattainable using data from exome sequencing or microarray studies. This highlights the importance of whole genome sequencing for investigating adaptation by natural selection.     

Identifying Signatures of Natural Selection in Tibetan and Andean Populations Using Dense Genome Scan Data

High-altitude hypoxia (reduced inspired oxygen tension due to decreased barometric pressure) exerts severe physiological stress on the human body. Two high-altitude regions where humans have lived for millennia are the Andean Altiplano and the Tibetan Plateau. Populations living in these regions exhibit unique circulatory, respiratory, and hematological adaptations to life at high altitude. Although these responses have been well characterized physiologically, their underlying genetic basis remains unknown. We performed a genome scan to identify genes showing evidence of adaptation to hypoxia. We looked across each chromosome to identify genomic regions with previously unknown function with respect to altitude phenotypes. In addition, groups of genes functioning in oxygen metabolism and sensing were examined to test the hypothesis that particular pathways have been involved in genetic adaptation to altitude. Applying four population genetic statistics commonly used for detecting signatures of natural selection, we identified selection-nominated candidate genes and gene regions in these two populations (Andeans and Tibetans) separately. The Tibetan and Andean patterns of genetic adaptation are largely distinct from one another, with both populations showing evidence of positive natural selection in different genes or gene regions. Interestingly, one gene previously known to be important in cellular oxygen sensing, EGLN1 (also known as PHD2), shows evidence of positive selection in both Tibetans and Andeans. However, the pattern of variation for this gene differs between the two populations. Our results indicate that several key HIF-regulatory and targeted genes are responsible for adaptation to high altitude in Andeans and Tibetans, and several different chromosomal regions are implicated in the putative response to selection. These data suggest a genetic role in high-altitude adaption and provide a basis for future genotype/phenotype association studies necessary to confirm the role of selection-nominated candidate genes and gene regions in adaptation to altitude.     

A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.     

Genomic Scan Reveals Loci under Altitude Adaptation in Tibetan and Dahe Pigs

High altitude environments are of particular interest in the studies of local adaptation as well as their implications in physiology and clinical medicine in human. Some Chinese pig breeds, such as Tibetan pig (TBP) that is well adapted to the high altitude and Dahe pig (DHP) that dwells at the moderate altitude, provide ideal materials to study local adaptation to altitudes. Yet, it is still short of in-depth analysis and understanding of the genetic adaptation to high altitude in the two pig populations. In this study we conducted a genomic scan for selective sweeps using F_ST to identify genes showing evidence of local adaptations in TBP and DHP, with Wuzhishan pig (WZSP) as the low-altitude reference. Totally, we identified 12 specific selective genes (CCBE1, F2RL1, AGGF1, ZFPM2, IL2, FGF5, PLA2G4A, ADAMTS9, NRBF2, JMJD1C, VEGFC and ADAM19) for TBP and six (OGG1, FOXM, FLT3, RTEL1, CRELD1 and RHOG) for DHP. In addition, six selective genes (VPS13A, GNA14, GDAP1, PARP8, FGF10 and ADAMTS16) were shared by the two pig breeds. Among these selective genes, three (VEGFC, FGF10 and ADAMTS9) were previously reported to be linked to the local adaptation to high altitudes in pigs, while many others were newly identified by this study. Further bioinformatics analysis demonstrated that majority of these selective signatures have some biological functions relevant to the altitude adaptation, for examples, response to hypoxia, development of blood vessels, DNA repair and several hematological involvements. These results suggest that the local adaptation to high altitude environments is sophisticated, involving numerous genes and multiple biological processes, and the shared selective signatures by the two pig breeds may provide an effective avenue to identify the common adaptive mechanisms to different altitudes.     

Induction of SENP1 in Endothelial Cells Contributes to Hypoxia-driven VEGF Expression and Angiogenesis

SENP1 (SUMO-specific protease 1) has been shown to be essential for the stability and activity of hypoxia-inducible factor 1 (HIF-1α) under hypoxia conditions. However, it is unknown how SENP1 activation and hypoxia signaling are coordinated in the cellular response to hypoxia. Here, we report the essential role of SENP1 in endothelial cells as a positive regulator of hypoxia-driven VEGF production and angiogenesis. SENP1 expression is increased in endothelial cells following exposure to hypoxia. Silencing of HIF-1α blocks SENP1 expression in cell response to hypoxia. Mutation of the hypoxia response element (HRE) on the Senp1 promoter abolishes its transactivation in response to hypoxia. Moreover, silencing of SENP1 expression decreases VEGF production and abrogates the angiogenic functions of endothelial cell. We also find that the elongated endothelial cells in embryonic brain section and vascular endothelial cells in embryonic renal glomeruli in Senp1^−/− mice are markedly reduced than those in wild-type. Thus, these results show that hypoxia implies a positive feedback loop mediated by SENP1. This feedback loop is important in VEGF production, which is essential for angiogenesis in endothelial cells.     

BALB/c-congenic ANP32B-deficient mice reveal a modifying locus that determines viability

We previously found that deletion of the multifunctional factor ANP32B (a.k.a. SSP29, APRIL, PAL31, PHAPI2) resulted in a severe but strain-specific defect resulting in perinatal lethality. The difficulty in generating an adult cohort of ANP32B-deficient animals limited our ability to examine adult phenotypes, particularly cancer-related phenotypes. We bred the Anp32b-null allele into the BALB/c and FVB/N genetic background. The BALB/c, but not the FVB/N, background provided sufficient frequency of adult Anp32b-null (Anp32b^−/−) animals. From these, we found no apparent oncogenic role for this protein in mammary tumorigenesis contrary to what was predicted based on human data. We also found runtism, pathologies in various organ systems, and an unusual clinical chemistry signature in the adult Anp32b^−/− mice. Intriguingly, genome-wide single-nucleotide polymorphism analysis suggested that our colony retained an unlinked C57BL/6J locus at high frequency. Breeding this locus to homozygosity demonstrated that it had a strong effect on Anp32b^−/− viability indicating that this locus contains a modifier gene of Anp32b with respect to development. This suggests a functionally important genetic interaction with one of a limited number of candidate genes, foremost among them being the variant histone gene H2afv. Using congenic breeding strategies, we have generated a viable ANP32B-deficient animal in a mostly pure background. We have used this animal to reliably exclude mouse ANP32B as an important oncogene in mammary tumorigenesis. Our further phenotyping strengthens the evidence that ANP32B is a widespread regulator of gene expression. These studies may also impact the choice of subsequent groups with respect to congenic breeding versus de novo zygote targeting strategies for background analyses in mouse genetics.     

Mitochondrial OXPHOS genes provides insights into genetics basis of hypoxia adaptation in anchialine cave shrimps

Cave shrimps from the genera Typhlatya, Stygiocaris and Typhlopatsa (TST complex) comprises twenty cave-adapted taxa, which mainly occur in the anchialine environment. Anchialine habitats may undergo drastic environmental fluctuations, including spatial and temporal changes in salinity, temperature, and dissolved oxygen content. Previous studies of crustaceans from anchialine caves suggest that they have possessed morphological, behavioral, and physiological adaptations to cope with the extreme conditions, similar to other cave-dwelling crustaceans. However, the genetic basis has not been thoroughly explored in crustaceans from anchialine habitats, which can experience hypoxic regimes. To test whether the TST shrimp-complex hypoxia adaptations matched adaptive evolution of mitochondrial OXPHOS genes. The 13 OXPHOS genes from mitochondrial genomes of 98 shrimps and 1 outgroup were examined. For each of these genes was investigated and compared to orthologous sequences using both gene (i.e. branch-site and Datamonkey) and protein (i.e. TreeSAAP) level approaches. Positive selection was detected in 11 of the 13 candidate genes, and the radical amino acid changes sites scattered throughout the entire TST complex phylogeny. Additionally, a series of parallel/convergent amino acid substitutions were identified in mitochondrial OXPHOS genes of TST complex shrimps, which reflect functional convergence or similar genetic mechanisms of cave adaptation. The extensive occurrence of positive selection is suggestive of their essential role in adaptation to hypoxic anchialine environment, and further implying that TST complex shrimps might have acquired a finely capacity for energy metabolism. These results provided some new insights into the genetic basis of anchialine hypoxia adaptation.     

Comparative analysis of fertility restoration genes for WA,Y, and DA cytoplasmic male sterility in rice

Rice chromosome single segment substitution line (SSSL) W23-19-06-06-11 with the genotype Rf3Rf3/Rf4Rf4, a strong restorer line for wild-abortive (WA) cytoplasmic male sterility (CMS), was recently identified from the SSSL library. To investigate the genetic mode of Rf genes and the genetic relationship among WA, yegong (Y), and dwarf-wildabortive (DA) CMS systems, the plants derived from three BC3F2 populations involving W23-19-06-06-11 and the three CMS lines, that carried the Rf3Rf3/Rf4Rf4, Rf3Rf3/rf4rf4, and rf3rf3/Rf4Rf4 genotypes and WA-, Y-, and DA-CMS cytoplasm, were selected and their pollen and spikelet fertility were evaluated. The results show that the genetic effect displayed a trend of Y-CMS > WA-CMS > DA-CMS in the genetic background of W23-19-06-06-11, the effect of Rf4 appeared to be slightly larger than that of Rf3, and their effects were additive for the three CMS systems. Two pairs of dominant genes governed the fertility restoration in pollen and spikelet in the W23-19-06-06-11 which indicates that the genetic mode of the Rf genes was a qualitative character for the three CMS systems.     

Molecular basis of cytoplasmic male sterility and fertility restoration in rice

Cytoplasmic male sterility (CMS) is a maternally inherited trait that causes dysfunctions in pollen and anther development. CMS is caused by the interaction between nuclear and mitochondrial genomes. A product of a CMS-causing gene encoded by the mitochondrial genome affects mitochondrial function and the regulation of nuclear genes, leading to male sterility. In contrast, the RESTORER OF FERTILITY gene (Rf gene) in the nuclear genome suppresses the expression of the CMS-causing gene and restores male fertility. An alloplasmic CMS line is often bred as a result of nuclear substitution, which causes the removal of functional Rf genes and allows the expression of a CMS-causing gene in mitochondria. The CMS/Rf system is an excellent model for understanding the genetic interactions and cooperative functions of mitochondrial and nuclear genomes in plants, and is also an agronomically important trait for hybrid seed production. In this review article, pollen and anther phenotypes of CMS, CMS-associated mitochondrial genes, Rf genes, and the mechanism that causes pollen abortion and its agronomical application for rice are described.     

Hypoxia adaptation and hemoglobin mutation in Tibetan chick embryo

Tibetan chick lives at high altitudes between 2600 and 4200 m with a high hatchability and low land breeds survive rarely with a hatchability of 3.0% under hypoxia of simulated 4200 m. Under hypoxia of whole 21 d, the hatchability of Tibetan chick and Recessive White Feather broiler differed with a greatest disparity from day 4 to 11 and also significantly in other stages except from day 1 to 3. Hypoxia in each stage did not reduce significantly survival rate of this stage except hatchability. These two results indicated that the hypoxia in the early stage had an adverse effect on the later stage. All exons encoding chick hemoglobins were sequenced to analyze gene polymorphism. The functional mutation Met-32(B13)-Leu, related with hypoxia, was found in α^D globin chain and the mutation frequency increased with increased altitude. In addition, under hypoxic conditions, the population with higher mutation frequency had a higher hatchability. The automated homology model building was carried out using crystal structure coordinates of chick HbD. The results indicated that the substitution Met-32(B13)-Leu provides a more hydrophobic environment which leads to higher stability of heme and oxygen affinity of hemoglobin. The occurrence of the mutation Met-32(B13)-Leu is related to the origin of Tibetan chick.     

Relevant genetic differentiation among Brazilian populations of Anastrepha fraterculus (Diptera,Tephritidae)

We used a population genetic approach to detect the presence of genetic diversity among six populations of Anastrepha fraterculus across Brazil. To this aim, we used Simple Sequence Repeat (SSR) markers, which may capture the presence of differentiative processes across the genome in distinct populations. Spatial analyses of molecular variance were used to identify groups of populations that are both genetically and geographically homogeneous while also being maximally differentiated from each other. The spatial analysis of genetic diversity indicates that the levels of diversity among the six populations vary significantly on an eco-geographical basis. Particularly, altitude seems to represent a differentiating adaptation, as the main genetic differentiation is detected between the two populations present at higher altitudes and the other four populations at sea level. The data, together with the outcomes from different cluster analyses, identify a genetic diversity pattern that overlaps with the distribution of the known morphotypes in the Brazilian area.     

Features of the population and genetic structure of Norway spruce (Picea abies (L.) Karst.) from the Ukrainian Carpathian Mountains and Polesie

Genetic differences between a lowland, island, and sparse population (2 ha) of Norway spruce (Picea abies (L.) Karst.) from Polesie and two populations at different altitudes in the Ukrainian Carpathian Mountains (in the range 1000–1500 m above sea level) are identified on the basis of an analysis of the variability of 20 allozyme loci. Despite the different sizes of the samples, different altitude location, and geographic remoteness, these populations are characterized by very low genetic differentiation (D _N = 0.002–0.005). Recombination variability analysis demonstrated that, in terms of the structure of double-locus genotypes, the Polesie population possesses lesser affinity than do the individual Carpathian populations with respect to each other.     

The Genomic Architecture of Adaptation to Larval Malnutrition Points to a Trade-off with Adult Starvation Resistance in Drosophila

Periods of nutrient shortage impose strong selection on animal populations. Experimental studies of genetic adaptation to nutrient shortage largely focus on resistance to acute starvation at adult stage; it is not clear how conclusions drawn from these studies extrapolate to other forms of nutritional stress. We studied the genomic signature of adaptation to chronic juvenile malnutrition in six populations of Drosophila melanogaster evolved for 150 generations on an extremely nutrient-poor larval diet. Comparison with control populations evolved on standard food revealed repeatable genomic differentiation between the two set of population, involving >3,000 candidate SNPs forming >100 independently evolving clusters. The candidate genomic regions were enriched in genes implicated in hormone, carbohydrate, and lipid metabolism, including some with known effects on fitness-related life-history traits. Rather than being close to fixation, a substantial fraction of candidate SNPs segregated at intermediate allele frequencies in all malnutrition-adapted populations. This, together with patterns of among-population variation in allele frequencies and estimates of Tajima’s D, suggests that the poor diet results in balancing selection on some genomic regions. Our candidate genes for tolerance to larval malnutrition showed a high overlap with genes previously implicated in acute starvation resistance. However, adaptation to larval malnutrition in our study was associated with reduced tolerance to acute adult starvation. Thus, rather than reflecting synergy, the shared genomic architecture appears to mediate an evolutionary trade-off between tolerances to these two forms of nutritional stress.