Divergent and parallel routes of biochemical adaptation in high-altitude passerine birds from the Qinghai-Tibet Plateau.

When different species experience similar selection pressures, the probability of evolving similar adaptive solutions may be influenced by legacies of evolutionary history, such as lineage-specific changes in genetic background. Here we test for adaptive convergence in hemoglobin (Hb) function among high-altitude passerine birds that are native to the Qinghai-Tibet Plateau, and we examine whether convergent increases in Hb-O2 affinity have a similar molecular basis in different species. We documented that high-altitude parid and aegithalid species from the Qinghai-Tibet Plateau have evolved derived increases in Hb-O2 affinity in comparison with their closest lowland relatives in East Asia. However, convergent increases in Hb-O2 affinity and convergence in underlying functional mechanisms were seldom attributable to the same amino acid substitutions in different species. Using ancestral protein resurrection and site-directed mutagenesis, we experimentally confirmed two cases in which parallel substitutions contributed to convergent increases in Hb-O2 affinity in codistributed high-altitude species. In one case involving the ground tit (Parus humilis) and gray-crested tit (Lophophanes dichrous), parallel amino acid replacements with affinity-enhancing effects were attributable to nonsynonymous substitutions at a CpG dinucleotide, suggesting a possible role for mutation bias in promoting recurrent changes at the same site. Overall, most altitude-related changes in Hb function were caused by divergent amino acid substitutions, and a select few were caused by parallel substitutions that produced similar phenotypic effects on the divergent genetic backgrounds of different species.

Evolution of beak morphology in the Ground Tit revealed by comparative transcriptomics.

BACKGROUND: Beak morphology exhibits considerable adaptive plasticity in birds, which results in highly varied or specialized forms in response to variations in ecology and life history. As the only parid species endemic to the Qinghai-Tibet Plateau, the Ground Tit (Parus humilis) has evolved a distinctly long and curved beak from other parids. An integration of morphometrics, phylogenetics, transcriptomics and embryology allows us to address the evolutionary and developmental mechanisms of the adaptive beak structure observed in the Ground Tit. RESULTS: A morphometric approach quantified that the Ground Tit has a comparatively longer and more decurved upper beaks than other parids. We estimated that the ancestor of the Ground Tit likely had a short straight upper beak similar to most current recognized parid species using an ancestral state reconstruction. This morphological specialization is considered an adaptation to its ground-oriented behavior on the high plateau. To identify genetic mechanisms behind this adaptive change, a comparative transcriptomic analysis was applied between the Ground Tit and its closely related species, the Great Tit (Parus major). We detected that 623 genes were significantly differentially expressed in embryonic upper beaks between the two species, 17 of which were functionally annotated to correlate with bone development and morphogenesis, although genes related to bone development were not found to undergo accelerated evolution in the Ground Tit. RT-qPCR validation confirmed differential expression of five out of eight genes that were selected from the 17 genes. Subsequent functional assays in chicken embryos demonstrated that two of these genes, FGF13 and ITGB3, may affect beak morphology by modulating levels of osteoblasts and osteoclasts. CONCLUSIONS: Our results provide preliminary evidence that development of the long decurved beak of the Ground Tit is likely regulated by transcriptional activities of multiple genes coordinating osteoblasts and osteoclasts. The integration of multiple approaches employed here sheds light on ecological and genetic mechanisms in the evolution of avian morphology.

The voltage-gated potassium channel subfamily KQT member 4 (KCNQ4) displays parallel evolution in echolocating bats.

Bats are the only mammals that use highly developed laryngeal echolocation, a sensory mechanism based on the ability to emit laryngeal sounds and interpret the returning echoes to identify objects. Although this capability allows bats to orientate and hunt in complete darkness, endowing them with great survival advantages, the genetic bases underlying the evolution of bat echolocation are still largely unknown. Echolocation requires high-frequency hearing that in mammals is largely dependent on somatic electromotility of outer hair cells. Then, understanding the molecular evolution of outer hair cell genes might help to unravel the evolutionary history of echolocation. In this work, we analyzed the molecular evolution of two key outer hair cell genes: the voltage-gated potassium channel gene KCNQ4 and CHRNA10, the gene encoding the α10 nicotinic acetylcholine receptor subunit. We reconstructed the phylogeny of bats based on KCNQ4 and CHRNA10 protein and nucleotide sequences. A phylogenetic tree built using KCNQ4 amino acid sequences showed that two paraphyletic clades of laryngeal echolocating bats grouped together, with eight shared substitutions among particular lineages. In addition, our analyses indicated that two of these parallel substitutions, M388I and P406S, were probably fixed under positive selection and could have had a strong functional impact on KCNQ4. Moreover, our results indicated that KCNQ4 evolved under positive selection in the ancestral lineage leading to mammals, suggesting that this gene might have been important for the evolution of mammalian hearing. On the other hand, we found that CHRNA10, a gene that evolved adaptively in the mammalian lineage, was under strong purifying selection in bats. Thus, the CHRNA10 amino acid tree did not show echolocating bat monophyly and reproduced the bat species tree. These results suggest that only a subset of hearing genes could underlie the evolution of echolocation. The present work continues to delineate the genetic bases of echolocation and ultrasonic hearing in bats.

Identification of Two Novel Variants of the DMD Gene in Chinese Families with Duchenne Muscular Dystrophy.

Background: Duchenne muscular dystrophy (DMD), an X-linked recessive neuromuscular disorder, is caused by pathogenic variants in the DMD gene encoding a large structural protein in muscle cells. Methods: Two probands, a 6-year old boy and a 1-month old infant, respectively, were clinically diagnosed with DMD based on elevated levels of creatine kinase and creatine kinase isoenzyme. CNVplex and whole exome sequencing (WES) were performed for causal variants, and Sanger sequencing was used for verification. Results: CNVplex found no large deletions or duplications in the DMD gene in both patients, but WES discovered a single-nucleotide deletion in exon 48 (NM_004006.2:c.6963del, p.Asp2322ThrfsTer16) in the proband of pedigree 1, and a nonsense mutation in exon 27 (NM_004006.2:c.3637A>T, p.K1213Ter) in the proband of pedigree 2. Conclusion: The results of our study expand the mutation spectrum of DMD and enrich our understanding of the clinical characteristics of DMD. Genetic counseling was provided for the two families involved in this study.

Molecular cloning and evolutionary analysis of GJB6 in mammals.

GJB6 plays a crucial role in hearing. In mammals, bats use ultrasonic echolocation for orientation and locating prey. To investigate the evolution of GJB6 in mammals, we cloned the full-length coding region of GJB6 from 16 species of bats and 4 other mammal species and compared them with orthologous sequences in 11 other mammals. The results show purifying selection on GJB6 in mammals, as well as in the bat lineage, which indicates an important role for GJB6 in mammal hearing. We also found one unique amino acid substitution shared by 16 species of bats and 10 shared by two species of artiodactyls. This positioned the artiodactyls at an abnormal location in the gene tree. In addition, the cytoplasmic loop and carboxy terminus were more variable than other domains in all the mammals. These results demonstrate that GJB6 is basically conserved in mammals but has undergone relatively rapid evolution in particular lineages and domains.

A comparative study of prenatal development in Miniopterus schreibersii fuliginosus, Hipposideros armiger and H. pratti.

BACKGROUND: Bats comprise the second largest order of mammals. However, there are far fewer morphological studies of post-implantation embryonic development than early embryonic development in bats. RESULTS: We studied three species of bats (Miniopterus schreibersii fuliginosus, Hipposideros armiger and H. pratti), representing the two suborders Yangochiroptera and Yinpterochiroptera. Using an established embryonic staging system, we identified the embryonic stages for M. schreibersii fuliginosus, H. armiger and H. pratti and described the morphological changes in each species, including the development of the complex and distinctive nose-leaves in H. armiger and H. pratti. Finally, we compared embryonic and fetal morphology of the three species in the present study with five other species for which information is available. CONCLUSION: As a whole, the organogenetic sequence of bat embryos is uniform and the embryos appear homoplastic before Stage 16. Morphological differentiation between species occurs mainly after embryonic Stage 16. Our study provides three new bat species for interspecific comparison of post-implantation embryonic development within the order Chiroptera and detailed data on the development of nose-leaves for bats in the superfamily Rhinolophoidea.

Digital gene expression tag profiling of bat digits provides robust candidates contributing to wing formation.

BACKGROUND: As the only truly flying mammals, bats use their unique wing - consisting of four elongated digits (digits II-V) connected by membranes - to power their flight. In addition to the elongated digits II-V, the forelimb contains one shorter digit (digit I) that is morphologically similar to the hindlimb digits. Here, we capitalized on the morphological variation among the bat forelimb digits to investigate the molecular mechanisms underlying digit elongation and wing formation. Using next generation sequencing technology, we performed digital gene expression tag profiling (DGE-tag profiling) of developing digits in a pooled sample of two Myotis ricketti and validated our sequencing results using real-time quantitative PCR (RT-qPCR) of gene expression in the developing digits of two Hipposideros armiger. RESULTS: Among hundreds of genes exhibiting significant differences in expression between the short and long digits, we highlight 14 genes most related to digit elongation. These genes include two Tbx genes (Tbx3 and Tbx15), five BMP pathway genes (Bmp3, RGMB, Smad1, Smad4 and Nog), four Homeobox genes (Hoxd8, Hoxd9, Hoxa1 and Satb1), and three other genes (Twist1, Tmeff2 and Enpp2) related to digit malformations or cell proliferation. In addition, our results suggest that Tbx4 and Pitx2 contribute to the morphological similarity and five genes (Acta1, Tnnc2, Atp2a1, Hrc and Myoz1) contribute to the functional similarity between the thumb and hindlimb digits. CONCLUSIONS: Results of this study not only implicate many developmental genes as robust candidates underlying digit elongation and wing formation in bats, but also provide a better understanding of the genes involved in autopodial development in general.

Adaptive evolution of 5'HoxD genes in the origin and diversification of the cetacean flipper.

The homeobox (Hox) genes Hoxd12 and Hoxd13 control digit patterning and limb formation in tetrapods. Both show strong expression in the limb bud during embryonic development, are highly conserved across vertebrates, and show mutations that are associated with carpal, metacarpal, and phalangeal deformities. The most dramatic evolutionary reorganization of the mammalian limb has occurred in cetaceans (whales, dolphins, and porpoises), in which the hind limbs have been lost and the forelimbs have evolved into paddle-shaped flippers. We reconstructed the phylogeny of digit patterning in mammals and inferred that digit number has changed twice in the evolution of the cetacean forelimb. First, the divergence of the early cetaceans from their even-toed relatives coincided with the reacquisition of the pentadactyl forelimb, whereas the ancestors of tetradactyl baleen whales (Mysticeti) later lost a digit again. To test whether the evolution of the cetacean forelimb is associated with positive selection or relaxation of Hoxd12 and Hoxd13, we sequenced these genes in a wide range of mammals. In Hoxd12, we found evidence of Darwinian selection associated with both episodes of cetacean forelimb reorganization. In Hoxd13, we found a novel expansion of a polyalanine tract in cetaceans compared with other mammals (17/18 residues vs. 14/15 residues, respectively), lengthening of which has previously been shown to be linked to synpolydactyly in humans and mice. Both genes also show much greater sequence variation among cetaceans than across other mammalian lineages. Our results strongly implicate 5'HoxD genes in the modulation of digit number, web forming, and the high morphological diversity of the cetacean manus.

Global Proteomic Analysis of Lysine Malonylation in Toxoplasma gondii.

Lysine malonylation (Kmal) is a new post-translational modification (PTM), which has been reported in several prokaryotic and eukaryotic species. Although Kmal can regulate many and diverse biological processes in various organisms, knowledge about this important PTM in the apicomplexan parasite Toxoplasma gondii is limited. In this study, we performed the first global profiling of malonylated proteins in T. gondii tachyzoites using affinity enrichment and Liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis. Three experiments performed in tandem revealed 294, 345, 352 Kmal sites on 203, 236, 230 malonylated proteins, respectively. Computational analysis showed the identified malonylated proteins to be localized in various subcellular compartments and involved in many cellular functions, particularly mitochondrial function. Additionally, one conserved Kmal motif with a strong bias for cysteine was detected. Taken together, these findings provide the first report of Kmal profile in T. gondii and should be an important resource for studying the physiological roles of Kmal in this parasite.

Small intestine proteomics coupled with serum metabolomics reveal disruption of amino acid metabolism in Chinese hamsters with type 2 diabetes mellitus.

Type 2 diabetes mellitus (T2DM) is a metabolic disorder characterized by hyperglycemia, with metabolic disturbances resulting from defects in insulin secretion, insulin resistance (IR), or both. Chinese hamsters have potential value as non-obese animal models of spontaneous T2DM for studying the pathogenesis and molecular characteristics of diabetes. In this study, the molecular characteristics of the Chinese hamster diabetes animal model were investigated through small intestine proteomics and serum metabolomics. A total of 213 differentially abundant proteins and 14 differentially abundant metabolites were identified through liquid chromatography-tandem mass spectrometry (LC-MS/MS) and gas chromatography-time of flight mass spectrometry (GC-TOF/MS) analysis, respectively. Annotation by bioinformatics analysis revealed that these differentially abundant proteins in the small intestine were commonly associated with abnormal glucose and lipid metabolism, IR, impaired insulin secretion, amino acid metabolism disorders, and inflammatory dysregulation. Moreover, differentially abundant metabolites in the serum were amino acids and were related to diabetic IR. Through the analysis of small intestine proteomics and serum metabolomics in the Chinese hamster diabetes model, we provide a preliminary understanding of the diabetic characteristics of this model from a molecular perspective. This study provides data incentivizing the popularization and application of Chinese hamsters in T2DM research. SIGNIFICANCE: Spontaneous rodent models of diabetes, such as Chinese hamsters, effectively summarizes the clinical characteristics of type 2 diabetes and has high applicative value for studying the pathophysiology of diabetes. In order to explore the potential value of the Chinese hamster diabetes animal model in the study of the T2DM molecular mechanism, we performed small intestine proteomic analysis and serum metabolomic analysis in Chinese hamsters for the first time. After an integrated analysis of proteomics and metabolomics, we have a preliminary understanding of the diabetic characteristics of this model from a molecular perspective. Further, we found that in the occurrence and development of T2DM, the metabolic abnormalities of this model are particularly prominent, especially the metabolism of amino acids. These findings not only provide basic data in support of the popularization and application of the current model in T2DM research, but also provide a new perspective for the exploration of mechanisms related to T2DM.

Molecular cloning and evolutionary analysis of hemoglobin alpha-chain genes in bats.

Bats are the only mammals with the capacity for powered flight. When flying, they need abundant energy and oxygen. According to previous works, the hemoglobin (Hb) oxygen loading function of bats is insensitive to variations in body temperature, although different bat species have different heat sensitivity. We cloned Hb alpha-chain sequences from eight bat species to investigate whether they have different characteristics. We found that Hb in the bat lineages is under purifying selection, which accords with the importance of its function in bats. Three turn regions in bat Hb, however, have distinct evolutionary rates compared with those of other mammals, and the codons in these regions have an accelerated rate of evolution. These codons are under divergent selection in bats. These changes in Hb may have occurred in response to the physiological requirements of the species concerned, as adaptations to different lifestyles.

Evolutionary relationships between bat coronaviruses and their hosts.

Recent studies have suggested that bats are the natural reservoir of a range of coronaviruses (CoVs), and that rhinolophid bats harbor viruses closely related to the severe acute respiratory syndrome (SARS) CoV, which caused an outbreak of respiratory illness in humans during 2002-2003. We examined the evolutionary relationships between bat CoVs and their hosts by using sequence data of the virus RNA-dependent RNA polymerase gene and the bat cytochrome b gene. Phylogenetic analyses showed multiple incongruent associations between the phylogenies of rhinolophid bats and their CoVs, which suggested that host shifts have occurred in the recent evolutionary history of this group. These shifts may be due to either virus biologic traits or host behavioral traits. This finding has implications for the emergence of SARS and for the potential future emergence of SARS-CoVs or related viruses.

Genetic responses to seasonal variation in altitudinal stress: whole-genome resequencing of great tit in eastern Himalayas.

Species that undertake altitudinal migrations are exposed to a considerable seasonal variation in oxygen levels and temperature. How they cope with this was studied in a population of great tit (Parus major) that breeds at high elevations and winters at lower elevations in the eastern Himalayas. Comparison of population genomics of high altitudinal great tits and those living in lowlands revealed an accelerated genetic selection for carbohydrate energy metabolism (amino sugar, nucleotide sugar metabolism and insulin signaling pathways) and hypoxia response (PI3K-akt, mTOR and MAPK signaling pathways) in the high altitudinal population. The PI3K-akt, mTOR and MAPK pathways modulate the hypoxia-inducible factors, HIF-1α and VEGF protein expression thus indirectly regulate hypoxia induced angiogenesis, erythropoiesis and vasodilatation. The strategies observed in high altitudinal great tits differ from those described in a closely related species on the Tibetan Plateau, the sedentary ground tit (Parus humilis). This species has enhanced selection in lipid-specific metabolic pathways and hypoxia-inducible factor pathway (HIF-1). Comparative population genomics also revealed selection for larger body size in high altitudinal great tits.

Multiple bursts of pancreatic ribonuclease gene duplication in insect-eating bats.

Pancreatic ribonuclease gene (RNASE1) was previously shown to have undergone duplication and adaptive evolution related to digestive efficiency in several mammalian groups that have evolved foregut fermentation, including ruminants and some primates. RNASE1 gene duplications thought to be linked to diet have also been recorded in some carnivores. Of all mammals, bats have evolved the most diverse dietary specializations, mainly including frugivory and insectivory. Here we cloned, sequenced and analyzed RNASE1 gene sequences from a range of bat species to determine whether their dietary adaptation is mirrored by molecular adaptation. We found that seven insect-eating members of the families Vespertilionidae and Molossidae possessed two or more duplicates, and we also detected three pseudogenes. Reconstructed RNASE1 gene trees based on both Bayesian and maximum likelihood methods supported independent duplication events in these two families. Selection tests revealed that RNASE1 gene duplicates have undergone episodes of positive selection indicative of functional modification, and lineage-specific tests revealed strong adaptive evolution in the Tadarida ß clade. However, unlike the RNASE1 duplicates that function in digestion in some mammals, the bat RNASE1 sequences were found to be characterized by relatively high isoelectric points, a feature previously suggested to promote defense against viruses via the breakdown of double-stranded RNA. Taken together, our findings point to an adaptive diversification of RNASE1 in these two bat families, although we find no clear evidence that this was driven by diet. Future experimental assays are needed to resolve the functions of these enzymes in bats.

Ground tit genome reveals avian adaptation to living at high altitudes in the Tibetan plateau.

The ground tit (Parus humilis) is endemic to the Tibetan plateau. It is a member of family Paridae but it was long thought to be related to the ground jays because of their morphological similarities. Here we present the ground tit's genome and re-sequence two tits and one ground jay, to clarify this controversially taxonomic status and uncover its genetic adaptations to the Tibetan plateau. Our results show that ground tit groups with two tits and it diverges from them between 7.7 and 9.9 Mya. Compared with other avian genomes, ground tit shows expansion in genes linked to energy metabolism and contractions in genes involved in immune and olfactory perception. We also found positively selected and rapidly evolving genes in hypoxia response and skeletal development. These results indicated that ground tit evolves basic strategies and 'tit-to-jay' change for coping with the life in an extreme environment.