Natural selection and adaptive evolution of leptin in the ochotona family driven by the cold environmental stress

Evidence of Intraspecific Adaptive Variation in the American Pika (Ochotona princeps) on a Continental Scale Using a Target Enrichment and Mitochondrial Genome Skimming Approach

Montane landscapes present an array of abiotic challenges that drive adaptive evolution amongst organisms. These adaptations can promote habitat specialisation, which may heighten the risk of extirpation from environmental change. For example, higher metabolic rates in an endothermic species may contribute to heightened cold tolerance, whilst simultaneously limiting heat tolerance. Here, using the climate-sensitive American pika (Ochotona princeps), we test for evidence of intraspecific adaptive variation amongst environmental gradients across the Intermountain West of North America. We leveraged results from previous studies on pika adaptation to generate a custom nuclear target enrichment design to sequence several hundred candidate genes related to cold, hypoxia and dietary detoxification. We also applied a 'genome skimming' approach to sequence mitochondrial DNA. Using genotype-environment association tests, we identified rare genomic variants associated with elevation and temperature variation amongst populations. Amongst mitochondrial genes, we identified intraspecific variation in selective signals and significant changes to the amino acid property equilibrium constant, which may relate to electron transport chain efficiency. These results illustrate a complex dynamic of adaptive variation amongst O. princeps where lineages and populations have adapted to unique regional conditions. Some of the clearest signals of selection were in a genetic lineage that includes pikas of the Great Basin region, which is also where recent localised extirpations have taken place and highlights the risk of losing adaptive alleles during environmental change.

Adaptive evolution and co-evolution of chloroplast genomes in Pteridaceae species occupying different habitats: overlapping residues are always highly mutated

BACKGROUND

The evolution of protein residues depends on the mutation rates of their encoding nucleotides, but it may also be affected by co-evolution with other residues. Chloroplasts function as environmental sensors, transforming fluctuating environmental signals into different physiological responses. We reasoned that habitat diversity may affect their rate and mode of evolution, which might be evidenced in the chloroplast genome. The Pteridaceae family of ferns occupy an unusually broad range of ecological niches, which provides an ideal system for analysis.

RESULTS

We conducted adaptive evolution and intra-molecular co-evolution analyses of Pteridaceae chloroplast DNAs (cpDNAs). The results indicate that the residues undergoing adaptive evolution and co-evolution were mostly independent, with only a few residues being simultaneously involved in both processes, and these overlapping residues tend to exhibit high mutations. Additionally, our data showed that Pteridaceae chloroplast genes are under purifying selection. Regardless of whether we grouped species by lineage (which corresponded with ecological niches), we determined that positively selected residues mainly target photosynthetic genes.

CONCLUSIONS

Our work provides evidence for the adaptive evolution of Pteridaceae cpDNAs, especially photosynthetic genes, to different habitats and sheds light on the adaptive evolution and co-evolution of proteins.

In silico analyses of leptin and leptin receptor of spotted snakehead Channa punctata

The present study, in addition to molecular characterization of leptin (lepa) and its receptor (lepr) of spotted snakehead Channa punctata, is focussed on physicochemical, structural, evolutionary and selection pressure analyses which are poorly elucidated in teleosts in spite of that existence of these genes is well reported in several fish species. The putative full-length Lep and Lepr of C. punctata showed conserved structural and functional domains, especially the residues responsible for structural integrity and signal transduction. Conversely, residues predicted essential for Lep-Lepr interaction displayed divergence between teleosts and tetrapods. Impact of substitutions/deletions predicted using protein variation effect analyser tool highlighted species specificity in ligand-receptor interaction. Physicochemical properties of ligand and receptor predicted for the first time in vertebrates revealed high aliphatic and instability indices for both Lepa and Lepr, indicating thermostability of proteins but their instability under ex vivo conditions. Positive grand average of hydropathy score of Lepa suggests its hydrophobic nature conjecturing existence of leptin binding proteins in C. punctata. In addition to disulphide bonding, a novel posttranslational modification (S-126 phosphorylation) was predicted in Lepa of C. punctata. In Lepr, disulphide bond formation and N-linked glycosylation near WSXWS motif in ECD, and phosphorylation at tyrosine residues in ICD were predicted. Leptin and its receptor sequence of C. punctata cladded with its homolog from C. striata and C. argus of order Anabantiformes. Leptin system of Anabantiformes was phylogenetically closer to that of Pleuronectiformes, Scombriformes and Perciformes. Selection pressure analysis showed higher incidence of negative selection in teleostean leptin genes indicating limited adaptation in their structure and function. However, evidence of pervasive and episodic diversifying selection laid a foundation of co-evolution of Lepa and Lepr in teleosts.

Leptin Gene Protects Against Cold Stress in Antarctic Toothfish

Leptin is a cytokine-like peptide, predominantly biosynthesized in adipose tissue, which plays an important role in regulating food intake, energy balance and reproduction in mammals. However, how it may have been modified to enable life in the chronic cold is unclear. Here, we identified a leptin-a gene (lepa) in the cold-adapted and neutrally buoyant Antarctic toothfish Dissostichus mawsoni that encodes a polypeptide carrying four α-helices and two cysteine residues forming in-chain disulfide bonds, structures shared by most vertebrate leptins. Quantitative RT-PCR confirmed that mRNA levels of the leptin-a gene of D. mawsoni (DM-lepa) were highest in muscle, followed by kidney and liver; detection levels were low in the gill, brain, intestine, and ovary tissues. Compared with leptin-a genes of fishes living in warmer waters, DM-lepa underwent rapid evolution and was subjected to positive selection. Over-expression of DM-lepa in the zebrafish cell line ZFL resulted in signal accumulation in the cytoplasm and significantly increased cell proliferation both at the normal culture temperature and under cold treatment. DM-lepa over-expression also reduced apoptosis under low-temperature stress and activated the STAT3 signaling pathway, in turn upregulating the anti-apoptotic proteins bcl2l1, bcl2a, myca and mdm2 while downregulating the pro-apoptotic baxa, p53 and caspase-3. These results demonstrate that DM-lepa, through STAT3 signaling, plays a protective role in cold stress by preventing apoptotic damage. Our study reveals a new role of lepa in polar fish.

Building a Chinese pan-genome of 486 individuals

Pan-genome sequence analysis of human population ancestry is critical for expanding and better defining human genome sequence diversity. However, the amount of genetic variation still missing from current human reference sequences is still unknown. Here, we used 486 deep-sequenced Han Chinese genomes to identify 276 Mbp of DNA sequences that, to our knowledge, are absent in the current human reference. We classified these sequences into individual-specific and common sequences, and propose that the common sequence size is uncapped with a growing population. The 46.646 Mbp common sequences obtained from the 486 individuals improved the accuracy of variant calling and mapping rate when added to the reference genome. We also analyzed the genomic positions of these common sequences and found that they came from genomic regions characterized by high mutation rate and low pathogenicity. Our study authenticates the Chinese pan-genome as representative of DNA sequences specific to the Han Chinese population missing from the GRCh38 reference genome and establishes the newly defined common sequences as candidates to supplement the current human reference.

Higher Expression Levels of Aquaporin Family of Proteins in the Kidneys of Arid-Desert Living Lepus yarkandensis

Lepus yarkandensis specifically lives in arid climate with rare precipitation of Tarim Basin in western China. Aquaporins (AQPs) are a family of channel proteins that facilitate water transportation across cell membranes. Kidney AQPs play vital roles in renal tubule water permeability and maintenance of body water homeostasis. This study aimed to investigate whether kidney AQPs exhibit higher expression in arid-desert living animals. Immunohistochemistry results revealed localization of AQP1 to the capillary endothelial cells in glomerulus and epithelial cells in proximal tubule and descending thin limbs, AQP2 to the apical plasma membrane of principal cells in the cortical collecting duct (CCD), outer medullary collecting duct (OMCD), and IMCD cells in the initial inner medullary collecting duct (IMCD1) and middle IMCD (IMCD2), and AQP3 and AQP4 to the basolateral plasma membrane of principal cells and IMCD cells in CCD, OMCD, IMCD1, and IMCD2 in L. yarkandensis kidneys. Quantitative real-time PCR analysis showed higher mRNA levels of AQP1, AQP2, AQP3, and AQP4 in L. yarkandensis kidneys compared with Oryctolagus cuniculus. Similar results were obtained by western blotting. Our results suggested that higher expression levels of AQP1, AQP2, AQP3, and AQP4 in L. yarkandensis kidneys favored for drawing more water from the tubular fluid.

Hibernation in bats (Mammalia: Chiroptera) did not evolve through positive selection of leptin

Temperature regulation is an indispensable physiological activity critical for animal survival. However, relatively little is known about the origin of thermoregulatory regimes in a phylogenetic context, or the genetic mechanisms driving the evolution of these regimes. Using bats as a study system, we examined the evolution of three thermoregulatory regimes (hibernation, daily heterothermy, and homeothermy) in relation to the evolution of leptin, a protein implicated in regulation of torpor bouts in mammals, including bats. A threshold model was used to test for a correlation between lineages with positively selected lep, the gene encoding leptin, and the thermoregulatory regimes of those lineages. Although evidence for episodic positive selection of lep was found, positive selection was not correlated with lineages of heterothermic bats, a finding that contradicts results from previous studies. Evidence from our ancestral state reconstructions suggests that the most recent common ancestor of bats used daily heterothermy and that the presence of hibernation is highly unlikely at this node. Hibernation likely evolved independently at least four times in bats-once in the common ancestor of Vespertilionidae and Molossidae, once in the clade containing Rhinolophidae and Rhinopomatidae, and again independently in the lineages leading to Taphozous melanopogon and Mystacina tuberculata. Our reconstructions revealed that thermoregulatory regimes never transitioned directly from hibernation to homeothermy, or the reverse, in the evolutionary history of bats. This, in addition to recent evidence that heterothermy is best described along a continuum, suggests that thermoregulatory regimes in mammals are best represented as an ordered continuous trait (homeothermy ← → daily torpor ← → hibernation) rather than as the three discrete regimes that evolve in an unordered fashion. These results have important implications for methodological approaches in future physiological and evolutionary research.

Evolution for extreme living: variation in mitochondrial cytochrome c oxidase genes correlated with elevation in pikas (genus Ochotona)

The genus Ochotona (pikas) is a clade of cold-tolerant lagomorphs that includes many high-elevation species. Pikas offer a unique opportunity to study adaptations and potential limitations of an ecologically important mammal to high-elevation hypoxia. We analyzed the evolution of 3 mitochondrial genes encoding the catalytic core of cytochrome c oxidase (COX) in 10 pika species occupying elevations from sea level to 5000 m. COX is an enzyme highly reliant on oxygen and essential for cell function. One amino acid property, the equilibrium constant (ionization of COOH), was found to be under selection in the overall protein complex. We observed a strong relationship between the net value change in this property and the elevation each species occupies, with higher-elevation species having potentially more efficient proteins. We also found evidence of selection in low-elevation species for potentially less efficient COX, perhaps trading efficiency for heat production in the absence of hypoxia. Our results suggest that different pika species may have evolved elevation-specific COX proteins, specialization that may indicate limitations in their ability to shift their elevational ranges in response to future climate change.

Evolutionary Genetics of Hypoxia and Cold Tolerance in Mammals

Low oxygen and fluctuant ambient temperature pose serious challenges to mammalian survival. Physiological adaptations in mammals to hypoxia and low temperatures have been intensively investigated, yet their underlying molecular mechanisms need further exploration. Independent invasions of high-altitude plateaus, subterranean burrows and marine environments by different mammals provide opportunities to conduct such analyses. Here, we focused on six genes in the hypoxia inducible factor (HIF) pathway and two non-shivering thermogenesis (NST)-related genes [PPAR co-activator 1 (PGC-1) and uncoupling protein 1 (UCP1)] in representative species of pikas and other mammals to understand whether these loci were targeted by natural selection during independent invasions to conditions characterized by hypoxia and temperature fluctuations by high-altitude, subterranean and marine mammals. Our analyses revealed pervasive positive selection signals in the HIF pathway genes of mammals occupying high-altitude, subterranean and aquatic ecosystems; however, the mechanisms underlying their independent adaptations to hypoxic environments varied by taxa, since different genes were positively selected in each taxon and expression levels of individual genes varied among species. Additionally, parallel amino acid substitutions were also detected in hypoxia-tolerant mammals, indicating that convergent evolution may play a role in their independent adaptations to hypoxic environments. However, divergent evolutionary histories of NST-related genes were noted, since significant evidence of positive selection was observed in PGC-1 and UCP1 in high-altitude species and subterranean rodents; however, UCP1 may have already lost its function in diving cetaceans, which may be related to the thick blubber layer of adipose and connective tissue in these mammals.

Genome-wide identification, evolution, and expression analysis of the KT/HAK/KUP family in pear

The K⁺ transporter/high-affinity K⁺/K⁺ uptake (KT/HAK/KUP) family, as one of the largest K⁺ transporter families in higher plants, plays an essential role in plant growth, mineral element absorption, salt stress tolerance, and other physiological processes. However, little is known about this family in pear (Pyrus). Here, we identified 20 K⁺ transporter genes in pear (P. bretschneideri) using genome-wide analysis. Their gene structure, chromosomal distribution, conserved motifs, phylogenetics, duplication events, and expression patterns were also examined. The results of phylogenetic analysis showed that PbrKT/HAK/KUP genes were clustered into three major groups (Groups I-III). Among the 20 PbrKT/HAK/KUP genes, 18 were mapped to nine chromosomes and two to scaffolds. Four WGD/segmental gene pairs were identified, indicating that WGD/segmental duplication may have contributed to the expansion of the KT/HAK/KUP family in pear. Among the four pairs of WGD/segmentally duplicated genes, both members of three pairs had been subjected to purifying selection, whereas the fourth pair had been subjected to positive selection. Furthermore, phenotypic experiments showed that the growth of pear seedlings was affected by potassium deficiency treatment. Expression patterns of 20 PbrKT/HAK/KUP genes in roots were further assayed with qRT-PCR. PbrHAK1 and PbrHAK12/16 were significantly expressed in response to K⁺ deficiency, suggesting that these genes are crucial for K⁺ uptake in pear, especially under the condition of K⁺ starvation. Our results provide a foundation for further study on the function of KT/HAK/KUP genes in pear.

Diet simplification selects for high gut microbial diversity and strong fermenting ability in high-altitude pikas

The gut microbiota in mammals plays a key role in host metabolism and adaptation. However, relatively little is known regarding to how the animals adapts to extreme environments through regulating gut microbial diversity and function. Here, we investigated the diet, gut microbiota, short-chain fatty acid (SCFA) profiles, and cellulolytic activity from two common pika (Ochotona spp.) species in China, including Plateau pika (Ochotona curzoniae) from the Qinghai-Tibet Plateau and Daurian pika (Ochotona daurica) from the Inner Mongolia Grassland. Despite a partial diet overlap, Plateau pikas harbored lower diet diversity than Daurian pikas. Some bacteria (e.g., Prevotella and Ruminococcus) associated with fiber degradation were enriched in Plateau pikas. They harbored higher gut microbial diversity, total SCFA concentration, and cellulolytic activity than Daurian pikas. Interestingly, cellulolytic activity was positively correlated with the gut microbial diversity and SCFAs. Gut microbial communities and SCFA profiles were segregated structurally between host species. PICRUSt metagenome predictions demonstrated that microbial genes involved in carbohydrate metabolism and energy metabolism were overrepresented in the gut microbiota of Plateau pikas. Our results demonstrate that Plateau pikas harbor a stronger fermenting ability for the plant-based diet than Daurian pikas via gut microbial fermentation. The enhanced ability for utilization of plant-based diets in Plateau pikas may be partly a kind of microbiota adaptation for more energy requirements in cold and hypoxic high-altitude environments.

Gut region influences the diversity and interactions of bacterial communities in pikas (Ochotona curzoniae and Ochotona daurica)

The mammalian microbial communities in the gastrointestinal tract (GIT) play important roles in host nutrition and health. However, we still lack an understanding of how these communities are organized across GIT in natural environments. Here, using 16S rRNA gene sequencing, we analyzed the bacterial community diversity, network interactions and ecosystem stability across five gut regions (mouth, stomach, small intestine, cecum and colon) emanating from two common pika species in China, including Plateau pikas (Ochotona curzoniae) inhabiting high-altitude regions, as well as Daurian pikas (O. daurica) occupying low-altitude areas. The relative abundances of dominant Bacteroidetes and Firmicutes exhibited an increasing trend from mouth to colon. Cecum and colon harbored higher bacterial diversity compared with other anatomical regions. Gut region significantly influenced the structure of bacterial communities in the GIT. Network analysis indicated that topological features showed marked variations among gut regions. Interestingly, the ecosystem stability of bacterial communities increased gradually from mouth to colon. Our results suggest that gut region influences the diversity, structure and network interactions of bacterial communities in pikas. For hindgut-fermenting herbivorous mammals, relatively higher bacterial diversity and ecosystem stability in the cecum may provide a favorable condition for the fermentation of indigestible plant polysaccharides.

Molecular cloning, characterization and evolutionary analysis of leptin gene in Chinese giant salamander, Andrias davidianus

Leptin is an important hormone possessing diverse physiological roles in mammals and teleosts. However, it has been characterized only in a few amphibian species, and its evolutions are still under debate. Here, the full length of the leptin (Adlep) cDNA of Chinese giant salamander (Andrias davidianus), an early diverging amphibian species, is characterized and according to the results of the primary sequence analysis, tertiary structure reconstruction and phylogenetic analysis is confirmed to be an ortholog of mammalian leptin. An intron was identified between the coding exons of A. davidianus leptin, which indicated that the leptin is present in the salamander genome and contains a conserved gene structure in vertebrates. Adlep is widely distributed but expression levels vary among different tissues, with highest expression levels in the muscle. Additionally, the leptin receptor and other genes were mapped to three known leptin signaling pathways, suggesting that the leptin signaling pathways are present in A. davidianus. Phylogenetic topology of leptins are consistent with the generally accepted evolutionary relationships of vertebrates, and multiple leptin members found in teleosts seem to be obtained through a Cluopeocephala-specific gene duplication event. Our results will lay a foundation for further investigations into the physiological roles of leptin in A. davidianus.

Beyond thermoregulation: metabolic function of cetacean blubber in migrating bowhead and beluga whales

The processes of lipid deposition and utilization, via the gene leptin (Lep), are poorly understood in taxa with varying degrees of adipose storage. This study examines how these systems may have adapted in marine aquatic environments inhabited by cetaceans. Bowhead (Balaena mysticetus) and beluga whales (Delphinapterus leucas) are ideal study animals-they possess large subcutaneous adipose stores (blubber) and undergo bi-annual migrations concurrent with variations in food availability. To answer long-standing questions regarding how (or if) energy and lipid utilization adapted to aquatic stressors, we quantified variations in gene transcripts critical to lipid metabolism related to season, age, and blubber depth. We predicted leptin tertiary structure conservation and assessed inter-specific variations in Lep transcript numbers between bowheads and other mammals. Our study is the first to identify seasonal and age-related variations in Lep and lipolysis in these cetaceans. While Lep transcripts and protein oscillate with season in adult bowheads reminiscent of hibernating mammals, transcript levels reach up to 10 times higher in bowheads than any other mammal. Data from immature bowheads are consistent with the hypothesis that short baleen inhibits efficient feeding. Lipolysis transcripts also indicate young Fall bowheads and those sampled during Spring months limit energy utilization. These novel data from rarely examined species expand the existing knowledge and offer unique insight into how the regulation of Lep and lipolysis has adapted to permit seasonal deposition and maintain vital blubber stores.

Pika Gut May Select for Rare but Diverse Environmental Bacteria

The composition of the mammalian gut bacterial communities can be influenced by the introduction of environmental bacteria in their respective habitats. However, there are no extensive studies examining the interactions between environmental bacteriome and gut bacteriome in wild mammals. Here, we explored the relationship between the gut bacterial communities of pika (Ochotona spp.) and the related environmental bacteria across host species and altitudinal sites using 16S rRNA gene sequencing. Plateau pikas (O. curzoniae) and Daurian pikas (O. daurica) were sampled at five different sites, and plant and soil samples were collected at each site as well. Our data indicated that Plateau pikas and Daurian pikas had distinct bacterial communities. The pika, plant and soil bacterial communities were also distinct. Very little overlap occurred in the pika core bacteria and the most abundant environmental bacteria. The shared OTUs between pikas and environments were present in the environment at relatively low abundance, whereas they were affiliated with diverse bacterial taxa. These results suggested that the pika gut may mainly select for low-abundance but diverse environmental bacteria in a host species-specific manner.

Effect of Hypoxia on Ldh-c Expression in Somatic Cells of Plateau Pika

Sperm specific lactate dehydrogenases (LDH-C₄) is a lactate dehydrogenase that catalyzes the conversion of pyruvate to lactate. In mammals, Ldh-c was originally thought to be expressed only in testes and spermatozoa. Plateau pika (Ochotona curzoniae), which belongs to the genus Ochotona of the Ochotonidea family, is a hypoxia-tolerant mammal living 3000–5000 m above sea level on the Qinghai-Tibet Plateau, an environment which is strongly hypoxic. Ldh-c is expressed not only in testes and sperm, but also in the somatic tissues of plateau pika. To reveal the effect of hypoxia on pika Ldh-c expression, we investigated the mRNA and protein level of Ldh-c as well as the biochemical index of anaerobic glycolysis in pika somatic tissues at the altitudes of 2200 m, 3200 m and 3900 m. Our results showed that mRNA and protein expression levels of Ldh-c in the tissues of pika’s heart, liver, brain and skeletal muscle were increased significantly from 2200 m to 3200 m, but had no difference from 3200 m to 3900 m; the activities of LDH and the contents of lactate showed no difference from 2200 m to 3200 m, but were increased significantly from 3200 m to 3900 m. Hypoxia up-regulated and maintained the expression levels of Ldh-c in the pika somatic cells. Under the hypoxia condition, plateau pikas increased anaerobic glycolysis in somatic cells by LDH-C₄, and that may have reduced their dependence on oxygen and enhanced their adaptation to the hypoxic environment.

Enzymatic Kinetic Properties of the Lactate Dehydrogenase Isoenzyme C₄ of the Plateau Pika (Ochotona curzoniae)

Testis-specific lactate dehydrogenase (LDH-C₄) is one of the lactate dehydrogenase (LDH) isozymes that catalyze the terminal reaction of pyruvate to lactate in the glycolytic pathway. LDH-C₄ in mammals was previously thought to be expressed only in spermatozoa and testis and not in other tissues. Plateau pika (Ochotona curzoniae) belongs to the genus Ochotona of the Ochotonidea family. It is a hypoxia-tolerant species living in remote mountain areas at altitudes of 3000-5000 m above sea level on the Qinghai-Tibet Plateau. Surprisingly, Ldh-c is expressed not only in its testis and sperm, but also in somatic tissues of plateau pika. To shed light on the function of LDH-C₄ in somatic cells, Ldh-a, Ldh-b, and Ldh-c of plateau pika were subcloned into bacterial expression vectors. The pure enzymes of Lactate Dehydrogenase A₄ (LDH-A₄), Lactate Dehydrogenase B₄ (LDH-B₄), and LDH-C₄ were prepared by a series of expression and purification processes, and the three enzymes were identified by the method of sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and native polyacrylamide gel electrophoresis (PAGE). The enzymatic kinetics properties of these enzymes were studied by Lineweaver-Burk double-reciprocal plots. The results showed the Michaelis constant (Km) of LDH-C₄ for pyruvate and lactate was 0.052 and 4.934 mmol/L, respectively, with an approximate 90 times higher affinity of LDH-C₄ for pyruvate than for lactate. At relatively high concentrations of lactate, the inhibition constant (Ki) of the LDH isoenzymes varied: LDH-A₄ (Ki = 26.900 mmol/L), LDH-B₄ (Ki = 23.800 mmol/L), and LDH-C₄ (Ki = 65.500 mmol/L). These data suggest that inhibition of lactate by LDH-A₄ and LDH-B₄ were stronger than LDH-C₄. In light of the enzymatic kinetics properties, we suggest that the plateau pika can reduce reliance on oxygen supply and enhance its adaptation to the hypoxic environments due to increased anaerobic glycolysis by LDH-C₄.

Divergence of the bZIP Gene Family in Strawberry, Peach, and Apple Suggests Multiple Modes of Gene Evolution after Duplication

The basic leucine zipper (bZIP) transcription factors are the most diverse members of dimerizing transcription factors. In the present study, 50, 116, and 47 bZIP genes were identified in Malus domestica (apple), Prunus persica (peach), and Fragaria vesca (strawberry), respectively. Species-specific duplication was the main contributor to the large number of bZIPs observed in apple. After WGD in apple genome, orthologous bZIP genes corresponding to strawberry on duplicated regions in apple genome were retained. However, in peach ancestor, these syntenic regions were quickly lost or deleted. Maybe the positive selection contributed to the expansion of clade S to adapt to the development and environment stresses. In addition, purifying selection was mainly responsible for bZIP sequence-specific DNA binding. The analysis of orthologous pairs between chromosomes indicates that these orthologs derived from one gene duplication located on one of the nine ancient chromosomes in the Rosaceae. The comparative analysis of bZIP genes in three species provides information on the evolutionary fate of bZIP genes in apple and peach after they diverged from strawberry.

Genetic subdivision and candidate genes under selection in North American grey wolves

Previous genetic studies of the highly mobile grey wolf (Canis lupus) found population structure that coincides with habitat and phenotype differences. We hypothesized that these ecologically distinct populations (ecotypes) should exhibit signatures of selection in genes related to morphology, coat colour and metabolism. To test these predictions, we quantified population structure related to habitat using a genotyping array to assess variation in 42 036 single-nucleotide polymorphisms (SNPs) in 111 North American grey wolves. Using these SNP data and individual-level measurements of 12 environmental variables, we identified six ecotypes: West Forest, Boreal Forest, Arctic, High Arctic, British Columbia and Atlantic Forest. Next, we explored signals of selection across these wolf ecotypes through the use of three complementary methods to detect selection: FST /haplotype homozygosity bivariate percentilae, bayescan, and environmentally correlated directional selection with bayenv. Across all methods, we found consistent signals of selection on genes related to morphology, coat coloration, metabolism, as predicted, as well as vision and hearing. In several high-ranking candidate genes, including LEPR, TYR and SLC14A2, we found variation in allele frequencies that follow environmental changes in temperature and precipitation, a result that is consistent with local adaptation rather than genetic drift. Our findings show that local adaptation can occur despite gene flow in a highly mobile species and can be detected through a moderately dense genomic scan. These patterns of local adaptation revealed by SNP genotyping likely reflect high fidelity to natal habitats of dispersing wolves, strong ecological divergence among habitats, and moderate levels of linkage in the wolf genome.

Evidence of adaptive evolution of alpine pheasants to high-altitude environment from mitogenomic perspective

Adaptive evolutions to high-altitude adaptation have been intensively studied in mammals. However, considering the additional vertebrate groups, new perception regarding selection challenged by high-altitude stress on mitochondrial genome can be gained. To test this hypothesis, we compiled and analyzed the mitochondrial genomes of 5 alpine pheasants and 12 low-altitude species in Phasianidae. The results that evolutionary rates of ATP6 and ND6 showing significant fluctuation among branches when involved with five alpine pheasants revealed both genes might have implications with adapting to highland environment. The radical physico-chemical property changes identified by the modified MM01 model, including composition (C) and equilibrium constant (ionization of COOH) (Pk') in ATP6 and beta-structure tendencies (Pβ), Pk', and long-range non-bonded energy (El) in ND6, suggested that minor overall adjustments in size, protein conformation and relative orientation of reaction interfaces have been optimized to provide the ideal environments for electron transfer, proton translocation and generation of adenosine triphosphate (ATP). Additionally, three unique substitution sites were identified under selection in ND6, which could be potentially important adaptive changes contributing to cellular energy production. Our findings suggested that adaptive evolution may occur in alpine pheasants, which are an important complement to the knowledge of genetic mechanisms against the high-altitude environment in non-mammal animals.

Testis-Specific Lactate Dehydrogenase (LDH-C4) in Skeletal Muscle Enhances a Pika's Sprint-Running Capacity in Hypoxic Environment

LDH-C4 is a lactate dehydrogenase that catalyzes the conversion of pyruvate to lactate. In mammals, ldh-c was originally thought to be expressed only in testis and spermatozoa. Plateau pika (Ochotona curzoniae), which belongs to the genus Ochotona of the Ochotonidea family, is a hypoxia tolerant mammal living 3000-5000 m above sea level on the Qinghai-Tibet Plateau, an environment which is strongly hypoxic. Ldh-c is expressed not only in testis and sperm but also in somatic tissues of plateau pika. In this study, the effects of N-propyl oxamate and N-isopropyl oxamate on LDH isozyme kinetics were compared to screens for a selective inhibitor of LDH-C4. To reveal the role and physiological mechanism of LDH-C4 in skeletal muscle of plateau pika, we investigated the effect of N-isopropyl oxamate on the pika exercise tolerance as well as the physiological mechanism. Our results show that Ki of N-propyl oxamate and N-isopropyl oxamate for LDH-A4, LDH-B4, and LDH-C4 were 0.094 mmol/L and 0.462 mmol/L, 0.119 mmol/L and 0.248 mmol/L, and 0.015 mmol/L and 0.013 mmol/L, respectively. N-isopropyl oxamate is a powerful selective inhibitor of plateau pika LDH-C4. In our exercise tolerance experiment, groups treated with inhibitors had significantly lower swimming times than the uninhibited control group. The inhibition rates of LDH, LD, and ATP were 37.12%, 66.27%, and 32.42%, respectively. Our results suggested that ldh-c is expressed in the skeletal muscle of plateau pika, and at least 32.42% of ATP in the skeletal muscle is catalyzed by LDH-C4 by anaerobic glycolysis. This suggests that pika has reduced dependence on oxygen and enhanced adaptation to hypoxic environment due to increased anaerobic glycolysis by LDH-C4 in skeletal muscle. LDH-C4 in plateau pika plays the crucial role in anaerobic glycolysis and generates ATP rapidly since this is the role of LDH-A4 in most species on plain land, which provide evidence that the native humans and animals in Qinghai-Tibet plateau can adapt to the hypoxia environment.

Intermittent cold exposure results in visceral adipose tissue "browning" in the plateau pika (Ochotona curzoniae)

The plateau pika has developed tolerance to cold and hypoxia in order to adapt to living in the extreme environment of the Qinghai-Tibetan Plateau. One mammalian mechanism for cold adaptation is thermogenesis by brown adipose tissue (BAT), but the degree to which pika exploits this mechanism or how it may be modified by the additional stresses of high altitude is not known. Intermittent Cold Exposure (ICE) is an approachable method to study cold adaptation in rodents. To investigate the role of adipose tissue in the adaptation of pika to cold temperatures, we have studied pika during ICE. We find that pika kept in warm temperatures has little classical brown fat, but "browning" of white adipose tissues is observed rapidly upon cold exposure. This is demonstrated by the increased expression of several markers of brown fat differentiation including uncoupling protein 1 (UCP-1). Surprisingly, this occurs mainly in visceral rather than epididymal adipose tissue. In addition, ICE increases the expression of several general adipose differentiation markers at both the mRNA and protein levels. These substantial changes in the distribution of fat are accomplished without changes in weight or blood levels of glucose and triglycerides, suggesting that the adaptable changes are coordinated and self-compensated. Together, our results demonstrate that ICE promotes recruitment of BAT in pika, and unlike small mammals in at lower altitudes, pika can activate visceral WAT to adapt to cold stress without major changes overall energy balance.

VEGF189 expression is highly related to adaptation of the plateau pika (Ochotona curzoniae) inhabiting high altitudes

The plateau pika (Ochotona curzonia) has adapted to high-altitude hypoxia during evolution. Higher microvessel density in specific tissues and a blunted hypoxic pulmonary vasoconstriction response are the critical components of this adaptation. VEGF, vascular endothelial growth factor, has proved to be a key regulator of angiogenesis in response to tissue hypoxia and to play an important role in vascular vasodilation. However, the role of VEGF in adaptation to high-altitude hypoxia in the plateau pika remains unknown. In this study, we cloned cDNAs for VEGF165 and VEGF189 and examined their expression in pikas inhabiting altitudes of 3200 and 4750 m. Phylogenetic analysis reveals that pika VEGF165 and VEGF189 are evolutionarily conserved. Real-time PCR analysis demonstrates that VEGF165 and VEGF189 display tissue and altitude-specific expression patterns. Interestingly, we found that the levels of VEGF189 mRNA are significantly higher than those of VEGF165 in the brain and muscle tissues of the pika, which is different from what was previously observed in sea-level mammals. VEGF189 mRNA levels in brain, muscle, and lung of the pika increased with increased habitat altitude, whereas VEGF165 shows less change. Our study suggests an important role for VEGF189 in adaptation to hypoxia by the plateau pika in the high-altitude environment.

The complete mitochondrial genome sequence of the Tibetan red fox (Vulpes vulpes montana)

In this study, the complete mitochondrial genome of the Tibetan red fox (Vulpes Vulpes montana) was sequenced for the first time using blood samples obtained from a wild female red fox captured from Lhasa in Tibet, China. Qinghai--Tibet Plateau is the highest plateau in the world with an average elevation above 3500 m. Sequence analysis showed it contains 12S rRNA gene, 16S rRNA gene, 22 tRNA genes, 13 protein-coding genes and 1 control region (CR). The variable tandem repeats in CR is the main reason of the length variability of mitochondrial genome among canide animals.

The complete mitochondrial genome sequence of the Tibetan wolf (Canis lupus laniger)

In this study, the complete mitochondrial genome of the Tibetan wolf (Canis lupus laniger) was sequenced using blood samples obtained from a wild female Tibetan wolf captured from Lhasa in Tibet, China. Qinghai-Tibet Plateau, with an average elevation above 3500 m, is the highest plateau in the world. Sequence analysis showed that its structure is in accordance with other Canidae species, but GTG is used as the start codon in ND4L gene which is different from many canide animals.

Leptin gene in rabbit: cloning and expression in mammary epithelial cells during pregnancy and lactation

Leptin is known as a cytokine mostly produced by fat cells and implicated in regulation of energy metabolism and food intake but has also been shown to be involved in many physiological mechanisms such as tissue metabolism and cell differentiation and proliferation. In particular, leptin influences the development of mammary gland. Although leptin expression in mammary gland has been studied in several species, no data are available in the rabbit. Leptin transcripts in this species have been described as being encoded by only two exons rather than three as in other species. Our focus was to clone and sequence the rabbit leptin cDNA and to prepare the recombinant biologically active protein for validation of the proper sequence and then to describe leptin expression in rabbit mammary gland during different stages of pregnancy and lactation. The leptin sequence obtained was compared with those of other species, and genome alignment demonstrated that the rabbit leptin gene is also encoded by three exons. Additionally, we analyzed the expression of leptin during pregnancy and lactation. Leptin mRNA was weakly expressed throughout pregnancy, whereas mRNA levels were higher during lactation, with a significant increase between days 3 and 16. Leptin transcripts and protein were localized in luminal epithelial cells, thus indicating that leptin synthesis occurs in this compartment. Therefore, mammary synthesized leptin may constitute a major regulator of mammary gland development by acting locally as an autocrine and/or paracrine factor. Furthermore, our results support the possible physiological role of leptin in newborns through consumption of milk.

Identification of cda gene in bighead carp and its expression in response to microcystin-LR

Microcystin-LR (MCLR) is a widespread cyanotoxin, which can influence genes transcription and cause nucleic acid damage in different organisms. To identify MCLR induced transcriptionally changed hepatic genes in bighead carp by subtractive suppression hybridization, we obtained the cDNA fragment of cda. Then we cloned its full-length cDNA, which encodes a cytidine deaminase (CDA). 3D structure prediction showed that the 3D structure and amino acid residues related to function sites of bighead carp CDA were highly conserved. Bighead carp CDA shared high identities with other CDA sequences, and evolved closely to non-mammalian CDAs. Bighead carp expressed cda in all tested tissues under normal situation, and changed its expression profile in a time inversely dependent and dose dependent manner to MCLR, so as to protect itself from MCLR induced toxic damage. These indicated that cda might be involved in anti-MCLR response, especially in the regulation of cytidine and dexocytidine metabolism pathway.

Phocid seal leptin: tertiary structure and hydrophobic receptor binding site preservation during distinct leptin gene evolution

The cytokine hormone leptin is a key signalling molecule in many pathways that control physiological functions. Although leptin demonstrates structural conservation in mammals, there is evidence of positive selection in primates, lagomorphs and chiropterans. We previously reported that the leptin genes of the grey and harbour seals (phocids) have significantly diverged from other mammals. Therefore we further investigated the diversification of leptin in phocids, other marine mammals and terrestrial taxa by sequencing the leptin genes of representative species. Phylogenetic reconstruction revealed that leptin diversification was pronounced within the phocid seals with a high dN/dS ratio of 2.8, indicating positive selection. We found significant evidence of positive selection along the branch leading to the phocids, within the phocid clade, but not over the dataset as a whole. Structural predictions indicate that the individual residues under selection are away from the leptin receptor (LEPR) binding site. Predictions of the surface electrostatic potential indicate that phocid seal leptin is notably different to other mammalian leptins, including the otariids. Cloning the grey seal leptin binding domain of LEPR confirmed that this was structurally conserved. These data, viewed in toto, support a hypothesis that phocid leptin divergence is unlikely to have arisen by random mutation. Based upon these phylogenetic and structural assessments, and considering the comparative physiology and varying life histories among species, we postulate that the unique phocid diving behaviour has produced this selection pressure. The Phocidae includes some of the deepest diving species, yet have the least modified lung structure to cope with pressure and volume changes experienced at depth. Therefore, greater surfactant production is required to facilitate rapid lung re-inflation upon surfacing, while maintaining patent airways. We suggest that this additional surfactant requirement is met by the leptin pulmonary surfactant production pathway which normally appears only to function in the mammalian foetus.

Adaptive evolution in the glucose transporter 4 gene Slc2a4 in Old World fruit bats (family: Pteropodidae)

Frugivorous and nectarivorous bats are able to ingest large quantities of sugar in a short time span while avoiding the potentially adverse side-effects of elevated blood glucose. The glucose transporter 4 protein (GLUT4) encoded by the Slc2a4 gene plays a critical role in transmembrane skeletal muscle glucose uptake and thus glucose homeostasis. To test whether the Slc2a4 gene has undergone adaptive evolution in bats with carbohydrate-rich diets in relation to their insect-eating sister taxa, we sequenced the coding region of the Slc2a4 gene in a number of bat species, including four Old World fruit bats (Pteropodidae) and three New World fruit bats (Phyllostomidae). Our molecular evolutionary analyses revealed evidence that Slc2a4 has undergone a change in selection pressure in Old World fruit bats with 11 amino acid substitutions detected on the ancestral branch, whereas, no positive selection was detected in the New World fruit bats. We noted that in the former group, amino acid replacements were biased towards either Serine or Isoleucine, and, of the 11 changes, six were specific to Old World fruit bats (A133S, A164S, V377F, V386I, V441I and G459S). Our study presents preliminary evidence that the Slc2a4 gene has undergone adaptive changes in Old World fruit bats in relation to their ability to meet the demands of a high sugar diet.

Adaptive evolution of Leptin in heterothermic bats

Heterothermy (hibernation and daily torpor) is a key strategy that animals use to survive in harsh conditions and is widely employed by bats, which are found in diverse habitats and climates. Bats comprise more than 20% of all mammals and although heterothermy occurs in divergent lineages of bats, suggesting it might be an ancestral condition, its evolutionary history is complicated by complex phylogeographic patterns. Here, we use Leptin, which regulates lipid metabolism and is crucial for thermogenesis of hibernators, as molecular marker and combine physiological, molecular and biochemical analyses to explore the possible evolutionary history of heterothermy in bat. The two tropical fruit bats examined here were homeothermic; in contrast, the two tropical insectivorous bats were clearly heterothermic. Molecular evolutionary analyses of the Leptin gene revealed positive selection in the ancestors of all bats, which was maintained or further enhanced the lineages comprising mostly heterothermic species. In contrast, we found evidence of relaxed selection in homeothermic species. Biochemical assays of bat Leptin on the activity on adipocyte degradation revealed that Leptin in heterothermic bats was more lipolytic than in homeothermic bats. This shows that evolutionary sequence changes in this protein are indeed functional and support the interpretation of our physiological results and the molecular evolutionary analyses. Our combined data strongly support the hypothesis that heterothermy is the ancestral state of bats and that this involved adaptive changes in Leptin. Subsequent loss of heterothermy in some tropical lineages of bats likely was associated with range and dietary shifts.

Evidence for positive selection on the leptin gene in Cetacea and Pinnipedia

The leptin gene has received intensive attention and scientific investigation for its importance in energy homeostasis and reproductive regulation in mammals. Furthermore, study of the leptin gene is of crucial importance for public health, particularly for its role in obesity, as well as for other numerous physiological roles that it plays in mammals. In the present work, we report the identification of novel leptin genes in 4 species of Cetacea, and a comparison with 55 publicly available leptin sequences from mammalian genome assemblies and previous studies. Our study provides evidence for positive selection in the suborder Odontoceti (toothed whales) of the Cetacea and the family Phocidae (earless seals) of the Pinnipedia. We also detected positive selection in several leptin gene residues in these two lineages. To test whether leptin and its receptor evolved in a coordinated manner, we analyzed 24 leptin receptor gene (LPR) sequences from available mammalian genome assemblies and other published data. Unlike the case of leptin, our analyses did not find evidence of positive selection for LPR across the Cetacea and Pinnipedia lineages. In line with this, positively selected sites identified in the leptin genes of these two lineages were located outside of leptin receptor binding sites, which at least partially explains why co-evolution of leptin and its receptor was not observed in the present study. Our study provides interesting insights into current understanding of the evolution of mammalian leptin genes in response to selective pressures from life in an aquatic environment, and leads to a hypothesis that new tissue specificity or novel physiologic functions of leptin genes may have arisen in both odontocetes and phocids. Additional data from other species encompassing varying life histories and functional tests of the adaptive role of the amino acid changes identified in this study will help determine the factors that promote the adaptive evolution of the leptin genes in marine mammals.

Nonsynonymous natural genetic polymorphisms in the bovine leptin gene affect biochemical and biological characteristics of the mature hormone

Leptin (LEP) is a cytokine-like hormone proven to be involved in diverse biological processes. In livestock, it regulates feed intake, BW homeostasis, and energy balance, among other traits. Natural nonsynonymous genetic polymorphisms in the ovine leptin (oLEP) alter the biochemical and physiological characteristics of its gene products. Here we studied in vitro and in vivo the biochemical and physiological characteristics of recombinant hormones representing the oLEP and bovine leptin (bLEP) reference sequences of wild-type (WT) leptins (GenBank accession No. U84247 and U50365, respectively), oLEP and bLEP recombinant muteins carrying the R4C mutation, and oLEP recombinant hormones carrying the A59V and Q62R mutations, which were detected in bLEP. All proteins were purified to homogeneity as monomers and formed 1:1 molar ratio complexes with the chicken leptin-binding domain (LBD). Surface plasmon resonance experiments revealed that all protein variants exhibit reduced (P < 0.05) affinity to chicken (ch) and human (h) LBD compared with the WT oLEP and bLEP recombinant proteins. The ovine and bovine R4C muteins exhibited significantly (P < 0.05) greater induction of cell proliferation in a Baf/3 cell line bioassay, despite lower affinity toward both hLBD and chLBD. Intra-third cerebral ventricle infusion of oLEP and its 3 muteins in sheep resulted in reduced feed intake. However, the 3 tested muteins had a decreased (P < 0.05) inhibitory effect than the WT LEP. It was concluded that natural genetic polymorphisms in the bLEP are associated with variation in the biochemical and physiological properties of the protein.

Biochemical and in vitro biological significance of natural sequence variation in the ovine leptin gene

The hormone leptin is involved in diverse biological processes, including regulation of food intake, body-weight homeostasis and energy balance. Sequence variation in the bovine leptin gene has been found to be associated with variations in carcass fat content and average daily gain, as well as in milk yield, milk somatic cell count and several traits governing reproduction. We sequenced genomic DNA and cDNA samples of individuals from three divergent sheep breeds and revealed synonymous as well as novel non-synonymous allelic variation at the third exon of the ovine leptin gene (oLEP) as compared to the sequence published at Accession No. U84247 (reference sequence). In addition, two alternatively spliced oLEP transcripts were found in the abdominal fat tissue. The biochemical and the in vitro biological significance of the sequence variation in the oLEP was examined by generating recombinant oLEP-protein variants namely: p.Q28del, p.N78S, p.R84Q, p.P99Q, p.V123L and p.R138Q, carrying the corresponding sequence variation. Surface plasmon resonance experiments revealed, in most cases, reduced affinity of the oLEP protein variants examined, to human leptin-binding domain (hLBD), relative to the reference variant, being 0.75, 0.60, 0.60, 0.89, 0.92 and 1.03, respectively. In competitive binding assays between biotinylated oLEP and the recombinant leptin protein variants, p.N78S and p.R84Q variants exhibited the lowest affinity to hLBD (0.18 and 0.41, respectively) as compared to the reference hormone. We then tested the protein variants' ability to induce proliferation in Baf-3 cells stably expressing the long form of the human leptin receptor: significant differences in proliferative activity were only found for p.N78S (1.8-fold higher) and p.R138Q (4.2-fold lower) relative to the reference oLEP variant.

Functional evolution of leptin of Ochotona curzoniae in adaptive thermogenesis driven by cold environmental stress

Background

Environmental stress can accelerate the directional selection and evolutionary rate of specific stress-response proteins to bring about new or altered functions, enhancing an organism's fitness to challenging environments. Plateau pika (Ochotona curzoniae), an endemic and keystone species on Qinghai-Tibetan Plateau, is a high hypoxia and low temperature tolerant mammal with high resting metabolic rate and non-shivering thermogenesis to cope in this harsh plateau environment. Leptin is a key hormone related to how these animals regulate energy homeostasis. Previous molecular evolutionary analysis helped to generate the hypothesis that adaptive evolution of plateau pika leptin may be driven by cold stress.

Methodology/Principal Findings

To test the hypothesis, recombinant pika leptin was first purified. The thermogenic characteristics of C57BL/6J mice injected with pika leptin under warm (23±1°C) and cold (5±1°C) acclimation is investigated. Expression levels of genes regulating adaptive thermogenesis in brown adipose tissue and the hypothalamus are compared between pika leptin and human leptin treatment, suggesting that pika leptin has adaptively and functionally evolved. Our results show that pika leptin regulates energy homeostasis via reduced food intake and increased energy expenditure under both warm and cold conditions. Compared with human leptin, pika leptin demonstrates a superior induced capacity for adaptive thermogenesis, which is reflected in a more enhanced β-oxidation, mitochondrial biogenesis and heat production. Moreover, leptin treatment combined with cold stimulation has a significant synergistic effect on adaptive thermogenesis, more so than is observed with a single cold exposure or single leptin treatment.

Conclusions/Significance

These findings support the hypothesis that cold stress has driven the functional evolution of plateau pika leptin as an ecological adaptation to the Qinghai-Tibetan Plateau.

Leptin in teleost fishes: an argument for comparative study

All organisms face tradeoffs with regard to how limited energy resources should be invested. When is it most favorable to grow, to reproduce, how much lipid should be allocated to storage in preparation for a period of limited resources (e.g., winter), instead of being used for growth or maturation? These are a few of the high consequence fitness "decisions" that represent the balance between energy acquisition and allocation. Indeed, for animals to make favorable decisions about when to grow, eat, or reproduce, they must integrate signals among the systems responsible for energy acquisition, storage, and demand. We make the argument that leptin signaling is a likely candidate for an integrating system. Great progress has been made understanding the leptin system in mammals, however our understanding in fishes has been hampered by difficulty in cloning fish orthologs of mammalian proteins and (we assert), underutilization of the comparative approach.

Mitogenomic analysis of Chinese snub-nosed monkeys: Evidence of positive selection in NADH dehydrogenase genes in high-altitude adaptation

Chinese snub-nosed monkeys belong to the genus Rhinopithecus and are limited in distribution to six isolated mountainous areas in the temperate regions of Central and Southwest China. Compared to the other members of the subfamily Colobinae (or leaf-eating monkeys), these endangered primates are unique in being adapted to a high altitude environment and display a remarkable ability to tolerate low temperatures and hypoxia. They thus offer an interesting organismal model of adaptation to extreme environmental stress. Mitochondria generate energy by oxidative phosphorylation (OXPHOS) and play important roles in oxygen usage and energy metabolism. We analyzed the mitochondrial genomes of two Chinese snub-nosed monkey species and eight other colobines in the first attempt to understand the genetic basis of high altitude adaptation in non-human primates. We found significant evidence of positive selection in one Chinese snub-nosed monkey, Rhinopithecus roxellana, which is suggestive of adaptive change related to high altitude and cold weather stress. In addition, our study identified two potentially important adaptive amino acid residues (533 and 3307) in the NADH2 and NADH6 genes, respectively. Surprisingly, no evidence for positive selection was found in Rhinopithecus bieti (the other Chinese snub-nosed monkey analyzed). This finding is intriguing, especially considering that R. bieti inhabits a higher altitudinal distribution than R. roxellana. We hypothesize that a different adaptive genetic basis to high altitude survival exists in R. bieti from those seen in other mammals, and that positive selection and functionally associated mutations in this species may be detected in nuclear genes related to energy and oxygen metabolism. More information on the structure, function, and evolution of mitochondrial and nuclear genomes in Chinese snub-nosed monkeys is required to reveal the molecular mechanisms that underlie adaptations to high altitude survival in non-human primates.

Evolution of leptin structure and function

Leptin, the protein product of the obese(ob or Lep) gene, is a hormone synthesized by adipocytes that signals available energy reserves to the brain, and thereby influences development, growth, metabolism and reproduction. In mammals, leptin functions as an adiposity signal: circulating leptin fluctuates in proportion to fat mass, and it acts on the hypothalamus to suppress food intake. Orthologs of mammalian Lep genes were recently isolated from several fish and two amphibian species, and here we report the identification of two Lep genes in a reptile, the lizard Anolis carolinensis. While vertebrate leptins show large divergence in their primary amino acid sequence, they form similar tertiary structures, and may have similar potencies when tested in vitro on heterologous leptin receptors (LepRs). Leptin binds to LepRs on the plasma membrane, activating several intracellular signaling pathways. Vertebrate LepRs signal via the Janus kinase (Jak) and signal transducer and activator of transcription (STAT) pathway. Three tyrosine residues located within the LepR cytoplasmic domain are phosphorylated by Jak2 and are required for activation of SH2-containing tyrosine phosphatase-2, STAT5 and STAT3 signaling. These tyrosines are conserved from fishes to mammals, demonstrating their critical role in signaling by the LepR. Leptin is anorexigenic in representatives of all vertebrate classes, suggesting that its role in energy balance is ancient and has been evolutionarily conserved. In addition to its integral role as a regulator of appetite and energy balance, leptin exerts pleiotropic actions in development, physiology and behavior.

Molecular cloning and gene expression analysis of the leptin receptor in the Chinese mitten crab Eriocheir sinensis

Background

Leptin is an adipocyte-derived hormone with multiple functions that regulates energy homeostasis and reproductive functions. Increased knowledge of leptin receptor function will enhance our understanding of the physiological roles of leptin in animals.

Methodology/Principal Findings

In the present study, a full-length leptin receptor (lepr) cDNA, consisting of 1,353 nucleotides, was cloned from Chinese mitten crab (Eriocheir sinensis) using rapid amplification of cDNA ends (RACE) following the identification of a single expressed sequence tag (EST) clone in a cDNA library. The lepr cDNA consisted of a 22-nucleotide 5′-untranslated region (5′ UTR), a 402-nucleotide open reading frame (ORF) and a 929-nucleotide 3′ UTR. Multiple sequence alignments revealed that Chinese mitten crab lepr shared a conserved vacuolar protein sorting 55 (Vps55) domain with other species. Chinese mitten crab lepr expression was determined in various tissues and at three different reproductive stages using quantitative real-time RT-PCR. Lepr expression was highest in the intestine, thoracic ganglia, gonad, and accessory gonad, moderate in hepatopancreas and cranial ganglia, and low in muscle, gill, heart, haemocytes, and stomach. Furthermore, lepr expression was significantly higher in the intestine, gonad and thoracic ganglia in immature crabs relative to precocious and mature crabs. In contrast, lepr expression was significantly lower in the hepatopancreas of immature crabs relative to mature crabs.

Conclusions/Significance

We are the first to identify the lepr gene and to determine its gene expression patterns in various tissues and at three different reproductive stages in Chinese mitten crab. Taken together, our results suggest that lepr may be involved in the nutritional regulation of metabolism and reproduction in Chinese mitten crabs.

Adaptive evolution and inherent tolerance to extreme thermal environments

BACKGROUND

When introduced to novel environments, the ability for a species to survive and rapidly proliferate corresponds with its adaptive potential. Of the many factors that can yield an environment inhospitable to foreign species, phenotypic response to variation in the thermal climate has been observed within a wide variety of species. Experimental evolution studies using bacteriophage model systems have been able to elucidate mutations, which may correspond with the ability of phage to survive modest increases/decreases in the temperature of their environment.

RESULTS

Phage PhiX174 was subjected to both elevated (50 degrees C) and extreme (70 degrees C+) temperatures for anywhere from a few hours to days. While no decline in the phage's fitness was detected when it was exposed to 50 degrees C for a few hours, more extreme temperatures significantly impaired the phage; isolates that survived these heat treatments included the acquisition of several mutations within structural genes. As was expected, long-term treatment of elevated and extreme temperatures, ranging from 50-75 degrees C, reduced the survival rate even more. Isolates which survived the initial treatment at 70 degrees C for 24 or 48 hours exhibited a significantly greater tolerance to subsequent heat treatments.

CONCLUSIONS

Using the model organism PhiX174, we have been able to study adaptive evolution on the molecular level under extreme thermal changes in the environment, which to-date had yet to be thoroughly examined. Under both acute and extended thermal selection, we were able to observe mutations that occurred in response to excessive external pressures independent of concurrently evolving hosts. Even though its host cannot tolerate extreme temperatures such as the ones tested here, this study confirms that PhiX174 is capable of survival.

Mitochondrial genome analysis of Ochotona curzoniae and implication of cytochrome c oxidase in hypoxic adaptation

Pikas originated in Asia and are small lagomorphs native to cold climates. The plateau pika, Ochotona curzoniae is a keystone species on the Qinghai-Tibet Plateau and an ideal animal model for hypoxic adaptation studies. Altered mitochondrial function, especially cytochrome c oxidase activity, is an important factor in modulation of energy generation and expenditure during cold and hypoxia adaptation. In this study, we determined the complete nucleotide sequence of the O. curzoniae mitochondrial genome. The plateau pika mitochondrial DNA is 17,131bp long and encodes the complete set of 37 proteins typical for vertebrates. Phylogenetic analysis based on concatenated heavy-strand encoded protein-coding genes revealed that pikas are closer to rabbit and hare than to rat. This suggests that rabbit or hare would be a good control animal for pikas in cold and hypoxia adaptation studies. Fifteen novel mitochondrial DNA-encoded amino acid changes were identified in the pikas, including three in the subunits of cytochrome c oxidase. These amino acid substitutions potentially function in modulation of mitochondrial complexes and electron transport efficiency during cold and hypoxia adaptation.