Comparative genomics reveals convergent evolution between the bamboo-eating giant and red pandas

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Impact of DNA methylation on digestive and metabolic gene expression in red pandas (Ailurus fulgens) during the transition from milk to bamboo diet

BACKGROUND

DNA methylation plays a crucial role in species development and environmental adaptation. In mammals, there are significant dietary changes from infancy to adulthood. Notably, the red panda transitions from milk consumption as juveniles to a bamboo-based diet as adults, with significant alterations in food characteristics and nutritional content. However, the regulatory role of DNA methylation in this process remains unclear. In this study, we investigate the regulatory role of DNA methylation on the expression of digestive and metabolic genes in the liver and pancreas during the red panda's dietary transition from suckling stage to adulthood.

RESULTS

Our findings reveal significant differences in DNA methylation patterns before and after dietary transition, highlighting the specific alterations in the methylation profiles of genes involved in lipid, carbohydrate, and amino acid metabolism. We found that perilipin-4 (PLIN4) is hypomethylated and highly expressed in the liver of adult red pandas, facilitating lipid droplet formation and storage, crucial for adapting to the low-fat content in bamboo. In contrast, genes like lipoprotein lipase (LPL), crucial for lipid breakdown, exhibited hypermethylated with low-expression patterns, reflecting a reduced lipid metabolism capacity in adults. Carbohydrate metabolism-related genes like ADH4 and FAM3C are hypomethylated and highly expressed in adults, enhancing glycogen production and glucose utilization. Genes involved in protein metabolism like CTSZ and GLDC, exhibit hypomethylated with high-expression and hypermethylated with low-expression patterns in the pancreas of adults, respectively, contributing to protein metabolism balance post-weaning.

CONCLUSION

This study reveals the regulatory role of DNA methylation in the dietary transition of red pandas from milk to bamboo and provides methylation evidence for the molecular regulation of adaptive expression of digestive and metabolic genes in red pandas with specialized diets.

Diversity, functions, and antibiotic resistance genes of bacteria and fungi are examined in the bamboo plant phyllosphere that serve as food for the giant pandas

The phyllosphere of bamboo is rich in microorganisms that can disrupt the intestinal microbiota of the giant pandas that consume them, potentially leading to their death. In the present study, the abundance, diversity, biological functions (e.g., KEGG and CAZyme), and antibiotic resistance genes (ARGs) of bacteria and fungi in two bamboo species phyllosphere (Chimonobambusa szechuanensis, CS; Bashania fangiana, BF) in Daxiangling Nature Reserve (an important part of the Giant Panda National Park) were investigated respectively by amplicon sequencing of the whole 16S rRNA and ITS1-ITS2 genes on PacBio Sequel and whole-metagenome shotgun sequencing on Illumina NovaSeq 6000 platform. The results suggested that there were respectively 18 bacterial and 34 fungi biomarkers between the phyllosphere of the two species of bamboo. Beta diversity of bacteria and fungi communities exited between the two bamboos according to the (un)weighted UniFrac distance matrix. Moreover, the functional analysis showed that the largest relative abundance was found in the genes related to metabolism and global and overview maps. Glycoside hydrolases (GHs) and glycosyl transferases (GTs) have a higher abundance in two bamboo phyllospheres. Co-occurrence network modeling suggested that bacteria and fungi communities in CS phyllosphere employed a much more complex metabolic network than that in BF, and the abundance of multidrug, tetracycline, and glycopeptide resistance genes was higher and closely correlated with other ARGs. This study references the basis for protecting bamboo resources foraged by wild giant pandas and predicts the risk of antibiotic resistance in bamboo phyllosphere bacterial and fungal microbiota in the Giant Panda National Park, China.

Cross-kingdom regulation of gene expression in giant pandas via plant-derived miRNA

Giant pandas (Ailuropoda melanoleuca) belong to the order Carnivora, but they mainly feed on bamboo, and their unique dietary adaptability has always been the focus of research. Recent research indicates that plant-derived microRNAs (miRNAs) can be delivered to animal organisms via exosomes and exert cross-kingdom regulatory effects on gene expression. To explore the role of plant-derived miRNAs in the dietary adaptation of giant pandas, we collected peripheral blood samples from three groups of pandas: juvenile females, adult females, and adult males-and extracted exosomes from the blood for small RNA sequencing. Additionally, three types of bamboo (shoots, stems, and leaves) consumed by the pandas were sampled for miRNA sequencing. Through comparative analysis, we identified 57 bamboo-derived miRNAs in the extracellular exosomes of giant panda peripheral blood. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses of the target genes for these miRNAs revealed their involvement in various pathways, including taste and olfactory signal transduction, digestion and absorption, and hormonal signal transduction. Furthermore, we found that plant-derived miRNAs can modulate dopamine metabolism in giant pandas, thereby influencing their food preferences. This study shows that plant-derived miRNAs can enter the bloodstream of giant pandas and exert cross-kingdom regulatory effects, potentially playing a vital role in their dietary adaptation process.

Genomic insights into marine environment adaptation and conservation of the threatened olive ridley turtle (Lepidochelys olivacea)

Sea turtles are marine flagship species and most of them are currently in a threatened state. Long-term surviving in the ocean has driven significant morphological and physiological changes for this group, which makes them an ideal model for studying adaptive evolution of marine environments. Herein, we present a chromosome-scale genome of Lepidochelys olivacea with a genome size of 2.22 Gb and a contig N50 of 97.3 Mb. Comparative genomic analyses uncovered a suite of adaptive changes in genes related to olfaction, vision, virus defense, and longevity, which may help explain the genetic underpinnings of its marine environment adaptation. We also observed that the genome-wide heterozygosity of L. olivacea was low (6.45e-4), consistent with its prolonged population decline. Overall, our study provides valuable genetic resources for understanding evolutionary adaptations to aquatic environment and for the conservation of this threatened species.

Population Genomics of Premature Termination Codons in Cavefish With Substantial Trait Loss

Loss-of-function alleles are a pertinent source of genetic variation with the potential to contribute to adaptation. Cave-adapted organisms exhibit striking loss of ancestral traits such as eyes and pigment, suggesting that loss-of-function alleles may play an outsized role in these systems. Here, we leverage 141 whole genome sequences to evaluate the evolutionary history and adaptive potential of single nucleotide premature termination codons (PTCs) in Mexican tetra. We find that cave populations contain significantly more PTCs at high frequency than surface populations. We also find that PTCs occur more frequently in genes with inherent relaxed evolutionary constraint relative to the rest of the genome. Using SLiM to simulate PTC evolution in a cavefish population, we show that the smaller population size and increased genetic drift is sufficient to account for the observed increase in PTC frequency in cave populations without positive selection. Using CRISPR-Cas9, we show that mutation of one of these genes, pde6c, produces phenotypes in surface Mexican tetra that mimic cave-derived traits. Finally, we identify a small subset of candidate genes that contain high-frequency PTCs in cave populations, occur within selective sweeps, and may contribute to beneficial traits such as reduced energy expenditure, suggesting that a handful of PTCs may be adaptive. Overall, our work provides a rare characterization of PTCs across wild populations and finds that they may have an important role in loss-of-function phenotypes, contributing to a growing body of literature showing genome evolution through relaxed constraint in subterranean organisms.

Modulating vertebrate physiology by genomic fine-tuning of GPCR functions

G protein-coupled receptors (GPCRs) play a crucial role as membrane receptors, facilitating the communication of eukaryotic species with their environment and regulating cellular and organ interactions. Consequently, GPCRs hold immense potential in contributing to adaptation to ecological niches and responding to environmental shifts. Comparative analyses of vertebrate genomes reveal patterns of GPCR gene loss, expansion, and signatures of selection. Integrating these genomic data with insights from functional analyses of gene variants enables the interpretation of genotype-phenotype correlations. This review underscores the involvement of GPCRs in adaptive processes, presenting numerous examples of how alterations in GPCR functionality influence vertebrate physiology or, conversely, how environmental changes impact GPCR functions. The findings demonstrate that modifications in GPCR function contribute to adapting to aquatic, arid, and nocturnal habitats, influencing camouflage strategies, and specializing in particular dietary preferences. Furthermore, the adaptability of GPCR functions provides an effective mechanism in facilitating past, recent, or ongoing adaptations in animal domestication and human evolution and should be considered in therapeutic strategies and drug development.

Convergent relaxation of molecular constraint in herbivores reveals the changing role of liver and kidney functions across mammalian diets

Mammalia comprises a great diversity of diet types and associated adaptations. An understanding of the genomic mechanisms underlying these adaptations may offer insights for improving human health. Comparative genomic studies of diet that employ taxonomically restricted analyses or simplified diet classifications may suffer reduced power to detect molecular convergence associated with diet evolution. Here, we use a quantitative carnivory score-indicative of the amount of animal protein in the diet-for 80 mammalian species to detect significant correlations between the relative evolutionary rates of genes and changes in diet. We have identified six genes-ACADSB, CLDN16, CPB1, PNLIP, SLC13A2, and SLC14A2-that experienced significant changes in evolutionary constraint alongside changes in carnivory score, becoming less constrained in lineages evolving more herbivorous diets. We further consider the biological functions associated with diet evolution and observe that pathways related to amino acid and lipid metabolism, biological oxidation, and small molecule transport experienced reduced purifying selection as lineages became more herbivorous. Liver and kidney functions show similar patterns of constraint with dietary change. Our results indicate that these functions are important for the consumption of animal matter and become less important with the evolution of increasing herbivory. So, genes expressed in these tissues experience a relaxation of evolutionary constraint in more herbivorous lineages.

Comparison Study on the Trophic Niche of Red Pandas Using Stable Isotope Analysis

The red panda (Ailurus), a rare and endangered mammal native to the Himalayan-Hengduan Mountains, has a specialized bamboo diet. Combining morphological and genomic evidence, red pandas have been classified as Ailurus fulgens and Ailurus styani. However, previous studies focused on ecological aspects such as foraging behaviors, habitat use and threats within specific distributions; hence, there is still a gap in quantitative comparative studies on the trophic niches of these two species. In this study, we calculated and compared the isotopic trophic niche widths of A. styani and A. fulgens by measuring carbon and nitrogen stable isotopes in hair keratin of 102 red pandas, then conducting a fecal analysis to explore the diet of a population. The results showed that (1) there was no significant difference in δ13C values between A. styani (-24.73 ± 1.26‱, N = 86) and A. fulgens (-24.43 ± 1.20‱, N = 16); however, A. styani had 1.4‱ higher average δ1⁵N value than that of A. fulgens, and the isotopic niche widths converted by Bayesian ellipse model were A. styani SEAC = 4.40‱2, A. fulgens SEAC = 4.70‱2, which might be caused by their lower and narrower altitude distribution. (2) Niche widths of three genetic populations in A. styani were XXL-LS (5.16‱2), EH-GLG (3.40‱2), QL (2.31‱2), and the comparison of the fecal composition suggested their diet diversity, which might imply a degree of competitive pressure with sympatric animals. (3) Sichuan red pandas occupied a lower trophic position compared to other herbivores but higher than giant pandas in the same ecosystem, and the niche width of red pandas was two times larger than that of sympatric giant pandas, which was related to the foraging strategies and microhabitat selection of the giant and red panda. This research provides basic stable isotopic data for red pandas, offers scientific support for the differentiated and refined conservation and management of the red panda and its habitat, and develops the application of stable isotope analysis in the ecological studies of endangered species.

Avian Island Radiations Shed Light on the Dynamics of Adaptive and Nonadaptive Radiation

Understanding the mechanisms underlying species formation and differentiation is a central goal of evolutionary biology and a formidable challenge. This understanding can provide valuable insights into the origins of the astonishing diversity of organisms living on our planet. Avian evolutionary radiations on islands have long fascinated biologists as they provide the ideal variation to study the ecological and evolutionary forces operating on the continuum between incipient lineages to complete speciation. In this review, we summarize the key insights gained from decades of research on adaptive and nonadaptive radiations of both extant and extinct insular bird species. We present a new comprehensive global list of potential avian radiations on oceanic islands, based on published island species checklists, taxonomic studies, and phylogenetic analyses. We demonstrate that our understanding of evolutionary processes is being greatly enhanced through the use of genomic tools. However, to advance the field, it is critical to complement this information with a solid understanding of the ecological and behavioral traits of both extinct and extant avian island species.

Convergent Evolution of Armor: Thermal Resistance in Deep-Sea Hydrothermal Vent Crustaceans

Organisms occupy diverse ecological niches worldwide, each with characteristics finely evolved for their environments. Crustaceans residing in deep-sea hydrothermal vents, recognized as one of Earth's extreme environments, may have adapted to withstand severe conditions, including elevated temperatures and pressure. This study compares the exoskeletons of two vent crustaceans (bythograeid crab Austinograea sp. and squat lobster Munidopsis lauensis) with four coastal species (Asian paddle crabs, blue crab, hermit crab, and mantis shrimp) to identify traits influenced by vent environments. The goal was to identify distinctive exoskeletal characteristics commonly observed in vent crustaceans, resulting from their exposure to severe abiotic factors, including elevated temperatures and pressures, found in vent environments. Results show that the exoskeletons of vent crustaceans demonstrated significantly enhanced thermal stability compared to coastal species. These vent crustaceans consistently featured exoskeletons characterized by a reduced proportion of volatile components, such as water, and an increased proportion of CaCO3, compared with coastal crustaceans. Furthermore, vent crustaceans lacked carotenoid pigments that had low heat resistance. However, no apparent differences were observed in the mechanical properties. Our findings suggest that the similar composition of exoskeletons in vent crustaceans evolved convergently to withstand high temperatures.

Distinct fecal microbiome between wild and habitat-housed captive polar bears (Ursus maritimus): Impacts of captivity and dietary shifts

Understanding the gut microbiome of polar bears can shed light on the effects of climate change-induced prolonged ice-free seasons on their health and nutritional status as a sentinel species. The fecal microbiome of habitat-housed captive polar bears who had consumed a high protein diet long-term was compared with that of the wild population. Individual differences, season, year and dietary inclusion of a brown seaweed (Fucus spiralis; part of the natural diet of wild polar bears), as a representation for nutritional change, were investigated for their effects on the fecal microbiome of captive polar bears. Microbial variations among fecal samples from wild and captive polar bears were investigated using 16s rRNA gene based metataxonomic profiling. The captive bears exhibited more diverse fecal microbiota than wild bears (p<0.05). The difference was due to significantly increased Firmicutes, Campilobacterota and Fusobacteriota, decreased Actinobacteriota (p<0.05), and absent Bdellovibrionota and Verrucomicrobiota in the captive bears. Compared with other factors, individual variation was the main driver of differences in fecal microbial composition in the captive bears. Seaweed consumption did not alter microbial diversity or composition, but this did not rule out dietary influences on the hosts. This is the first study, to the best of our knowledge, comparing the fecal microbiota of captive and wild polar bears and it reveals distinct differences between the two groups, which could result from many factors, including available food sources and the ratio of dietary macronutrients. Our findings provide preliminary insights into climate-change induced dietary shifts in polar bears related to climate-associated habitat change.

The Chinese longsnout catfish genome provides novel insights into the feeding preference and corresponding metabolic strategy of carnivores

Fish show variation in feeding habits to adapt to complex environments. However, the genetic basis of feeding preference and the corresponding metabolic strategies that differentiate feeding habits remain elusive. Here, by comparing the whole genome of a typical carnivorous fish (Leiocassis longirostris Günther) with that of herbivorous fish, we identify 250 genes through both positive selection and rapid evolution, including taste receptor taste receptor type 1 member 3 (tas1r3) and trypsin We demonstrate that tas1r3 is required for carnivore preference in tas1r3-deficient zebrafish and in a diet-shifted grass carp model. We confirm that trypsin correlates with the metabolic strategies of fish with distinct feeding habits. Furthermore, marked alterations in trypsin activity and metabolic profiles are accompanied by a transition of feeding preference in tas1r3-deficient zebrafish and diet-shifted grass carp. Our results reveal a conserved adaptation between feeding preference and corresponding metabolic strategies in fish, and provide novel insights into the adaptation of feeding habits over the evolution course.

Comparative genomics sheds new light on the convergent evolution of infrared vision in snakes

Infrared vision is a highly specialized sensory system that evolved independently in three clades of snakes. Apparently, convergent evolution occurred in the transient receptor potential ankyrin 1 (TRPA1) proteins of infrared-sensing snakes. However, this gene can only explain how infrared signals are received, and not the transduction and processing of those signals. We sequenced the genome of Xenopeltis unicolor, a key outgroup species of pythons, and performed a genome-wide analysis of convergence between two clades of infrared-sensing snakes. Our results revealed pervasive molecular adaptation in pathways associated with neural development and other functions, with parallel selection on loci associated with trigeminal nerve structural organization. In addition, we found evidence of convergent amino acid substitutions in a set of genes, including TRPA1 and TRPM2. The analysis also identified convergent accelerated evolution in non-coding elements near 12 genes involved in facial nerve structural organization and optic nerve development. Thus, convergent evolution occurred across multiple dimensions of infrared vision in vipers and pythons, as well as amino acid substitutions, non-coding elements, genes and functions. These changes enabled independent groups of snakes to develop and use infrared vision.

Description and molecular characterisation of Babesia ailuropodae n. sp., a new piroplasmid species infecting giant pandas

BACKGROUND

Babesia spp. are protozoan parasites that infect the red blood cells of domesticated animals, wildlife and humans. A few cases of giant pandas (a flagship species in terms of wildlife conservation) infected with a putative novel Babesia sp. have been reported. However, comprehensive research on the morphological and molecular taxonomic classification of this novel Babesia sp. is still lacking. This study was designed to close this gap and formally describe this new Babesia sp. infecting giant pandas.

METHODS

Detailed morphological, molecular and phylogenetic analyses were conducted to characterise this Babesia sp. and to assess its systematic relationships with other Babesia spp. Blood samples from giant pandas infected with Babesia were subjected to microscopic examination. The 18S ribosomal RNA (18S rRNA), cytochrome b (cytb) and mitochondrial genome (mitogenome) of the new Babesia sp. were amplified, sequenced and assembled using DNA purified from blood samples taken from infected giant pandas. Based on the newly generated 18S rRNA, cytb and mitogenome sequences, phylogenetic trees were constructed.

RESULTS

Morphologically, the Babesia sp. from giant pandas exhibited various forms, including round to oval ring-shaped morphologies, resembling those found in other small canine Babesia spp. and displaying typical tetrads. Phylogenetic analyses with the 18S rRNA, cytb and mitogenome sequences revealed that the new Babesia sp. forms a monophyletic group, with a close phylogenetic relationship with the Babesia spp. that infect bears (Ursidae), raccoons (Procyonidae) and canids (Canidae). Notably, the mitogenome structure consisted of six ribosomal large subunit-coding genes (LSU1-6) and three protein-coding genes (cytb, cox3 and cox1) arranged linearly.

CONCLUSIONS

Based on coupled morphological and genetic analyses, we describe a novel species of the genus Babesia, namely, Babesia ailuropodae n. sp., which infects giant pandas.

Adaptive Expression and ncRNA Regulation of Genes Related to Digestion and Metabolism in Stomach of Red Pandas during Suckling and Adult Periods

Red pandas evolved from carnivores to herbivores and are unique within Carnivora. Red pandas and carnivorous mammals consume milk during the suckling period, while they consume bamboo and meat during the adult period, respectively. Red pandas and carnivorous mammal ferrets have a close phylogenetic relationship. To further investigate the molecular mechanisms of dietary changes and nutrient utilization in red pandas from suckling to adult, comparative analysis of the whole transcriptome was performed on stomach tissues from red pandas and ferrets during the suckling and adult periods. The main results are as follows: (1) we identified ncRNAs for the first time in stomach tissues of both species, and found significant expression changes of 109 lncRNAs and 106 miRNAs in red pandas and 756 lncRNAs and 109 miRNAs in ferrets between the two periods; (2) up-regulated genes related to amino acid transport regulated by lncRNA-miRNA-mRNA networks may efficiently utilize limited bamboo amino acids in adult red pandas, while up-regulated genes related to amino acid degradation regulated by lncRNAs may maintain the balance of amino acid metabolism due to larger daily intakes in adult ferrets; and (3) some up-regulated genes related to lipid digestion may contribute to the utilization of rich nutrients in milk for the rapid growth and development of suckling red pandas, while up-regulated genes associated with linoleic acid metabolism regulated by lncRNA-miRNA-mRNA networks may promote cholesterol decomposition to reduce health risks for carnivorous adult ferrets. Collectively, our study offers evidence of gene expression adaptation and ncRNA regulation in response to specific dietary changes and nutrient utilization in red pandas during suckling and adult periods.

Structural variant landscapes reveal convergent signatures of evolution in sheep and goats

BACKGROUND

Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution.

RESULTS

We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats.

CONCLUSIONS

Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.

Comparative transcriptome and methylome of polar bears, giant and red pandas reveal diet-driven adaptive evolution

Epigenetic regulation plays an important role in the evolution of species adaptations, yet little information is available on the epigenetic mechanisms underlying the adaptive evolution of bamboo-eating in both giant pandas (Ailuropoda melanoleuca) and red pandas (Ailurus fulgens). To investigate the potential contribution of epigenetic to the adaptive evolution of bamboo-eating in giant and red pandas, we performed hepatic comparative transcriptome and methylome analyses between bamboo-eating pandas and carnivorous polar bears (Ursus maritimus). We found that genes involved in carbohydrate, lipid, amino acid, and protein metabolism showed significant differences in methylation and expression levels between the two panda species and polar bears. Clustering analysis of gene expression revealed that giant pandas did not form a sister group with the more closely related polar bears, suggesting that the expression pattern of genes in livers of giant pandas and red pandas have evolved convergently driven by their similar diets. Compared to polar bears, some key genes involved in carbohydrate metabolism and biological oxidation and cholesterol synthesis showed hypomethylation and higher expression in giant and red pandas, while genes involved in fat digestion and absorption, fatty acid metabolism, lysine degradation, resistance to lipid peroxidation and detoxification showed hypermethylation and low expression. Our study elucidates the special nutrient utilization mechanism of giant pandas and red pandas and provides some insights into the molecular mechanism of their adaptive evolution of bamboo feeding. This has important implications for the breeding and conservation of giant pandas and red pandas.

WGCCRR: a web-based tool for genome-wide screening of convergent indels and substitutions of amino acids

Summary

Genome-wide analyses of proteincoding gene sequences are being employed to examine the genetic basis of adaptive evolution in many organismal groups. Previous studies have revealed that convergent/parallel adaptive evolution may be caused by convergent/parallel amino acid changes. Similarly, detailed analysis of lineage-specific amino acid changes has shown correlations with certain lineage-specific traits. However, experimental validation remains the ultimate measure of causality. With the increasing availability of genomic data, a streamlined tool for such analyses would facilitate and expedite the screening of genetic loci that hold potential for adaptive evolution, while alleviating the bioinformatic burden for experimental biologists. In this study, we present a user-friendly web-based tool called WGCCRR (Whole Genome Comparative Coding Region Read) designed to screen both convergent/parallel and lineage-specific amino acid changes on a genome-wide scale. Our tool allows users to replicate previous analyses with just a few clicks, and the exported results are straightforward to interpret. In addition, we have also included amino acid indels that are usually neglected in previous work. Our website provides an efficient platform for screening candidate loci for downstream experimental tests.

Availability and Implementation

The tool is available at: https://fishevo.xmu.edu.cn/.

Accurate Detection of Convergent Mutations in Large Protein Alignments With ConDor

Evolutionary convergences are observed at all levels, from phenotype to DNA and protein sequences, and changes at these different levels tend to be correlated. Notably, convergent mutations can lead to convergent changes in phenotype, such as changes in metabolism, drug resistance, and other adaptations to changing environments. We propose a two-component approach to detect mutations subject to convergent evolution in protein alignments. The "Emergence" component selects mutations that emerge more often than expected, while the "Correlation" component selects mutations that correlate with the convergent phenotype under study. With regard to Emergence, a phylogeny deduced from the alignment is provided by the user and is used to simulate the evolution of each alignment position. These simulations allow us to estimate the expected number of mutations in a neutral model, which is compared to the observed number of mutations in the data studied. In Correlation, a comparative phylogenetic approach, is used to measure whether the presence of each of the observed mutations is correlated with the convergent phenotype. Each component can be used on its own, for example Emergence when no phenotype is available. Our method is implemented in a standalone workflow and a webserver, called ConDor. We evaluate the properties of ConDor using simulated data, and we apply it to three real datasets: sedge PEPC proteins, HIV reverse transcriptase, and fish rhodopsin. The results show that the two components of ConDor complement each other, with an overall accuracy that compares favorably to other available tools, especially on large datasets.

Ursids evolved dietary diversity without major alterations in metabolic rates

The diets of the eight species of ursids range from carnivory (e.g., polar bears, Ursus maritimus) to insectivory (e.g., sloth bears, Melursus ursinus), omnivory (e.g., brown bears, U. arctos), and herbivory (e.g., giant pandas, Ailuropoda melanoleuca). Dietary energy availability ranges from the high-fat, highly digestible, calorically dense diet of polar bears (~ 6.4 kcal digestible energy/g fresh weight) to the high-fiber, poorly digestible, calorically restricted diet (~ 0.7) of giant pandas. Thus, ursids provide the opportunity to examine the extent to which dietary energy drives evolution of energy metabolism in a closely related group of animals. We measured the daily energy expenditure (DEE) of captive brown bears in a relatively large, zoo-type enclosure and compared those values to previously published results on captive brown bears, captive and free-ranging polar bears, and captive and free-ranging giant pandas. We found that all three species have similar mass-specific DEE when travel distances and energy intake are normalized even though their diets differ dramatically and phylogenetic lineages are separated by millions of years. For giant pandas, the ability to engage in low-cost stationary foraging relative to more wide-ranging bears likely provided the necessary energy savings to become bamboo specialists without greatly altering their metabolic rate.

Conservation Genomics and Metagenomics of Giant and Red Pandas in the Wild

Giant pandas and red pandas are endangered species with similar specialized bamboo diet and partial sympatric distribution in China. Over the last two decades, the rapid development of genomics and metagenomics research on these species has enriched our knowledge of their biology, ecology, physiology, genetics, and evolution, which is crucial and useful for their conservation. We describe the evolutionary history, endangerment processes, genetic diversity, and population structure of wild giant pandas and two species of red pandas (Chinese and Himalayan red pandas). In addition, we explore how genomics and metagenomics studies have provided insight into the convergent adaptation of pandas to the specialized bamboo diet. Finally, we discuss how these findings are applied to effective conservation management of giant and red pandas in the wild and in captivity to promote the long-term persistence of these species.

A chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution

The electric catfish (Malapterurus electricus), belonging to the family Malapteruridae, order Siluriformes (Actinopterygii: Ostariophysi), is one of the six branches that has independently evolved electrical organs. We assembled a 796.75 Mb M. electricus genome and anchored 88.72% sequences into 28 chromosomes. Gene family analysis revealed 295 expanded gene families that were enriched on functions related to glutamate receptors. Convergent evolutionary analyses of electric organs among different lineage of electric fishes further revealed that the coding gene of rho guanine nucleotide exchange factor 4-like (arhgef4), which is associated with G-protein coupled receptor (GPCR) signaling pathway, underwent adaptive parallel evolution. Gene identification suggests visual degradation in catfishes, and an important role for taste in environmental adaptation. Our findings fill in the genomic data for a branch of electric fish and provide a relevant genetic basis for the adaptive evolution of Siluriformes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00197-8.

The study on the impact of sex on the structure of gut microbiota of bamboo rats in China

Introduction

Bamboo rats are rodents that eat bamboo, and their robust capacity for bamboo digestion is directly correlated with their gut flora. Chinese bamboo rat (Rhizomys sinensis) is a common bamboo rat in Chinese central and southern regions. As a single-stomach mammal, bamboo rats are a famous specificity bamboo-eating animal and their intestinal microbial composition may also play a key role in the digestion of cellulose and lignin. So, the gut microbiota of bamboo rat may play an important role in the adaptation of bamboo rats for digesting lignocellulose-based diet.

Methods

To study the microbiome differences of bamboo rats from different sexes, the microbial genomic DNA was extracted from each fecal sample and the V4 region of 16S rRNA genes was amplified and sequencing on an IlluminaHiSeq6000 platform. The operational taxonomic units (OTUs) were classified, the OTUs in different sexes was identified and compared at phylum and genus levels. For isolation and screening of cellulose degradation bacteria from bamboo rats, fresh feces from randomly selected bamboo rats were collected and used for the isolation and screening of cellulose degradation bacteria using Luria Bertani (LB) Agar medium containing Carboxymethyl cellulose. The cellulase activity, biochemical characterization and phylogenetic analysis of the purified bacteria strains were characterized.

Results and discussion

A total of 3,833 OTUs were classified. The total microbial diversity detected in the female and male rats was 3,049 OTUs and 3,452 OTUs, respectively. The Shannon index revealed significant differences between the two groups (p < 0.05), though they were all captive and had the same feeding conditions. At the phylum level, Firmicutes, Bacteroidota, and Proteobacteria were prominent in the microbial community. At the genus level, the microbial community was dominated by Lachnospiraceae, Lactobacillus, Bacteroides, and Prevotella, but there was a significant difference between the two groups of bamboo rats; ~90 bacteria genus in the female group was significantly higher than the male group. Among them, Bacteroides, Colidextribacter, and Oscillibacter were significantly higher genera, and the genera of Lachnoclostridium, Oscillibacter, and Papillibacter had the highest FC value among the male and female bamboo rats. The KEGG function annotation and different pathways analysis revealed that membrane transport, carbohydrate metabolism, and amino acid metabolism were the most enriched metabolic pathways in the two groups, and multiple sugar transport system permease protein (K02025 and K02026), RNA polymerase sigma-70 factor (K03088), and ATP-binding cassette (K06147) were the three different KEGG pathways (p < 0.05). Two cellulose degradation bacteria strains-Bacillus subtilis and Enterococcus faecalis-were isolated and characterized from the feces of bamboo rats.

Genomic signatures of convergent shifts to plunge-diving behavior in birds

Understanding the genetic basis of convergence at broad phylogenetic scales remains a key challenge in biology. Kingfishers (Aves: Alcedinidae) are a cosmopolitan avian radiation with diverse colors, diets, and feeding behaviors-including the archetypal plunge-dive into water. Given the sensory and locomotor challenges associated with air-water transitions, kingfishers offer a powerful opportunity to explore the effects of convergent behaviors on the evolution of genomes and phenotypes, as well as direct comparisons between continental and island lineages. Here, we use whole-genome sequencing of 30 diverse kingfisher species to identify the genomic signatures associated with convergent feeding behaviors. We show that species with smaller ranges (i.e., on islands) have experienced stronger demographic fluctuations than those on continents, and that these differences have influenced the dynamics of molecular evolution. Comparative genomic analyses reveal positive selection and genomic convergence in brain and dietary genes in plunge-divers. These findings enhance our understanding of the connections between genotype and phenotype in a diverse avian radiation.

The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

Comparative study of the digestion and metabolism related genes' expression changes during the postnatal food change in different dietary mammals

The changes in the expression of genes related to digestion and metabolism may be various in different dietary mammals from juvenile to adult, especially, the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens), which were once carnivores but have shifted to being specialized bamboo eaters, are unique features of their changes are more unclear. To elucidate the changing patterns of gene expression related to digestion and metabolism from juvenile to adult in different dietary mammals, we performed transcriptome analysis of the liver or pancreas in giant and red pandas, herbivorous rabbits (Oryctolagus cuniculus) and macaques (Macaca mulatta), carnivorous ferrets (Mustela putorius furo), and omnivorous mice (Mus musculus) from juvenile to adult. During the transition from juvenile to adulthood, giant and red pandas, as well as rabbits and macaques, show significant upregulation of key genes for carbohydrate metabolism, such as starch hydrolysis and sucrose metabolism, and unsaturated fatty acid metabolism, such as linoleic acid, while there is no significant difference in the expression of key genes for fatty acid β-oxidation. A large number of amino acid metabolism related genes were upregulated in adult rabbits and macaques compared to juveniles. While adult giant and red pandas mainly showed upregulation of key genes for arginine synthesis and downregulation of key genes for arginine and lysine degradation. In adult stages, mouse had significantly higher expression patterns in key genes for starch hydrolysis and sucrose metabolism, as well as lipid and protein metabolism. In contrast to general expectations, genes related to lipid, amino acid and protein metabolism were significantly higher expressed in adult group of ferrets, which may be related to their high metabolic levels. Our study elucidates the pattern of changes in the expression of genes related to digestion and metabolism from juvenile to adult in different dietary mammals, with giant and red pandas showing adaptations associated with specific nutritional limitations of bamboo.

Convergent molecular evolution of thermogenesis and circadian rhythm in Arctic ruminants

The muskox and reindeer are the only ruminants that have evolved to survive in harsh Arctic environments. However, the genetic basis of this Arctic adaptation remains largely unclear. Here, we compared a de novo assembled muskox genome with reindeer and other ruminant genomes to identify convergent amino acid substitutions, rapidly evolving genes and positively selected genes among the two Arctic ruminants. We found these candidate genes were mainly involved in brown adipose tissue (BAT) thermogenesis and circadian rhythm. Furthermore, by integrating transcriptomic data from goat adipose tissues (white and brown), we demonstrated that muskox and reindeer may have evolved modulating mitochondrion, lipid metabolism and angiogenesis pathways to enhance BAT thermogenesis. In addition, results from co-immunoprecipitation experiments prove that convergent amino acid substitution of the angiogenesis-related gene hypoxia-inducible factor 2alpha (HIF2A), resulting in weakening of its interaction with prolyl hydroxylase domain-containing protein 2 (PHD2), may increase angiogenesis of BAT. Altogether, our work provides new insights into the molecular mechanisms involved in Arctic adaptation.

Selection-driven trait loss in independently evolved cavefish populations

Laboratory studies have demonstrated that a single phenotype can be produced by many different genotypes; however, in natural systems, it is frequently found that phenotypic convergence is due to parallel genetic changes. This suggests a substantial role for constraint and determinism in evolution and indicates that certain mutations are more likely to contribute to phenotypic evolution. Here we use whole genome resequencing in the Mexican tetra, Astyanax mexicanus, to investigate how selection has shaped the repeated evolution of both trait loss and enhancement across independent cavefish lineages. We show that selection on standing genetic variation and de novo mutations both contribute substantially to repeated adaptation. Our findings provide empirical support for the hypothesis that genes with larger mutational targets are more likely to be the substrate of repeated evolution and indicate that features of the cave environment may impact the rate at which mutations occur.

Natural selection and convergent evolution of the HOX gene family in Carnivora

HOX genes play a central role in the development and regulation of limb patterns. For mammals in the order Carnivora, limbs have evolved in different forms, and there are interesting cases of phenotypic convergence, such as the pseudothumb of the giant and red pandas, and the flippers or specialized limbs of the pinnipeds and sea otter. However, the molecular bases of limb development remain largely unclear. Here, we studied the molecular evolution of the HOX9 ~ 13 genes of 14 representative species in Carnivora and explored the molecular evolution of other HOX genes. We found that only one limb development gene, HOXC10, underwent convergent evolution between giant and red pandas and was thus an important candidate gene related to the development of pseudothumbs. No signals of amino acid convergence and natural selection were found in HOX9 ~ 13 genes between pinnipeds and sea otter, but there was evidence of positive selection and rapid evolution in four pinniped species. Overall, few HOX genes evolve via natural selection or convergent evolution, and these could be important candidate genes for further functional validation. Our findings provide insights into potential molecular mechanisms of the development of specialized pseudothumbs and flippers (or specialized limbs).

Comparative Analysis of Olfactory Receptor Repertoires Sheds Light on the Diet Adaptation of the Bamboo-Eating Giant Panda Based on the Chromosome-Level Genome

The giant panda (Ailuropoda melanoleuca) is the epitome of a flagship species for wildlife conservation and also an ideal model of adaptive evolution. As an obligate bamboo feeder, the giant panda relies on the olfaction for food recognition. The number of olfactory receptor (OR) genes and the rate of pseudogenes are the main factors affecting the olfactory ability of animals. In this study, we used the chromosome-level genome of the giant panda to identify OR genes and compared the genome sequences of OR genes with five other Ursidae species (spectacled bear (Tremarctos ornatus), American black bear (Ursus americanus), brown bear (Ursus arctos), polar bear (Ursus maritimus) and Asian black bear (Ursus thibetanus)). The giant panda had 639 OR genes, including 408 functional genes, 94 partial OR genes and 137 pseudogenes. Among them, 222 OR genes were detected and distributed on 18 chromosomes, and chromosome 8 had the most OR genes. A total of 448, 617, 582, 521 and 792 OR genes were identified in the spectacled bear, American black bear, brown bear, polar bear and Asian black bear, respectively. Clustering analysis based on the OR protein sequences of the six species showed that the OR genes distributed in 69 families and 438 subfamilies based on sequence similarity, and the six mammals shared 72 OR gene subfamilies, while the giant panda had 31 unique OR gene subfamilies (containing 35 genes). Among the 35 genes, there are 10 genes clustered into 8 clusters with 10 known human OR genes (OR8J3, OR51I1, OR10AC1, OR1S2, OR1S1, OR51S1, OR4M1, OR4M2, OR51T1 and OR5W2). However, the kind of odor molecules can be recognized by the 10 known human OR genes separately, which needs further research. The phylogenetic tree showed that 345 (about 84.56%) functional OR genes were clustered as Class-II, while only 63 (about 15.44%) functional OR genes were clustered as Class-I, which required further and more in-depth research. The potential odor specificity of some giant panda OR genes was identified through the similarity to human protein sequences. Sequences similar to OR2B1, OR10G3, OR11H6 and OR11H7P were giant panda-specific lacking, which may be related to the transformation and specialization from carnivore to herbivore of the giant panda. Since our reference to flavoring agents comes from human research, the possible flavoring agents from giant panda-specific OR genes need further investigation. Moreover, the conserved motifs of OR genes were highly conserved in Ursidae species. This systematic study of OR genes in the giant panda will provide a solid foundation for further research on the olfactory function and variation of the giant panda.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Chromosome-Length Assembly of the Baikal Seal (Pusa sibirica) Genome Reveals a Historically Large Population Prior to Isolation in Lake Baikal

Pusa sibirica, the Baikal seal, is the only extant, exclusively freshwater, pinniped species. The pending issue is, how and when they reached their current habitat—the rift lake Baikal, more than three thousand kilometers away from the Arctic Ocean. To explore the demographic history and genetic diversity of this species, we generated a de novo chromosome-length assembly, and compared it with three closely related marine pinniped species. Multiple whole genome alignment of the four species compared with their karyotypes showed high conservation of chromosomal features, except for three large inversions on chromosome VI. We found the mean heterozygosity of the studied Baikal seal individuals was relatively low (0.61 SNPs/kbp), but comparable to other analyzed pinniped samples. Demographic reconstruction of seals revealed differing trajectories, yet remarkable variations in Ne occurred during approximately the same time periods. The Baikal seal showed a significantly more severe decline relative to other species. This could be due to the difference in environmental conditions encountered by the earlier populations of Baikal seals, as ice sheets changed during glacial–interglacial cycles. We connect this period to the time of migration to Lake Baikal, which occurred ~3–0.3 Mya, after which the population stabilized, indicating balanced habitat conditions.

Gene expressions between obligate bamboo-eating pandas and non-herbivorous mammals reveal converged specialized bamboo diet adaptation

BACKGROUND

It is inevitable to change the function or expression of genes during the environmental adaption of species. Both the giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens) belong to Carnivora and have developed similar adaptations to the same dietary switch to bamboos at the morphological and genomic levels. However, the genetic adaptation at the gene expression level is unclear. Therefore, we aimed to examine the gene expression patterns of giant and red panda convergent specialized bamboo-diets. We examined differences in liver and pancreas transcriptomes between the two panda species and other non-herbivorous species.

RESULTS

The clustering and PCA plots suggested that the specialized bamboo diet may drive similar expression shifts in these two species of pandas. Therefore, we focused on shared liver and pancreas DEGs (differentially expressed genes) in the giant and red panda relative to other non-herbivorous species. Genetic convergence occurred at multiple levels spanning carbohydrate metabolism, lipid metabolism, and lysine degradation. The shared adaptive convergence DEGs in both organs probably be an evolutionary response to the high carbohydrate, low lipid and lysine bamboo diet. Convergent expression of those nutrient metabolism-related genes in both pandas was an intricate process and subjected to multi-level regulation, including DNA methylation and transcription factor. A large number of lysine degradation and lipid metabolism related genes were hypermethylated in promoter regions in the red panda. Most genes related to carbohydrate metabolism had reduced DNA methylation with increased mRNA expression in giant pandas. Unlike the red panda, the core gene of the lysine degradation pathway (AASS) doesn't exhibit hypermethylation modification in the giant panda, and dual-luciferase reporter assay showed that transcription factor, NR3C1, functions as a transcriptional activator in AASS transcription through the binding to AASS promoter region.

CONCLUSIONS

Our results revealed the adaptive expressions and regulations of the metabolism-related genes responding to the unique nutrients in bamboo food and provided data accumulation and research hints for the future revelation of complex mechanism of two pandas underlying convergent adaptation to a specialized bamboo diet.

Similar adaptative mechanism but divergent demographic history of four sympatric desert rodents in Eurasian inland

Phenotypes associated with metabolism and water retention are thought to be key to the adaptation of desert species. However, knowledge on the genetic changes and selective regimes on the similar and divergent ways to desert adaptation in sympatric and phylogenetically close desert organisms remains limited. Here, we generate a chromosome level genome assembly for Northern three-toed jerboa (Dipus sagitta) and three other high-quality genome assemblies for Siberian jerboa (Orientallactaga sibirica), Midday jird (Meriones meridianus), and Desert hamster (Phodopus roborovskii). Genomic analyses unveil that desert adaptation of the four species mainly result from similar metabolic pathways, such as arachidonic acid metabolism, thermogenesis, oxidative phosphorylation, insulin related pathway, DNA repair and protein synthesis and degradation. However, the specific evolved genes in the same adaptative molecular pathway often differ in the four species. We also reveal similar niche selection but different demographic histories and sensitivity to climate changes, which may be related to the diversified genomic adaptative features. In addition, our study suggests that nocturnal rodents have evolved some specific adaptative mechanism to desert environments compared to large desert animals. Our genomic resources will provide an important foundation for further research on desert genetic adaptations.

GWAS Atlas: an updated knowledgebase integrating more curated associations in plants and animals

GWAS Atlas (https://ngdc.cncb.ac.cn/gwas/) is a manually curated resource of genome-wide genotype-to-phenotype associations for a wide range of species. Here, we present an updated implementation of GWAS Atlas by curating and incorporating more high-quality associations, with significant improvements and advances over the previous version. Specifically, the current release of GWAS Atlas incorporates a total of 278,109 curated genotype-to-phenotype associations for 1,444 different traits across 15 species (10 plants and 5 animals) from 830 publications and 3,432 studies. A collection of 6,084 lead SNPs of 439 traits and 486 experiment-validated causal variants of 157 traits are newly added. Moreover, 1,056 trait ontology terms are newly defined, resulting in 1,172 and 431 terms for Plant Phenotype and Trait Ontology and Animal Phenotype and Trait Ontology, respectively. Additionally, it is equipped with four online analysis tools and a submission platform, allowing users to perform data analysis and data submission. Collectively, as a core resource in the National Genomics Data Center, GWAS Atlas provides valuable genotype-to-phenotype associations for a diversity of species and thus plays an important role in agronomic trait study and molecular breeding.

Convergent Genomic Signatures of High-Altitude Adaptation among Six Independently Evolved Mammals

The species living in the Qinghai-Tibet Plateau provide an excellent model system for studying the relationship between molecular convergent evolution and adaptation. Distant species experiencing the same selection pressure (i.e., hypoxia, low temperature and strong ultraviolet radiation) are likely to evolve similar genetic adaptations independently. Here, we performed comparative genomics studies on six independently evolved high-altitude species. The results also showed that the convergent evolution of the six species was mainly reflected at the level of rapidly evolving genes, and the functions of these rapidly evolving genes were mainly related to hypoxia response and DNA damage repair. In addition, we found that high-altitude species had more gene family changes than their low-altitude relatives, except for the order Lagomorpha. The results also show that the convergence of the gene family contraction of high-altitude species is much greater than that of expansion, revealing a possible pattern of species in adapting to high-altitude. Furthermore, we detected a positive selection signature in four genes related to hypoxia response and ultraviolet radiation damage in these six species (FYCO1, ERBIN, SCAMP1 and CXCL10). Our study reveals that hypoxia response might play an important role in the adaptation of independently evolved species to a high-altitude environment, providing a basic perspective for further exploring the high-altitude adaptation mechanism of different related species in the future.

Diet evolution of carnivorous and herbivorous mammals in Laurasiatheria

Background

Laurasiatheria contains taxa with diverse diets, while the molecular basis and evolutionary history underlying their dietary diversification are less clear.

Results

In this study, we used the recently developed molecular phyloecological approach to examine the adaptive evolution of digestive system-related genes across both carnivorous and herbivorous mammals within Laurasiatheria. Our results show an intensified selection of fat and/or protein utilization across all examined carnivorous lineages, which is consistent with their high-protein and high-fat diets. Intriguingly, for herbivorous lineages (ungulates), which have a high-carbohydrate diet, they show a similar selection pattern as that of carnivorous lineages. Our results suggest that for the ungulates, which have a specialized digestive system, the selection intensity of their digestive system-related genes does not necessarily reflect loads of the nutrient components in their diets but appears to be positively related to the loads of the nutrient components that are capable of being directly utilized by the herbivores themselves. Based on these findings, we reconstructed the dietary evolution within Laurasiatheria, and our results reveal the dominant carnivory during the early diversification of Laurasiatheria. In particular, our results suggest that the ancestral bats and the common ancestor of ruminants and cetaceans may be carnivorous as well. We also found evidence of the convergent evolution of one fat utilization-related gene, APOB, across carnivorous taxa.

Conclusions

Our molecular phyloecological results suggest that digestive system-related genes can be used to determine the molecular basis of diet differentiations and to reconstruct ancestral diets.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02033-6.

The metabolic adaptation in wild vertebrates via omics approaches

Metabolism is the basis for sustaining life and essential to the adaptive evolution of organisms. With the development of high-throughput sequencing technology, genetic mechanisms of adaptive evolution, including metabolic adaptation, have been extensively resolved by omics approaches, but a deep understanding of genetic and epigenetic metabolic adaptation is still lacking. Exploring metabolic adaptations from genetic and epigenetic perspectives in wild vertebrates is vital to understanding species evolution, especially for the early stages of adaptative evolution. Herein, we summarize the advances in our understanding of metabolic adaptations via omics approaches in wild vertebrates based on three types of cases: extreme environment, periodically changing environment, and changes of species characteristics. We conclude that the understanding of the formation of metabolic adaptations at the genetic level alone can well identify the adaptive genetic variation that has developed during evolution, but cannot resolve the potential impact of metabolic adaptations on the adaptative evolution in the future. Thus, it seems imperative to include epigenomics and metabolomics in the study of adaptation, and that in the future genomic and epigenetic data should be integrated to understand the formation of metabolic adaptation of wild vertebrate organisms.

A comparative study of skin transcriptomes and histological observations for black and white hair colors of giant panda

The Giant pandas (Ailuropoda melanoleuca) are mammals belonging to the bear family, order Carnivora, and their characteristic hair color and distribution has been in the spotlight. In recent years, the gradual prevalence of skin diseases in giant pandas and even the discovery of albino individuals have made the study of the substrate of their skin hair distribution more and more urgent. In this study, by comparing the skin histology and transcriptomes for hairs of different color of giant pandas, we found that the melanin contents of hair follicles at the bases of black and white hairs differed, but the hair follicles at the base of white hairs also contained some amount of melanin. The transcriptome sequencing results showed that there were great differences in the expression of the transcriptome of the skin under different hair color blocks, in which the number of differentially expressed genes in the white skin was much smaller than that in the black skin. Transcriptomes for skin tissue samples for different hair colors revealed several enriched Kyoto encyclopedia of genes (KEGG) pathways that include tumor, cell adhesion and melanocyte growth-related signaling pathways. This study provides a theoretical basis for subsequent studies on hair color distribution and skin diseases in giant pandas.

Role of feeding specialization in taste receptor loss: insights from sweet and umami receptor evolution in Carnivora

Controversy and misunderstanding surround the role of feeding specialization in taste receptor loss in vertebrates. We refined and tested the hypothesis that this loss is caused by feeding specializations. Specifically, feeding specializations were proposed to trigger time-dependent process of taste receptor loss through deprivation of benefit of using the receptor's gustatory function. We propose that this process may be accelerated by abiotic environmental conditions or decelerated/stopped because of extragustatory functions of the receptor's protein(s). As test case we used evolution of the sweet (TAS1R2+TAS1R3) and umami (TAS1R1+TAS1R3) receptors in Carnivora (dogs, cats, and kin). We predicted these receptors' absence/presence using data on presence/absence of inactivating mutations in these receptors' genes and data from behavioral sweet/umami preference tests. We identified 20 evolutionary events of sweet (11) or umami (9) receptor loss. These events affected species with feeding specializations predicted to favor sweet/umami receptor loss (27 and 22 species, respectively). All species with feeding habits predicted to favor sweet/umami receptor retention (11 and 24, respectively) were found to retain that receptor. Six species retained the sweet (5) or umami (1) receptor despite feeding specialization predicted to favor loss of that receptor, which can be explained by the time dependence of sweet/umami receptor loss process and the possible decelerating effect of TAS1R extragustatory functions so that the sweet/umami receptor process is ongoing in these species. Our findings support the idea that feeding specialization leads to taste receptor loss and is the main if not only triggering factor for evolutionary loss of taste receptors.

An improved, chromosome-level genome of the giant panda (Ailuropoda melanoleuca)

BACKGROUND

The iconic giant panda (Ailuropoda melanoleuca), as both a flagship and umbrella species endemic to China, is a world famous symbol for wildlife conservation. The giant panda has several specific biological traits and holds a relatively small place in evolution. A high-quality genome of the giant panda is key to understanding the biology of this vulnerable species.

FINDINGS

We generated a 2.48-Gb chromosome-level genome (GPv1) of the giant panda named "Jing Jing" with a contig N50 of 28.56 Mb and scaffold N50 of 134.17 Mb, respectively. The total length of chromosomes (n = 21) was 2.39-Gb, accounting for 96.4% of the whole genome. Compared with the previously published four genomes of the giant panda, our genome is characterized by the highest completeness and the correct sequence orientation. A gap-free and 850 kb length of immunoglobulin heavy-chain gene cluster was manually annotated in close proximity to the telomere of chromosome 14. Additionally, we developed an algorithm to predict the centromere position of each chromosome. We also constructed a complete chromatin structure for "Jing Jing", which includes inter-chromosome interaction pattern, A/B compartment, topologically associated domain (TAD), TAD-clique and promoter-enhancer interaction (PEI).

CONCLUSIONS

We presented an improved chromosome-level genome and complete chromatin structure for the giant panda. This is a valuable resource for the future genetic and genomic studies on giant panda.

Animal-SNPAtlas: a comprehensive SNP database for multiple animals

Single-nucleotide polymorphisms (SNPs) as the most important type of genetic variation are widely used in describing population characteristics and play vital roles in animal genetics and breeding. Large amounts of population genetic variation resources and tools have been developed in human, which provided solid support for human genetic studies. However, compared with human, the development of animal genetic variation databases was relatively slow, which limits the genetic researches in these animals. To fill this gap, we systematically identified ∼ 499 million high-quality SNPs from 4784 samples of 20 types of animals. On that basis, we annotated the functions of SNPs, constructed high-density reference panels and calculated genome-wide linkage disequilibrium (LD) matrixes. We further developed Animal-SNPAtlas, a user-friendly database (http://gong_lab.hzau.edu.cn/Animal_SNPAtlas/) which includes high-quality SNP datasets and several support tools for multiple animals. In Animal-SNPAtlas, users can search the functional annotation of SNPs, perform online genotype imputation, explore and visualize LD information, browse variant information using the genome browser and download SNP datasets for each species. With the massive SNP datasets and useful tools, Animal-SNPAtlas will be an important fundamental resource for the animal genomics, genetics and breeding community.

Giant pandas' staple food bamboo phyllosphere fungal community and its influencing factors

Giant pandas have developed a series of foraging strategies to adapt to their special bamboo diets. Although bamboo is an important food resource for giant pandas in Liziping National Nature Reserve (Liziping NR), China, there are relatively few studies on their phyllosphere fungal community and its influencing factors. Herein, we used ITS1 amplification and metagenomic sequencing to analyze the phyllosphere fungi diversity and functions (KEGG, CAZyme, and antibiotic resistance gene) and explore the influencing factors for the three giant pandas foraging bamboo species (Arundinaria spanostachya, AS; Yushania lineolate, YL; and Fargesia ferax, FF) over different seasons (spring vs. autumn) in Liziping NR, China. We found that Ascomycota and Basidiomycota were the most dominant phyla in the bamboo phyllosphere. The alpha diversity (e.g., the Sobs index and Shannon index) was relatively higher in autumn samples than in spring samples, and the community structure differed significantly between the three bamboo species in spring and autumn. Some biotic and abiotic variables (e.g., the elevation and mean base diameter of bamboo) significantly influenced the abundance, diversity, and community structure of the bamboo phyllosphere fungal community. Moreover, the functional analysis showed the differences in the glycoside hydrolase community and antibiotic resistance gene (ARG) profile between spring and autumn samples. Co-occurrence network modeling suggested that AS phyllosphere fungal communities in autumn employed a much more complex network than that in spring, and the abundance of multidrug, tetracycline, and glycopeptide resistance genes was high and closely correlated with other ARGs. These results indicate that fungal community's abundance, diversity, and community structure are mainly affected by the season, host species, and elevation. The season and host species are major factors affecting the biological functions (KEGG and CAZyme), ARGs, and interactions between sympatric bacterial and fungal communities in bamboo phyllosphere. This integrated study can provide a reference basis for the seasonal management of bamboo resources foraged by wild giant pandas, and predict the risk of antibiotic resistance in bamboo phyllosphere fungal flora in Liziping NR (Xiaoxiangling mountains), China.

Comparative transcriptomics of high-altitude Vulpes and their low-altitude relatives

The harsh environment of Qinghai-Tibet Plateau (QTP) imposes strong selective stresses (e.g., hypoxia, high UV-radiation, and extreme temperature) to the native species, which have driven striking phenotypic and genetic adaptations. Although the mechanisms of high-altitude adaptation have been explored for many plateau species, how the phylogenetic background contributes to genetic adaption to high-altitude of Vulpes is largely unknown. In this study, we sequenced transcriptomic data across multiple tissues of two high-altitude Vulpes (Vulpes vulpes montana and Vulpes ferrilata) and their low-altitude relatives (Vulpes corsac and Vulpes lagopus) to search the genetic and gene expression changes caused by high-altitude environment. The results indicated that the positive selection genes (PSGs) identified by both high-altitude Vulpes are related to angiogenesis, suggesting that angiogenesis may be the result of convergent evolution of Vulpes in the face of hypoxic selection pressure. In addition, more PSGs were detected in V. ferrilata than in V. v. montana, which may be related to the longer adaptation time of V. ferrilata to plateau environment and thus more genetic changes. Besides, more PSGs associated with high-altitude adaptation were identified in V. ferrilata compared with V. v. montana, indicating that the longer the adaptation time to the high-altitude environment, the more genetic alterations of the species. Furthermore, the result of expression profiles revealed a tissue-specific pattern between Vulpes. We also observed that differential expressed genes in the high-altitude group exhibited species-specific expression patterns, revealed a convergent expression pattern of Vulpes in high-altitude environment. In general, our research provides a valuable transcriptomic resource for further studies, and expands our understanding of high-altitude adaptation within a phylogenetic context.

Ursids evolved early and continuously to be low-protein macronutrient omnivores

The eight species of bears world-wide consume a wide variety of diets. Some are specialists with extensive anatomical and physiological adaptations necessary to exploit specific foods or environments [e.g., polar bears (Ursus maritimus), giant pandas (Ailuropoda melanoleuca), and sloth bears (Melursus ursinus)], while the rest are generalists. Even though ursids evolved from a high-protein carnivore, we hypothesized that all have become low-protein macronutrient omnivores. While this dietary strategy has already been described for polar bears and brown bears (Ursus arctos), a recent study on giant pandas suggested their macronutrient selection was that of the ancestral high-protein carnivore. Consumption of diets with inappropriate macronutrient profiles has been associated with increased energy expenditure, ill health, failed reproduction, and premature death. Consequently, we conducted feeding and preference trials with giant pandas and sloth bears, a termite and ant-feeding specialist. Both giant pandas and sloth bears branched off from the ursid lineage a million or more years before polar bears and brown bears. We found that giant pandas are low-protein, high-carbohydrate omnivores, whereas sloth bears are low-protein, high-fat omnivores. The preference for low protein diets apparently occurred early in the evolution of ursids and may have been critical to their world-wide spread.