Convergent gene losses illuminate metabolic and physiological changes in herbivores and carnivores

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.

Genetic architecture of long-distance migration and population genomics of the endangered Japanese eel

The Japanese eel (Anguilla japonica), a flagship anguillid species for conservation, is known for its long-distance-oriented migration. However, our understanding of the genetic architecture underlying long-distance migration and population genomic characteristics of A. japonica is still limited. Here, we generated a high-quality chromosome-level genome assembly and conducted whole-genome resequencing of 218 individuals to explore these aspects. Strong signals of selection were found on genes involved in long-distance aerobic exercise and navigation, which might be associated with evolutionary adaptation to long-distance migrations. Low genetic diversity was detected, which might result from genetic drift associated with demographic declines. Both mitochondrial and nuclear genomic datasets supported the existence of a single panmictic population for Japanese eel, despite signals of single-generation selection. Candidate genes for local selection involved in functions like development and circadian rhythm. The findings can provide insights to adaptative evolution to long-distance migration and inform conservation efforts for A. japonica.

Dissection of the synthesis of polyunsaturated fatty acids in nematodes and Collembola of the soil fauna

It is becoming increasingly clear that not only unicellular, photoautotrophic eukaryotes, plants, and fungi, but also invertebrates are capable of synthesizing ω3 long-chain polyunsaturated fatty acids (LC-PUFA) de novo. However, the distribution of this anabolic capacity among different invertebrate groups and its implementation at the gene and protein level are often still unknown. This study investigated the PUFA pathways in common soil fauna, i.e. two nematode and two Collembola species. Of these, one species each (Panagrellus redivivus, Folsomia candida) was assumed to produce ω3 LC-PUFA de novo, while the others (Acrobeloides bodenheimeri, Isotoma caerulea) were supposed to be unable to do so. A highly labeled oleic acid (99 % 13C) was supplemented and the isotopic signal was used to trace its metabolic path. All species followed the main pathway of lipid biosynthesis. However, in A. bodenheimeri this terminated at arachidonic acid (ω6 PUFA), whereas the other three species continued the pathway to eicosapentaenoic acid (ω3 PUFA), including I. caerulea. For the nematode P. redivivus the identification and functional characterization of four new fatty acid desaturase (FAD) genes was performed. These genes encode the FAD activities Δ9, Δ6, and Δ5, respectively. Additionally, the Δ12 desaturase was analyzed, yet the observed activity of an ω3 FAD could not be attributed to a coding gene. In the Collembola F. candida, 11 potential first desaturases (Δ9) and 13 front-end desaturases (Δ6 or Δ5 FADs) have been found. Further sequence analysis indicates the presence of omega FADs, specifically Δ12, which are likely derived from Δ9 FADs.

Dietary type (carnivore, herbivore and omnivore) and animal species modulate the nutritional metabolome of terrestrial species

Ecometabolomics could be implemented as a powerful tool in molecular ecology studies, but it is necessary to know the baseline of certain metabolites and understand how different traits could affect the metabolome of the animals. Therefore, the main objective of this study was to provide values for the nutritional metabolome profile of different diet groups and animal species, as well as to study the differences in the metabolomic profile due to the effect of diet type and species. To achieve this goal, blood samples were taken from healthy animals (n = 43) of different species: lion (Panthera leo), jaguar (Panthera onca), chimpanzee (Pan troglodytes), bison (Bison bison), gazelle (Gazella cuvieri) and fallow deer (Dama dama), and with different types of diet (carnivore, herbivore and omnivore). Each blood sample was analysed to determine nutritional metabolites. The main results this study provides are the nutritional metabolic profile of these animals based on the type of diet and the animal species. A significant effect of the dietary type was found on nutritional metabolite levels, with those metabolites related to protein metabolism (total protein and creatine) being higher in carnivores. There is also an effect of the species on nutritional metabolites, observing a metabolome differentiation between lion and jaguar. In the case of herbivores, bison showed higher levels of uric acid and cholesterol, and lower urea levels than gazelle and fallow deer. More molecular ecology studies are needed to further the knowledge of the metabolism of these animals.

The Male Reproductive Toxicity Caused by 2-Naphthylamine Was Related to Testicular Immunity Disorders

2-naphthylamine (NAP) was classified as a group I carcinogen associated with bladder cancer. The daily exposure is mostly from cigarette and E-cigarette smoke. NAP can lead to testicular atrophy and interstitial tissue hyperplasia; however, the outcomes of NAP treatment on spermatogenesis and the associated mechanisms have not been reported. The study aimed to investigate the effect of NAP on spermatogenesis and sperm physiologic functions after being persistently exposed to NAP at 5, 20, and 40 mg/kg for 35 days. We found that sperm motility, progressive motility, sperm average path velocity, and straight-line velocity declined remarkably in the NAP (40 mg/kg) treated group, and the sperm deformation rate rose upon NAP administration. The testis immunity- and lipid metabolism-associated processes were enriched from RNA-sequence profiling. Plvap, Ccr7, Foxn1, Trim29, Sirpb1c, Cfd, and Lpar4 involved in testis immunity and Pnliprp1 that inhibit triglyceride and cholesterol absorption were confirmed to rise dramatically in the NAP-exposed group. The increased total cholesterol and CD68 levels were observed in the testis from the NAP-exposed group. Gpx5, serving as an antioxidant in sperm plasma, and Semg1, which contributes to sperm progressive motility, were both down-regulated. We concluded that the short-term exposure to NAP caused reproductive toxicity, primarily due to the inflammatory abnormality in the testis.

In-depth correlation analysis between tear glucose and blood glucose using a wireless smart contact lens

Tears have emerged as a promising alternative to blood for diagnosing diabetes. Despite increasing attempts to measure tear glucose using smart contact lenses, the controversy surrounding the correlation between tear glucose and blood glucose still limits the clinical usage of tears. Herein, we present an in-depth investigation of the correlation between tear glucose and blood glucose using a wireless and soft smart contact lens for continuous monitoring of tear glucose. This smart contact lens is capable of quantitatively monitoring the tear glucose levels in basal tears excluding the effect of reflex tears which might weaken the relationship with blood glucose. Furthermore, this smart contact lens can provide an unprecedented level of continuous tear glucose data acquisition at sub-minute intervals. These advantages allow the precise estimation of lag time, enabling the establishment of the concept called 'personalized lag time'. This demonstration considers individual differences and is successfully applied to both non-diabetic and diabetic humans, as well as in animal models, resulting in a high correlation.

Functional analysis of CYP71AV1 reveals the evolutionary landscape of artemisinin biosynthesis

Artemisinin biosynthesis, unique to Artemisia annua, is suggested to have evolved from the ancestral costunolide biosynthetic pathway commonly found in the Asteraceae family. However, the evolutionary landscape of this process is not fully understood. The first oxidase in artemisinin biosynthesis, CYP71AV1, also known as amorpha-4,11-diene oxidase (AMO), has specialized from ancestral germacrene A oxidases (GAOs). Unlike GAO, which exhibits catalytic promiscuity toward amorpha-4,11-diene, the natural substrate of AMO, AMO has lost its ancestral activity on germacrene A. Previous studies have suggested that the loss of the GAO copy in A. annua is responsible for the abolishment of the costunolide pathway. In the genome of A. annua, there are two copies of AMO, each of which has been reported to be responsible for the different product profiles of high- and low-artemisinin production chemotypes. Through analysis of their tissue-specific expression and comparison of their sequences with those of other GAOs, it was discovered that one copy of AMO (AMOHAP) exhibits a different transcript compared to the reported artemisinin biosynthetic genes and shows more sequence similarity to other GAOs in the catalytic regions. Furthermore, in a subsequent in vitro enzymatic assay, the recombinant protein of AMOHAP unequivocally demonstrated GAO activity. This result clearly indicates that AMOHAP is a GAO rather than an AMO and that its promiscuous activity on amorpha-4,11-diene has led to its misidentification as an AMO in previous studies. In addition, the divergent expression pattern of AMOHAP compared to that of the upstream germacrene A synthase may have contributed to the abolishment of costunolide biosynthesis in A. annua. Our findings reveal a complex evolutionary landscape in which the emergence of a new metabolic pathway replaces an ancestral one.

Co-option of an Astacin Metalloprotease Is Associated with an Evolutionarily Novel Feeding Morphology in a Predatory Nematode

Animals consume a wide variety of food sources to adapt to different environments. However, the genetic mechanisms underlying the acquisition of evolutionarily novel feeding morphology remain largely unknown. While the nematode Caenorhabditis elegans feeds on bacteria, the satellite species Pristionchus pacificus exhibits predatory feeding behavior toward other nematodes, which is an evolutionarily novel feeding habit. Here, we found that the astacin metalloprotease Ppa-NAS-6 is required for the predatory killing by P. pacificus. Ppa-nas-6 mutants were defective in predation-associated characteristics, specifically the tooth morphogenesis and tooth movement during predation. Comparison of expression patterns and rescue experiments of nas-6 in P. pacificus and C. elegans suggested that alteration of the spatial expression patterns of NAS-6 may be vital for acquiring predation-related traits. Reporter analysis of the Ppa-nas-6 promoter in C. elegans revealed that the alteration in expression patterns was caused by evolutionary changes in cis- and trans-regulatory elements. This study suggests that the co-option of a metalloprotease is involved in an evolutionarily novel feeding morphology.

Duplications of Human Longevity-Associated Genes Across Placental Mammals

Natural selection has shaped a wide range of lifespans across mammals, with a few long-lived species showing negligible signs of ageing. Approaches used to elucidate the genetic mechanisms underlying mammalian longevity usually involve phylogenetic selection tests on candidate genes, detections of convergent amino acid changes in long-lived lineages, analyses of differential gene expression between age cohorts or species, and measurements of age-related epigenetic changes. However, the link between gene duplication and evolution of mammalian longevity has not been widely investigated. Here, we explored the association between gene duplication and mammalian lifespan by analyzing 287 human longevity-associated genes across 37 placental mammals. We estimated that the expansion rate of these genes is eight times higher than their contraction rate across these 37 species. Using phylogenetic approaches, we identified 43 genes whose duplication levels are significantly correlated with longevity quotients (False Discovery Rate (FDR) < 0.05). In particular, the strong correlation observed for four genes (CREBBP, PIK3R1, HELLS, FOXM1) appears to be driven mainly by their high duplication levels in two ageing extremists, the naked mole rat (Heterocephalus glaber) and the greater mouse-eared bat (Myotis myotis). Further sequence and expression analyses suggest that the gene PIK3R1 may have undergone a convergent duplication event, whereby the similar region of its coding sequence was independently duplicated multiple times in both of these long-lived species. Collectively, this study identified several candidate genes whose duplications may underlie the extreme longevity in mammals, and highlighted the potential role of gene duplication in the evolution of mammalian long lifespans.

Genome-wide analysis of cuticle protein family genes in rice stem borer Chilo suppressalis: Insights into their role in environmental adaptation and insecticidal stress response

Insect cuticle plays a key role in insect survival, adaptation and prosperity by serving as the exoskeleton and the first barrier against environmental stresses. As the major components of insect cuticle, the diverse structural cuticle proteins (CPs) contribute to variation in physical properties and functions of cuticle. However, the roles of CPs in cuticular versatility, especially in the stress response or adaption, remain incompletely understood. In this study, we performed a genome-wide analysis of CP superfamily in the rice-boring pest Chilo suppressalis. A total of 211 CP genes were identified and their encoding proteins were classified into eleven families and three subfamilies (RR1, RR2, and RR3). The comparative genomic analysis of CPs revealed that C. suppressalis had fewer CP genes compared to other lepidopteran species, which largely resulted from a less expansion of his-rich RR2 genes involved in cuticular sclerotization, suggesting long-term boring life of C. suppressalis inside rice hosts might evolutionarily prefer cuticular elasticity rather than cuticular sclerotization. We also investigated the response pattern of all CP genes under insecticidal stresses. >50 % CsCPs were upregulated at least 2-fold under insecticidal stresses. Notably, the majority of the highly upregulated CsCPs formed gene pairs or gene clusters on chromosomes, indicating the rapid response of adjacent CsCPs to insecticidal stress. Most high-response CsCPs encoded AAPA/V/L motifs that are related to cuticular elasticity and >50 % of the sclerotization-related his-rich RR2 genes were also upregulated. These results suggested the potential roles of CsCPs in balancing the elasticity and sclerotization of cuticles, which is essential for the survival and adaptation of plant borers including C. suppressalis. Our study provides valuable information for further developing cuticle-based strategies of both pest management and biomimetic applications.

Next Generation Sequencing Revolutionizes Organismal Biology Research in Bats

The advent of next generation sequencing technologies (NGS) has greatly accelerated our understanding of critical aspects of organismal biology from non-model organisms. Bats form a particularly interesting group in this regard, as genomic data have helped unearth a vast spectrum of idiosyncrasies in bat genomes associated with bat biology, physiology, and evolution. Bats are important bioindicators and are keystone species to many eco-systems. They often live in proximity to humans and are frequently associated with emerging infectious diseases, including the COVID-19 pandemic. Nearly four dozen bat genomes have been published to date, ranging from drafts to chromosomal level assemblies. Genomic investigations in bats have also become critical towards our understanding of disease biology and host-pathogen coevolution. In addition to whole genome sequencing, low coverage genomic data like reduced representation libraries, resequencing data, etc. have contributed significantly towards our understanding of the evolution of natural populations, and their responses to climatic and anthropogenic perturbations. In this review, we discuss how genomic data have enhanced our understanding of physiological adaptations in bats (particularly related to ageing, immunity, diet, etc.), pathogen discovery, and host pathogen co-evolution. In comparison, the application of NGS towards population genomics, conservation, biodiversity assessment, and functional genomics has been appreciably slower. We reviewed the current areas of focus, identifying emerging topical research directions and providing a roadmap for future genomic studies in bats.

Reduction of Paraoxonase Expression Followed by Inactivation across Independent Semiaquatic Mammals Suggests Stepwise Path to Pseudogenization

Convergent adaptation to the same environment by multiple lineages frequently involves rapid evolutionary change at the same genes, implicating these genes as important for environmental adaptation. Such adaptive molecular changes may yield either change or loss of protein function; loss of function can eliminate newly deleterious proteins or reduce energy necessary for protein production. We previously found a striking case of recurrent pseudogenization of the Paraoxonase 1 (Pon1) gene among aquatic mammal lineages-Pon1 became a pseudogene with genetic lesions, such as stop codons and frameshifts, at least four times independently in aquatic and semiaquatic mammals. Here, we assess the landscape and pace of pseudogenization by studying Pon1 sequences, expression levels, and enzymatic activity across four aquatic and semiaquatic mammal lineages: pinnipeds, cetaceans, otters, and beavers. We observe in beavers and pinnipeds an unexpected reduction in expression of Pon3, a paralog with similar expression patterns but different substrate preferences. Ultimately, in all lineages with aquatic/semiaquatic members, we find that preceding any coding-level pseudogenization events in Pon1, there is a drastic decrease in expression, followed by relaxed selection, thus allowing accumulation of disrupting mutations. The recurrent loss of Pon1 function in aquatic/semiaquatic lineages is consistent with a benefit to Pon1 functional loss in aquatic environments. Accordingly, we examine diving and dietary traits across pinniped species as potential driving forces of Pon1 functional loss. We find that loss is best associated with diving activity and likely results from changes in selective pressures associated with hypoxia and hypoxia-induced inflammation.

Integrating gene annotation with orthology inference at scale

Annotating coding genes and inferring orthologs are two classical challenges in genomics and evolutionary biology that have traditionally been approached separately, limiting scalability. We present TOGA (Tool to infer Orthologs from Genome Alignments), a method that integrates structural gene annotation and orthology inference. TOGA implements a different paradigm to infer orthologous loci, improves ortholog detection and annotation of conserved genes compared with state-of-the-art methods, and handles even highly fragmented assemblies. TOGA scales to hundreds of genomes, which we demonstrate by applying it to 488 placental mammal and 501 bird assemblies, creating the largest comparative gene resources so far. Additionally, TOGA detects gene losses, enables selection screens, and automatically provides a superior measure of mammalian genome quality. TOGA is a powerful and scalable method to annotate and compare genes in the genomic era.

Establishment of killer whale (Orcinus orca) primary fibroblast cell cultures and their transcriptomic responses to pollutant exposure

Populations of killer whale (Orcinus orca) contain some of the most polluted animals on Earth. Yet, the knowledge on effects of chemical pollutants is limited in this species. Cell cultures and in vitro exposure experiments are pertinent tools to study effects of pollutants in free-ranging marine mammals. To investigate transcriptional responses to pollutants in killer whale cells, we collected skin biopsies of killer whales from the Northern Norwegian fjords and successfully established primary fibroblast cell cultures from the dermis of 4 out of 5 of them. Cells from the individual with the highest cell yield were exposed to three different concentrations of a mixture of persistent organic pollutants (POPs) that reflects the composition of the 10 most abundant POPs found in Norwegian killer whales (p,p'-DDE, trans-nonachlor, PCB52, 99, 101, 118, 138, 153, 180, 187). Transcriptional responses of 13 selected target genes were studied using digital droplet PCR, and whole transcriptome responses were investigated utilizing RNA sequencing. Among the target genes analysed, CYP1A1 was significantly downregulated in the cells exposed to medium (11.6 µM) and high (116 µM) concentrations of the pollutant mixture, while seven genes involved in endocrine functions showed a non-significant tendency to be upregulated at the highest exposure concentration. Bioinformatic analyses of RNA-seq data indicated that 13 and 43 genes were differentially expressed in the cells exposed to low and high concentrations of the mixture, respectively, in comparison to solvent control. Subsequent pathway and functional analyses of the differentially expressed genes indicated that the enriched pathways were mainly related to lipid metabolism, myogenesis and glucocorticoid receptor regulation. The current study results support previous correlative studies and provide cause-effect relationships, which is highly relevant for chemical and environmental management.

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.

Adaptive Evolution of the OAS Gene Family Provides New Insights into the Antiviral Ability of Laurasiatherian Mammals

Many mammals risk damage from virus invasion due to frequent environmental changes. The oligoadenylate synthesis (OAS) gene family, which is an important component of the immune system, provides an essential response to the antiviral activities of interferons by regulating immune signal pathways. However, little is known about the evolutionary characteristics of OASs in Laurasiatherian mammals. Here, we examined the evolution of the OAS genes in 64 mammals to explore the accompanying molecular mechanisms of the antiviral ability of Laurasiatherian mammals living in different environments. We found that OAS2 and OAS3 were found to be pseudogenes in Odontoceti species. This may be related to the fact that they live in water. Some Antilopinae, Caprinae, and Cervidae species lacked the OASL gene, which may be related to their habitats being at higher altitudes. The OASs had a high number of positive selection sites in Cetartiodactyla, which drove the expression of strong antiviral ability. The OAS gene family evolved in Laurasiatherian mammals at different rates and was highly correlated with the species' antiviral ability. The gene evolution rate in Cetartiodactyla was significantly higher than that in the other orders. Compared to other species of the Carnivora family, the higher selection pressure on the OAS gene and the absence of positive selection sites in Canidae may be responsible for its weak resistance to rabies virus. The OAS gene family was relatively conserved during evolution. Conserved genes are able to provide better maintenance of gene function. The rate of gene evolution and the number of positively selected sites combine to influence the resistance of a species to viruses. The positive selection sites demonstrate the adaptive evolution of the OAS gene family to the environment. Adaptive evolution combined with conserved gene function improves resistance to viruses. Our findings offer insights into the molecular and functional evolution of the antiviral ability of Laurasian mammals.

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 eggshell-matrix protein gene OC-17 is functionally lost in the viviparous Chinese crocodile lizard

Thickness reduction or loss of the calcareous eggshell is one of major phenotypic changes in the transition from oviparity to viviparity. Whether the reduction of eggshells in viviparous squamates is associated with specific gene losses is unknown. Taking advantage of a newly generated high-quality genome of the viviparous Chinese crocodile lizard (Shinisaurus crocodilurus), we found that ovocleidin-17 gene (OC-17), which encodes an eggshell matrix protein that is essential for calcium deposition in eggshells, is not intact in the crocodile lizard genome. Only OC-17 transcript fragments were found in the oviduct transcriptome, and no OC-17 peptides were identified in the eggshell proteome of crocodile lizards. In contrast, OC-17 was present in the eggshells of the oviparous Mongolia racerunner (Eremias argus). Although the loss of OC-17 is not common in viviparous species, viviparous squamates show fewer intact eggshell-specific proteins than oviparous squamates. Our study implies that functional loss of eggshell-matrix protein genes may be involved in the reduction of eggshells during the transition from oviparity to viviparity in the crocodile lizard.

Ambush predation and the origin of euprimates

Primates of modern aspect (euprimates) are characterized by a suite of characteristics (e.g., convergent orbits, grasping hands and feet, reduced claws, and leaping), but the selective pressures responsible for the evolution of these euprimate characteristics have long remained controversial. Here, we used a molecular phyloecological approach to determine the diet of the common ancestor of living primates (CALP), and the results showed that the CALP had increased carnivory. Given the carnivory of the CALP, along with the general observation that orbital convergence is largely restricted to ambush predators, our study suggests that the euprimate characteristics could have been more specifically adapted for ambush predation. In particular, our behavior experiment further shows that nonclaw climbing can significantly reduce noises, which could benefit the ancestral euprimates' stalking to ambush their prey in trees. Therefore, our study suggests that the distinctive euprimate characteristics may have evolved as their specialized adaptation for ambush predation in arboreal environments.

Vision-related convergent gene losses reveal SERPINE3's unknown role in the eye

Despite decades of research, knowledge about the genes that are important for development and function of the mammalian eye and are involved in human eye disorders remains incomplete. During mammalian evolution, mammals that naturally exhibit poor vision or regressive eye phenotypes have independently lost many eye-related genes. This provides an opportunity to predict novel eye-related genes based on specific evolutionary gene loss signatures. Building on these observations, we performed a genome-wide screen across 49 mammals for functionally uncharacterized genes that are preferentially lost in species exhibiting lower visual acuity values. The screen uncovered several genes, including SERPINE3, a putative serine proteinase inhibitor. A detailed investigation of 381 additional mammals revealed that SERPINE3 is independently lost in 18 lineages that typically do not primarily rely on vision, predicting a vision-related function for this gene. To test this, we show that SERPINE3 has the highest expression in eyes of zebrafish and mouse. In the zebrafish retina, serpine3 is expressed in Müller glia cells, a cell type essential for survival and maintenance of the retina. A CRISPR-mediated knockout of serpine3 in zebrafish resulted in alterations in eye shape and defects in retinal layering. Furthermore, two human polymorphisms that are in linkage with SERPINE3 are associated with eye-related traits. Together, these results suggest that SERPINE3 has a role in vertebrate eyes. More generally, by integrating comparative genomics with experiments in model organisms, we show that screens for specific phenotype-associated gene signatures can predict functions of uncharacterized genes.

Gene losses in the common vampire bat illuminate molecular adaptations to blood feeding

Vampire bats are the only mammals that feed exclusively on blood. To uncover genomic changes associated with this dietary adaptation, we generated a haplotype-resolved genome of the common vampire bat and screened 27 bat species for genes that were specifically lost in the vampire bat lineage. We found previously unknown gene losses that relate to reduced insulin secretion (FFAR1 and SLC30A8), limited glycogen stores (PPP1R3E), and a unique gastric physiology (CTSE). Other gene losses likely reflect the biased nutrient composition (ERN2 and CTRL) and distinct pathogen diversity of blood (RNASE7) and predict the complete lack of cone-based vision in these strictly nocturnal bats (PDE6H and PDE6C). Notably, REP15 loss likely helped vampire bats adapt to high dietary iron levels by enhancing iron excretion, and the loss of CYP39A1 could have contributed to their exceptional cognitive abilities. These findings enhance our understanding of vampire bat biology and the genomic underpinnings of adaptations to blood feeding.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes, as well as searching for species-specific structural variants (SVs). Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observe species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This strongly argues for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

Gene losses may contribute to subterranean adaptations in naked mole-rat and blind mole-rat

Background

Naked mole-rats (Heterocephalus glaber, NMRs) and blind mole-rats (Spalax galili, BMRs) are representative subterranean rodents that have evolved many extraordinary traits, including hypoxia tolerance, longevity, and cancer resistance. Although multiple candidate loci responsible for these traits have been uncovered by genomic studies, many of them are limited to functional changes to amino acid sequence and little is known about the contributions of other genetic events. To address this issue, we focused on gene losses (unitary pseudogenes) and systematically analyzed gene losses in NMRs and BMRs, aiming to elucidate the potential roles of pseudogenes in their adaptation to subterranean lifestyle.

Results

We obtained the pseudogene repertoires in NMRs and BMRs, as well as their respective aboveground relatives, guinea pigs and rats, on a genome-wide scale. As a result, 167, 139, 341, and 112 pseudogenes were identified in NMRs, BMRs, guinea pigs, and rats, respectively. Functional enrichment analysis identified 4 shared and 2 species-specific enriched functional groups (EFGs) in subterranean lineages. Notably, the pseudogenes in these EFGs might be associated with either regressive (e.g., visual system) or adaptive (e.g., altered DNA damage response) traits. In addition, several pseudogenes including TNNI3K and PDE5A might be associated with specific cardiac features observed in subterranean lineages. Interestingly, we observed 20 convergent gene losses in NMRs and BMRs. Given that the functional investigations of these genes are generally scarce, we provided functional evidence that independent loss of TRIM17 in NMRs and BMRs might be beneficial for neuronal survival under hypoxia, supporting the positive role of eliminating TRIM17 function in hypoxia adaptation. Our results also suggested that pseudogenes, together with positively selected genes, reinforced subterranean adaptations cooperatively.

Conclusions

Our study provides new insights into the molecular underpinnings of subterranean adaptations and highlights the importance of gene losses in mammalian evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01243-0.

Proteomic Analysis of the Intestinal Resistance to Thyroid Hormone Mouse Model With Thyroid Hormone Receptor Alpha Mutations

Thyroid hormone is critical during the development of vertebrates and affects the function of many organs and tissues, especially the intestine. Triiodothyronine (T₃) is the active form and can bind to thyroid hormone nuclear receptors (TRs) to play a vital role in the development of vertebrates. The resistance to thyroid hormone α, as seen in patients, has been mimicked by the Thra^E403X mutation. To investigate the mechanisms underlying the effect of TRα1 on intestinal development, the present study employed proteomic analysis to identify differentially expressed proteins (DEPs) in the distal ileum between homozygous Thra^E403X/E403X and wild-type Thra^+/+ mice. A total of 1,189 DEPs were identified, including 603 upregulated and 586 downregulated proteins. Proteomic analysis revealed that the DEPs were highly enriched in the metabolic process, the developmental process, the transporter of the nutrients, and the intestinal immune system-related pathway. Of these DEPs, 20 proteins were validated by parallel reaction monitoring analysis. Our intestinal proteomic results provide promising candidates for future studies, as they suggest novel mechanisms by which TRα1 may influence intestinal development, such as the transport of intestinal nutrients and the establishment of innate and adaptive immune barriers of the intestine.

Carotenoid modifying enzymes in metazoans

Carotenoids constitute an essential dietary component of animals and other non-carotenogenic species which use these pigments in both their modified and unmodified forms. Animals utilize uncleaved carotenoids to mitigate light damage and oxidative stress and to signal fitness and health. Carotenoids also serve as precursors of apocarotenoids including retinol, and its retinoid metabolites, which carry out essential functions in animals by forming the visual chromophore 11-cis-retinaldehyde. Retinoids, such as all-trans-retinoic acid, can also act as ligands of nuclear hormone receptors. The fact that enzymes and biochemical pathways responsible for the metabolism of carotenoids in animals bear resemblance to the ones in plants and other carotenogenic species suggests an evolutionary relationship. We will explore some of the modes of transmission of carotenoid genes from carotenogenic species to metazoans. This apparent relationship has been successfully exploited in the past to identify and characterize new carotenoid and retinoid modifying enzymes. We will review approaches used to identify putative animal carotenoid enzymes, and we will describe methods used to functionally validate and analyze the biochemistry of carotenoid modifying enzymes encoded by animals.

Reconstruction of evolutionary changes in fat and toxin consumption reveals associations with gene losses in mammals: a case study for the lipase inhibitor PNLIPRP1 and the xenobiotic receptor NR1I3

The inactivation of ancestral protein-coding genes (gene loss) can be associated with phenotypic modifications. Within placental mammals, repeated losses of PNLIPRP1 (gene inhibiting fat digestion) occurred preferentially in strictly herbivorous species, while repeated NR1I3 losses (gene involved in detoxification) occurred preferentially in strictly carnivorous species. It was hypothesized that lower fat contents of herbivorous diets and lower toxin contents of carnivorous diets cause relaxed selection pressure on these genes resulting in the accumulation of mutations and ultimately to convergent gene losses. However, since herbivorous and carnivorous diets differ vastly in their composition, a fine-grained analysis is required for hypothesis testing. We generated a trait matrix recording diet and semi-quantitative estimates of fat and toxin consumption for 52 placental species. By including data from 31 fossil taxa, we reconstructed the ancestral diets in major lineages (grundplan reconstruction). We found support that PNLIPRP1 loss is primarily associated with low levels of fat intake and not simply with herbivory/carnivory. In particular, PNLIPRP1 loss also occurred in carnivorous lineages feeding on a fat-poor diet, suggesting that the loss of this gene may be beneficial for occupying ecological niches characterized by fat-poor food resources. Similarly, we demonstrated that carnivorous species are indeed less exposed to diet-related toxins suggesting that the loss of NR1I3 and related genes (NR1I2, UGT1A6) resulted from relaxed selection pressure. This study illustrates the need of detailed phenotype studies to obtain a deeper understanding of factors underlying gene losses and to progress in understanding genomic causes of phenotypic variation in mammals.

Phenotyping in the era of genomics: MaTrics—a digital character matrix to document mammalian phenotypic traits

A new and uniquely structured matrix of mammalian phenotypes, MaTrics (Mammalian Traits for Comparative Genomics) in a digital form is presented. By focussing on mammalian species for which genome assemblies are available, MaTrics provides an interface between mammalogy and comparative genomics.

MaTrics was developed within a project aimed to find genetic causes of phenotypic traits of mammals using Forward Genomics. This approach requires genomes and comprehensive and recorded information on homologous phenotypes that are coded as discrete categories in a matrix. MaTrics is an evolving online resource providing information on phenotypic traits in numeric code; traits are coded either as absent/present or with several states as multistate. The state record for each species is linked to at least one reference (e.g., literature, photographs, histological sections, CT scans, or museum specimens) and so MaTrics contributes to digitalization of museum collections. Currently, MaTrics covers 147 mammalian species and includes 231 characters related to structure, morphology, physiology, ecology, and ethology and available in a machine actionable NEXUS-format*. Filling MaTrics revealed substantial knowledge gaps, highlighting the need for phenotyping efforts. Studies based on selected data from MaTrics and using Forward Genomics identified associations between genes and certain phenotypes ranging from lifestyles (e.g., aquatic) to dietary specializations (e.g., herbivory, carnivory). These findings motivate the expansion of phenotyping in MaTrics by filling research gaps and by adding taxa and traits. Only databases like MaTrics will provide machine actionable information on phenotypic traits, an important limitation to genomics. MaTrics is available within the data repository Morph·D·Base (www.morphdbase.de).

Noninsect-Based Diet Leads to Structural and Functional Changes of Acidic Chitinase in Carnivora

Acidic chitinase (Chia) digests the chitin of insects in the omnivorous stomach and the chitinase activity in carnivorous Chia is significantly lower than that of the omnivorous enzyme. However, mechanistic and evolutionary insights into the functional changes in Chia remain unclear. Here we show that a noninsect-based diet has caused structural and functional changes in Chia during the course of evolution in Carnivora. By creating mouse-dog chimeric Chia proteins and modifying the amino acid sequences, we revealed that F214L and A216G substitutions led to the dog enzyme activation. In 31 Carnivora, Chia was present as a pseudogene with stop codons in the open reading frame (ORF) region. Importantly, the Chia proteins of skunk, meerkat, mongoose, and hyena, which are insect-eating species, showed high chitinolytic activity. The cat Chia pseudogene product was still inactive even after ORF restoration. However, the enzyme was activated by matching the number and position of Cys residues to an active form and by introducing five meerkat Chia residues. Mutations affecting the Chia conformation and activity after pseudogenization have accumulated in the common ancestor of Felidae due to functional constraints. Evolutionary analysis indicates that Chia genes are under relaxed selective constraint in species with noninsect-based diets except for Canidae. These results suggest that there are two types of inactivating processes in Carnivora and that dietary changes affect the structure and activity of Chia.

Comparative transcriptome analysis revealed omnivorous adaptation of the small intestine of Melinae

As the main digestive organ, the small intestine plays a vital role in the digestion of animals. At present, most of the research on animal feeding habits focuses on carnivores and herbivores. However, the mechanism of feeding and digestion in omnivores remains unclear. This study aims to reveal the molecular basis of the omnivorous adaptive evolution of Melinae by comparing the transcriptome of the small intestines of Asian Badgers (Meles leucurus) and Northern Hog Badgers (Arctonyx albogularis). We obtained high-quality small intestinal transcriptome data from these two species. Key genes and signalling pathways were analysed through Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and other databases. Research has mainly found that orthologous genes related to six enzymes have undergone adaptive evolution. In addition, the study also found three digestion-related pathways (cGMP-PKG, cAMP, and Hippo). They are related to the digestion and absorption of nutrients, the secretion of intestinal fluids, and the transport of food through the small intestine, which may help omnivorous animals adapt to an omnivorous diet. Our study provides insight into the adaptation of Melinae to omnivores and affords a valuable transcriptome resource for future research.

Creating an atlas to visualize the biodiversity of the mammalian gut microbiota

Given the fundamental role played by the intestinal microbial community in influencing host health, it is not surprising that recent decades have been marked by increased efforts to determine the taxonomic composition of the human gut microbiota and its associated functions. Despite their generally accepted importance, these large-scale human-centered studies prevent an exhaustive overview of those mechanisms and factors that contribute to the mammalian gut microbiota assembly. However, Next-Generation Sequencing techniques and associated bioinformatic tools provide an exciting opportunity to rapidly expand our knowledge on the intestinal microbial communities associated with members of the Mammalia class. These non-human-focused studies established that dietary, host phylogeny, host physiology and anthropogenic influences represent the main factors driving the selection of a specific gut microbial consortium in mammals. The current review is aimed at providing a comprehensive overview on the impact that the above-mentioned factors exert on the assembly of the mammalian gut microbiota.

Variation in the Evolution and Sequences of Proglucagon and the Receptors for Proglucagon-Derived Peptides in Mammals

The mammalian proglucagon gene (Gcg) encodes three glucagon like sequences, glucagon, glucagon-like peptide-1 (GLP-1), and glucagon-like peptide-2 that are of similar length and share sequence similarity, with these hormones having cell surface receptors, glucagon receptor (Gcgr), GLP-1 receptor (Glp1r), and GLP-2 receptor (Glp2r), respectively. Gcgr, Glp1r, and Glp2r are all class B1 G protein-coupled receptors (GPCRs). Despite their sequence and structural similarity, analyses of sequences from rodents have found differences in patterns of sequence conservation and evolution. To determine whether these were rodent-specific traits or general features of these genes in mammals I analyzed coding and protein sequences for proglucagon and the receptors for proglucagon-derived peptides from the genomes of 168 mammalian species. Single copy genes for each gene were found in almost all genomes. In addition to glucagon sequences within Hystricognath rodents (e.g., guinea pig), glucagon sequences from a few other groups (e.g., pangolins and some bats) as well as changes in the proteolytic processing of GLP-1 in some bats are suggested to have functional effects. GLP-2 sequences display increased variability but accepted few substitutions that are predicted to have functional consequences. In parallel, Glp2r sequences display the most rapid protein sequence evolution, and show greater variability in amino acids at sites involved in ligand interaction, however most were not predicted to have a functional consequence. These observations suggest that a greater diversity in biological functions for proglucagon-derived peptides might exist in mammals.

Lagomorpha as a Model Morphological System

Due to their global distribution, invasive history, and unique characteristics, European rabbits are recognizable almost anywhere on our planet. Although they are members of a much larger group of living and extinct mammals [Mammalia, Lagomorpha (rabbits, hares, and pikas)], the group is often characterized by several well-known genera (e.g., Oryctolagus, Sylvilagus, Lepus, and Ochotona). This representation does not capture the extraordinary diversity of behavior and form found throughout the order. Model organisms are commonly used as exemplars for biological research, but there are a limited number of model clades or lineages that have been used to study evolutionary morphology in a more explicitly comparative way. We present this review paper to show that lagomorphs are a strong system in which to study macro- and micro-scale patterns of morphological change within a clade that offers underappreciated levels of diversity. To this end, we offer a summary of the status of relevant aspects of lagomorph biology.

Evolutionary and phylogenetic insights from a nuclear genome sequence of the extinct, giant, "subfossil" koala lemur Megaladapis edwardsi

No endemic Madagascar animal with body mass >10 kg survived a relatively recent wave of extinction on the island. From morphological and isotopic analyses of skeletal "subfossil" remains we can reconstruct some of the biology and behavioral ecology of giant lemurs (primates; up to ∼160 kg) and other extraordinary Malagasy megafauna that survived into the past millennium. Yet, much about the evolutionary biology of these now-extinct species remains unknown, along with persistent phylogenetic uncertainty in some cases. Thankfully, despite the challenges of DNA preservation in tropical and subtropical environments, technical advances have enabled the recovery of ancient DNA from some Malagasy subfossil specimens. Here, we present a nuclear genome sequence (∼2× coverage) for one of the largest extinct lemurs, the koala lemur Megaladapis edwardsi (∼85 kg). To support the testing of key phylogenetic and evolutionary hypotheses, we also generated high-coverage nuclear genomes for two extant lemurs, Eulemur rufifrons and Lepilemur mustelinus, and we aligned these sequences with previously published genomes for three other extant lemurs and 47 nonlemur vertebrates. Our phylogenetic results confirm that Megaladapis is most closely related to the extant Lemuridae (typified in our analysis by E. rufifrons) to the exclusion of L. mustelinus, which contradicts morphology-based phylogenies. Our evolutionary analyses identified significant convergent evolution between M. edwardsi and an extant folivore (a colobine monkey) and an herbivore (horse) in genes encoding proteins that function in plant toxin biodegradation and nutrient absorption. These results suggest that koala lemurs were highly adapted to a leaf-based diet, which may also explain their convergent craniodental morphology with the small-bodied folivore Lepilemur.

Recapitulating Evolutionary Divergence in a Single Cis-Regulatory Element Is Sufficient to Cause Expression Changes of the Lens Gene Tdrd7

Mutations in cis-regulatory elements play important roles for phenotypic changes during evolution. Eye degeneration in the blind mole rat (BMR; Nannospalax galili) and other subterranean mammals is significantly associated with widespread divergence of eye regulatory elements, but the effect of these regulatory mutations on eye development and function has not been explored. Here, we investigate the effect of mutations observed in the BMR sequence of a conserved noncoding element upstream of Tdrd7, a pleiotropic gene required for lens development and spermatogenesis. We first show that this conserved element is a transcriptional repressor in lens cells and that the BMR sequence partially lost repressor activity. Next, we recapitulated evolutionary changes in this element by precisely replacing the endogenous regulatory element in a mouse line by the orthologous BMR sequence with CRISPR-Cas9. Strikingly, this repressor replacement caused a more than 2-fold upregulation of Tdrd7 in the developing lens; however, increased mRNA level does not result in a corresponding increase in TDRD7 protein nor an obvious lens phenotype, possibly explained by buffering at the posttranscriptional level. Our results are consistent with eye degeneration in subterranean mammals having a polygenic basis where many small-effect mutations in different eye-regulatory elements collectively contribute to phenotypic differences.

Molecular Evolution of Ecological Specialisation: Genomic Insights from the Diversification of Murine Rodents

Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.

Molecular phyloecology suggests a trophic shift concurrent with the evolution of the first birds

Birds are characterized by evolutionary specializations of both locomotion (e.g., flapping flight) and digestive system (toothless, crop, and gizzard), while the potential selection pressures responsible for these evolutionary specializations remain unclear. Here we used a recently developed molecular phyloecological method to reconstruct the diets of the ancestral archosaur and of the common ancestor of living birds (CALB). Our results suggest a trophic shift from carnivory to herbivory (fruit, seed, and/or nut eater) at the archosaur-to-bird transition. The evolutionary shift of the CALB to herbivory may have essentially made them become a low-level consumer and, consequently, subject to relatively high predation risk from potential predators such as gliding non-avian maniraptorans, from which birds descended. Under the relatively high predation pressure, ancestral birds with gliding capability may have then evolved not only flapping flight as a possible anti-predator strategy against gliding predatory non-avian maniraptorans but also the specialized digestive system as an evolutionary tradeoff of maximizing foraging efficiency and minimizing predation risk. Our results suggest that the powered flight and specialized digestive system of birds may have evolved as a result of their tropic shift-associated predation pressure.

Study of non-metric characters of the skull to determine the epigenetic variability in populations of the European wildcat (Felis silvestris silvestris) and domestic cats (Felis catus)

We studied the variability of non-metric cranial traits, mainly foramina, of European wildcats (Felis silvestris silvestris) and domestic cats (Felis catus) from Germany based on 28 non-metric traits in 211 skulls. The domestic cats were grouped together as a statistical population. The wildcats were divided into two populations: Harz and Hesse, which were further subdivided, based on traffic infrastructure, natural landscape, and in the Harz, on time period. Epigenetic variability, epigenetic distance and the fluctuating asymmetry were calculated to assess genetic variability, possible depressions and population stability. The epigenetic variability Iev of the wildcat groups ranged from 0.27 (Hesse II) to 0.40 (Harz I). The difference in Iev between all specimens from Harz and Hesse respectively was less (Iev = 0.37 Harz and 0.31 Hesse). Compared to other studies these values are not assumed to indicate genetic depression. The epigenetic distance between the wildcat samples is 0.0774 overall, and in each case higher between sub-groups of the Harz and Hesse than between groups within these regions, respectively. The significant epigenetic distance between Harz and Hesse might indicate—at least past formerly—restricted connectivity between these regions. The fluctuating asymmetry for wildcats in total is 11.74% and in the sub-groups it ranges from 8.47 to 16.14%. These values are below 20% are at the lower range known from populations of other mammal species. The use of fluctuating asymmetry had also been discussed critically in its usefulness to assess viability of populations.

Sensitivity of turtles to anticoagulant rodenticides: Risk assessment for green sea turtles (Chelonia mydas) in the Ogasawara Islands and comparison of warfarin sensitivity among turtle species

Although anticoagulant rodenticides (ARs) are effectively used for the control of invasive rodents, nontarget species are also frequently exposed to ARs and secondary poisonings occur widely. However, little data is available on the effects of ARs, especially on marine organisms. To evaluate the effects of ARs on marine wildlife, we chose green sea turtles (Chelonia mydas), which are one of the most common marine organisms around the Ogasawara islands, as our primary study species. The sensitivity of these turtles to ARs was assessed using both in vivo and in vitro approaches. We administered 4 mg/kg of warfarin sodium either orally or intravenously to juvenile green sea turtles. The turtles exhibited slow pharmacokinetics, and prolongation of prothrombin time (PT) was observed only with intravenous warfarin administration. We also conducted an in vitro investigation using liver microsomes from green sea turtles, and two other turtle species (softshell turtle and red-eared slider) and rats. The cytochrome P450 metabolic activity in the liver of green sea turtles was lower than in rats. Additionally, vitamin K epoxide reductase (VKOR), which is the target enzyme of ARs, was inhibited by warfarin in the turtles at lower concentration levels than in rats. These data indicate that turtles may be more sensitive to ARs than rats. We expect that these findings will be helpful for sea turtle conservation following accidental AR-broadcast incidents.

Convergent Losses of TLR5 Suggest Altered Extracellular Flagellin Detection in Four Mammalian Lineages

Toll-like receptors (TLRs) play an important role for the innate immune system by detecting pathogen-associated molecular patterns. TLR5 encodes the major extracellular receptor for bacterial flagellin and frequently evolves under positive selection, consistent with coevolutionary arms races between the host and pathogens. Furthermore, TLR5 is inactivated in several vertebrates and a TLR5 stop codon polymorphism is widespread in human populations. Here, we analyzed the genomes of 120 mammals and discovered that TLR5 is convergently lost in four independent lineages, comprising guinea pigs, Yangtze river dolphin, pinnipeds, and pangolins. Validated inactivating mutations, absence of protein-coding transcript expression, and relaxed selection on the TLR5 remnants confirm these losses. PCR analysis further confirmed the loss of TLR5 in the pinniped stem lineage. Finally, we show that TLR11, encoding a second extracellular flagellin receptor, is also absent in these four lineages. Independent losses of TLR5 and TLR11 suggest that a major pathway for detecting flagellated bacteria is not essential for different mammals and predicts an impaired capacity to sense extracellular flagellin.

A cross-sectional study of owner-reported health in Canadian and American cats fed meat- and plant-based diets

Background

Cats, being obligate carnivores, have unique dietary requirements for nutrients most commonly found in dietary ingredients of animal origin. As such, feeding a diet devoid of animal-derived ingredients has been postulated as a possible cause of nutrient imbalances and adverse health outcomes. A small proportion of cat owners feed strictly plant-based diets to the cats in their care, yet the health and wellness of cats fed these diets has not been well documented.

Results

A total of 1325 questionnaires were complete enough for inclusion. The only exclusion criterion was failure to answer all questions. Most cats, 65% (667/1026), represented in the survey were fed a meat-based diet and 18.2% (187/1026) were fed a plant-based diet, with the rest fed either a combination of plant-based with meat-based (69/1026, 6.7%) or indeterminable (103/1026, 10%). Cat age ranged from 4 months to 23 years, with a median of 7 years, and was not associated with diet type. No differences in reported lifespan were detected between diet types. Fewer cats fed plant-based diets reported to have gastrointestinal and hepatic disorders. Cats fed plant-based diets were reported to have more ideal body condition scores than cats fed a meat-based diet. More owners of cats fed plant-based diets reported their cat to be in very good health.

Conclusions

Cat owner perception of the health and wellness of cats does not appear to be adversely affected by being fed a plant-based diet. Contrary to expectations, owners perceived no body system or disorder to be at particular risk when feeding a plant-based diet to cats. This study collected information from cat owners and is subject to bias, as well as methodological limitations. Further research is warranted to determine if these results are replicable in a prospective investigation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-021-02754-8.

Evolution of the mammalian insulin (Ins) gene; Changes in proteolytic processing

Disruption of insulin signaling in humans leads to diabetes yet changes in insulin function is tolerated in some species. Taking advantage of the large number of publicly available mammalian genome sequences I identified insulin gene (Ins) in the genomes of 151 of 156 mammalian species with well-annotated genomes, of which 141 had complete Ins coding sequences. Complete Ins coding sequences were identified from 8 additional species that lack complete genomes. Duplicated Ins genes were found in 12 rodents (9 with complete genomes) resulting in the identification of a total of 161 complete mammalian Ins coding sequences. While all 161 proinsulin protein sequences were predicted to have functional signal peptides, which should allow secretion of the hormone, unexpectedly, substitutions were found at prohormone convertase processing sites in sequences from 6 species, 2 from Chiroptera (Myotis brandtii and M. lucifugus) and 4 from Afrotheria (Chrysochloris asiatica, Echinops telfairi, Elephantulus edwardii, and Orycteropus afer). Both basic residues at the C-peptide-A-chain junction in the bats M. brandtii and M. lucifugas are replaced, which should prevent processing. Replacements of a single basic residue are found at the B-chain-C-peptide junction, in the two bats, and at the C-peptide-A-chain junction, in 4 species of Afrotheria, processing sites that suggest impaired processing. In addition, a large number of substitutions at sites that interact with the insulin receptor were found in the insulin sequences from M. brandtii and M. lucifugas suggesting a change in biological function.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

GBA3: a polymorphic pseudogene in humans that experienced repeated gene loss during mammalian evolution

The gene encoding the cytosolic β-glucosidase GBA3 shows pseudogenization due to a truncated allele (rs358231) that is polymorphic in humans. Since this enzyme is involved in the transformation of many plant β-glycosides, this particular case of gene loss may have been influenced by dietary adaptations during evolution. In humans, apart from the inactivating allele, we found that GBA3 accumulated additional damaging mutations, implying an extensive GBA3 loss. The allelic distribution of loss-of-function alleles revealed significant differences between human populations which can be partially related with their staple diet. The analysis of mammalian orthologs disclosed that GBA3 underwent at least nine pseudogenization events. Most events of pseudogenization occurred in carnivorous lineages, suggesting a possible link to a β-glycoside poor diet. However, GBA3 was also lost in omnivorous and herbivorous species, hinting that the physiological role of GBA3 is not fully understood and other unknown causes may underlie GBA3 pseudogenization. Such possibility relies upon a putative role in sialic acid biology, where GBA3 participates in a cellular network involving NEU2 and CMAH. Overall, our data shows that the recurrent loss of GBA3 in mammals is likely to represent an evolutionary endpoint of the relaxation of selective constraints triggered by diet-related factors.

PseudoChecker: an integrated online platform for gene inactivation inference

The rapid expansion of high-quality genome assemblies, exemplified by ongoing initiatives such as the Genome-10K and i5k, demands novel automated methods to approach comparative genomics. Of these, the study of inactivating mutations in the coding region of genes, or pseudogenization, as a source of evolutionary novelty is mostly overlooked. Thus, to address such evolutionary/genomic events, a systematic, accurate and computationally automated approach is required. Here, we present PseudoChecker, the first integrated online platform for gene inactivation inference. Unlike the few existing methods, our comparative genomics-based approach displays full automation, a built-in graphical user interface and a novel index, PseudoIndex, for an empirical evaluation of the gene coding status. As a multi-platform online service, PseudoChecker simplifies access and usability, allowing a fast identification of disruptive mutations. An analysis of 30 genes previously reported to be eroded in mammals, and 30 viable genes from the same lineages, demonstrated that PseudoChecker was able to correctly infer 97% of loss events and 95% of functional genes, confirming its reliability. PseudoChecker is freely available, without login required, at http://pseudochecker.ciimar.up.pt.

First assessment of pollutant exposure in two balaenopterid whale populations sampled in the Svalbard Archipelago, Norway

Pollutant concentrations are poorly known for the largest animals on Earth, blue whales Balaenoptera musculus and fin whales Balaenoptera physalus. In this study, concentrations of persistent organic pollutants (POPs) were determined in blubber biopsies and stable isotope values for nitrogen (δ¹⁵N) and carbon (δ¹³C) were measured using skin biopsies for 18 blue whales and 12 fin whales sampled in waters surrounding the Svalbard Archipelago, Norway. The samples were collected in summer during the period 2014-2018. POPs were dominated by DDTs, PCBs and toxaphenes, with median concentrations in blue/fin whales being 208/341, 127/275 and 133/233 ng/g lipid weight, respectively. Linear models indicated that pollutant concentrations were 1.6-3 times higher in fin whales than in blue whales, which is likely related to the higher trophic positions of fin whales, as indicated by their higher δ¹⁵N. Lower δ¹³C in fin whales suggests that they feed at higher latitudes than blue whales; these values were not correlated with pollutant concentrations. Pollutant levels were approximately twice as high in males compared to females (intraspecifically), which indicates that females of these species offload pollutants to their offspring during gestation and lactation, similar to many other mammalian species. Pollutant concentrations in balaenopterid whales from Svalbard waters were generally much lower than in conspecific whales from the Mediterranean Sea or the Gulf of California, but higher than those in conspecifics from the Antarctic Peninsula.

Novel phylogenetic methods are needed for understanding gene function in the era of mega-scale genome sequencing

Ongoing large-scale genome sequencing projects are forecasting a data deluge that will almost certainly overwhelm current analytical capabilities of evolutionary genomics. In contrast to population genomics, there are no standardized methods in evolutionary genomics for extracting evolutionary and functional (e.g. gene-trait association) signal from genomic data. Here, we examine how current practices of multi-species comparative genomics perform in this aspect and point out that many genomic datasets are under-utilized due to the lack of powerful methodologies. As a result, many current analyses emphasize gene families for which some functional data is already available, resulting in a growing gap between functionally well-characterized genes/organisms and the universe of unknowns. This leaves unknown genes on the 'dark side' of genomes, a problem that will not be mitigated by sequencing more and more genomes, unless we develop tools to infer functional hypotheses for unknown genes in a systematic manner. We provide an inventory of recently developed methods capable of predicting gene-gene and gene-trait associations based on comparative data, then argue that realizing the full potential of whole genome datasets requires the integration of phylogenetic comparative methods into genomics, a rich but underutilized toolbox for looking into the past.

CRISPR-Cas9-Mutated Pregnane X Receptor (pxr) Retains Pregnenolone-induced Expression of cyp3a65 in Zebrafish (Danio rerio) Larvae

Pregnane X receptor (PXR; NR1I2) is a nuclear receptor that regulates transcriptional responses to drug or xenobiotic exposure, including induction of CYP3A transcription, in many vertebrate species. PXR is activated by a wide range of ligands that differ across species, making functional studies on its role in the chemical defensome most relevant when approached in a species-specific manner. Knockout studies in mammals have shown a requirement for PXR in ligand-dependent activation of CYP3A expression or reporter gene activity. Morpholino knockdown of Pxr in zebrafish indicated a similar requirement. Here, we report on the generation of 2 zebrafish lines each carrying a heritable deletion in the pxr coding region, predicted to result in loss of a functional gene product. To our surprise, larvae homozygous for either of the pxr mutant alleles retain their ability to induce cyp3a65 mRNA expression following exposure to the established zebrafish Pxr ligand, pregnenolone. Thus, zebrafish carrying pxr alleles with deletions in either the DNA binding or the ligand-binding domains did not yield a loss-of-function phenotype, suggesting that a compensatory mechanism is responsible for cyp3a65 induction. Alternative possibilities are that Pxr is not required for the induction of selected genes, or that truncated yet functional mutant Pxr is sufficient for the downstream transcriptional effects. It is crucial that we develop a better understanding for the role of Pxr in this important biomedical test species. This study highlights the potential for compensatory mechanisms to avoid deleterious effects arising from gene mutations.

A genome alignment of 120 mammals highlights ultraconserved element variability and placenta-associated enhancers

BACKGROUND

Multiple alignments of mammalian genomes have been the basis of many comparative genomic studies aiming at annotating genes, detecting regions under evolutionary constraint, and studying genome evolution. A key factor that affects the power of comparative analyses is the number of species included in a genome alignment.

RESULTS

To utilize the increased number of sequenced genomes and to provide an accessible resource for genomic studies, we generated a mammalian genome alignment comprising 120 species. We used this alignment and the CESAR method to provide protein-coding gene annotations for 119 non-human mammals. Furthermore, we illustrate the utility of this alignment by 2 exemplary analyses. First, we quantified how variable ultraconserved elements (UCEs) are among placental mammals. Leveraging the high taxonomic coverage in our alignment, we estimate that UCEs contain on average 4.7%-15.6% variable alignment columns. Furthermore, we show that the center regions of UCEs are generally most constrained. Second, we identified enhancer sequences that are only conserved in placental mammals. We found that these enhancers are significantly associated with placenta-related genes, suggesting that some of these enhancers may be involved in the evolution of placental mammal-specific aspects of the placenta.

CONCLUSION

The 120-mammal alignment and all other data are available for analysis and visualization in a genome browser at https://genome-public.pks.mpg.de/and for download at https://bds.mpi-cbg.de/hillerlab/120MammalAlignment/.