Phylogenomic analysis resolves the interordinal relationships and rapid diversification of the laurasiatherian mammals

Phylogenomics of Afrotherian mammals and improved resolution of extant Paenungulata

Molecular investigations have gathered a diverse set of mammals-predominantly African natives like elephants, hyraxes, and aardvarks-into a clade known as Afrotheria. Nevertheless, the precise phylogenetic relationships among these species remain contentious. Here, we sourced orthologous markers and ultraconserved elements to discern the interordinal connections among Afrotherian mammals. Our phylogenetic analyses bolster the common origin of Afroinsectiphilia and Paenungulata, and propose Afrosoricida as the closer relative to Macroscelidea rather than Tubulidentata, while also challenging the notion of Sirenia and Hyracoidea as sister taxa. The approximately unbiased test and the gene concordance factor uniformly recognized the alliance of Proboscidea with Hyracoidea as the dominant topology within Paenungulata. Investigation into sites with extremly high phylogenetic signal unveiled their potential to intensify conflicts in the Paenungulata topology. Subsequent exploration suggested that incomplete lineage sorting was predominantly responsible for the observed contentious relationships, whereas introgression exerted a subsidiary influence. The divergence times estimated in our study hint at the Cretaceous-Paleogene (K-Pg) extinction event as a catalyst for Afrotherian diversification. Overall, our findings deliver a tentative but insightful overview of Afrotheria phylogeny and divergence, elucidating these relationships through the lens of phylogenomics.

Evolution and conservation genetics of pangolins

Pangolins (Pholidota, Manidae) are classified as an evolutionarily distinct and globally endangered mammal due to their unique morphology (nail-like scales and a myrmecophagous diet) and being the victim of heavy poaching and worldwide trafficking. As such, pangolins serve as a textbook example for studying the special phenotypic evolutionary adaptations and conservation genetics of an endangered species. Recent years have demonstrated significant advancements in the fields of molecular genetics and genomics, which have translated to a series of important research achievements and breakthroughs concerning the evolution and conservation genetics of pangolins. This review comprehensively presents the hitherto advances in phylogeny, adaptive evolution, conservation genetics, and conservation genomics that are related to pangolins, which will provide an ample understanding of their diversity, molecular adaptation mechanisms, and evolutionary potentials. In addition, we highlight the priority of investigating species/population diversity among pangolins and suggest several avenues of research that are highly relevant for future pangolin conservation.

Comparative Genomics Reveal Phylogenetic Relationship and Chromosomal Evolutionary Events of Eight Cervidae Species

Cervidae represents a family that is not only rich in species diversity but also exhibits a wide range of karyotypes. The controversies regarding the phylogeny and classification of Cervidae still persist. The flourishing development of the genomic era has made it possible to address these issues at the genomic level. Here, the genomes of nine species were used to explore the phylogeny and chromosomal evolutionary events of Cervidae. By conducting whole-genome comparisons, we identified single-copy orthologous genes across the nine species and constructed a phylogenetic tree based on the single-copy orthologous genes sequences, providing new insights into the phylogeny of Cervidae, particularly the phylogenetic relationship among sika deer, red deer, wapiti and Tarim red deer. Gene family analysis revealed contractions in the olfactory receptor gene family and expansions in the histone gene family across eight Cervidae species. Furthermore, synteny analysis was used to explore the chromosomal evolutionary events of Cervidae species, revealing six chromosomal fissions during the evolutionary process from Bovidae to Cervidae. Notably, specific chromosomal fusion events were found in four species of Cervus, and a unique chromosomal fusion event was identified in Muntiacus reevesi. Our study further completed the phylogenetic relationship within the Cervidae and demonstrated the feasibility of inferring species phylogeny at the whole-genome level. Additionally, our findings on gene family evolution and the chromosomal evolutionary events in eight Cervidae species lay a foundation for comprehensive research of the evolution of Cervidae.

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Improved mammalian family phylogeny using gap-rare multiple sequence alignment: A timetree of extant placentals and marsupials

The timing of mammalian diversification in relation to the Cretaceous-Paleogene (KPg) mass extinction continues to be a subject of substantial debate. Previous studies have either focused on limited taxonomic samples with available whole-genome data or relied on short sequence alignments coupled with extensive species samples. In the present study, we improved an existing dataset from the landmark study of Meredith et al. (2011) by filling in missing fragments and further generated another dataset containing 120 taxa and 98 exonic markers. Using these two datasets, we then constructed phylogenies for extant mammalian families, providing improved resolution of many conflicting relationships. Moreover, the timetrees generated, which were calibrated using appropriate molecular clock models and multiple fossil records, indicated that the interordinal diversification of placental mammals initiated before the Late Cretaceous period. Additionally, intraordinal diversification of both extant placental and marsupial lineages accelerated after the KPg boundary, supporting the hypothesis that the availability of numerous vacant ecological niches subsequent to the mass extinction event facilitated rapid diversification. Thus, our results support a scenario of placental radiation characterized by both basal cladogenesis and active interordinal divergences spanning from the Late Cretaceous into the Paleogene.

Pangolin Genomes Offer Key Insights and Resources for the World's Most Trafficked Wild Mammals

Pangolins form a group of scaly mammals that are trafficked at record numbers for their meat and purported medicinal properties. Despite their conservation concern, knowledge of their evolution is limited by a paucity of genomic data. We aim to produce exhaustive genomic resources that include 3,238 orthologous genes and whole-genome polymorphisms to assess the evolution of all eight extant pangolin species. Robust orthologous gene-based phylogenies recovered the monophyly of the three genera and highlighted the existence of an undescribed species closely related to Southeast Asian pangolins. Signatures of middle Miocene admixture between an extinct, possibly European, lineage and the ancestor of Southeast Asian pangolins, provide new insights into the early evolutionary history of the group. Demographic trajectories and genome-wide heterozygosity estimates revealed contrasts between continental versus island populations and species lineages, suggesting that conservation planning should consider intraspecific patterns. With the expected loss of genomic diversity from recent, extensive trafficking not yet realized in pangolins, we recommend that populations be genetically surveyed to anticipate any deleterious impact of the illegal trade. Finally, we produce a complete set of genomic resources that will be integral for future conservation management and forensic endeavors for pangolins, including tracing their illegal trade. These comprise the completion of whole-genomes for pangolins through the hybrid assembly of the first reference genome for the giant pangolin (Smutsia gigantea) and new draft genomes (∼43x-77x) for four additional species, as well as a database of orthologous genes with over 3.4 million polymorphic sites.

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.

Novel insights into Pinus species plastids genome through phylogenetic relationships and repeat sequence analysis

Pinus is one of the most economical and ecological important conifers, model specie for studying sequence divergence and molecular phylogeney of gymnosperms. The less availability of information for genome resources enable researchers to conduct evolutionary studies of Pinus species. To improve understanding, we firstly reported, previously released chloroplast genome of 72 Pinus species, the sequence variations, phylogenetic relationships and genome divergence among Pinus species. The results displayed 7 divergent hotspot regions (trnD-GUC, trnY-GUA, trnH-GUG, ycf1, trnL-CAA, trnK-UUU and trnV-GAC) in studied Pinus species, which holds potential to utilized as molecular genetic markers for future phylogenetic studies in Pinnus species. In addition, 3 types of repeats (tandem, palindromic and dispersed) were also studied in Pinus species under investigation. The outcome showed P. nelsonii had the highest, 76 numbers of repeats, while P. sabiniana had the lowest, 13 13 numbers of repeats. It was also observed, constructed phylogenetic tree displayed division into two significant diverged clades: single needle (soft pine) and double-needle (hard pine). Theoutcome of present investigation, based on the whole chloroplast genomes provided novel insights into the molecular based phylogeny of the genus Pinus which holds potential for its utilization in future studies focusing genetic diversity in Pinnus species.

Quantification of Urinary Sex Steroids in the Big Brown Bat (Eptesicus fuscus)

AbstractBats (order Chiroptera) are the second largest group of mammals, diverging ~52.5 million years ago. Many species exhibit an unusual reproductive cycle and extreme longevity without reproductive senescence, yet steroid profiles exist for few bats. Big brown bats (Eptesicus fuscus) are temperate insectivores found throughout North America. They mate promiscuously in fall, store sperm during winter hibernation, and have delayed ovulation and fertilization in spring. Here, we report the first urinary steroid profile in bats by quantifying 17β-estradiol (E₂) in captive male and female E. fuscus across their reproductive cycle. Male bats had higher urinary E₂ levels than females, and adults had higher levels than yearlings following creatinine adjustment for hydration. In nonpregnant females, several seasonal differences in creatinine-adjusted and unadjusted urinary E₂ levels were observed. Urinary E₂ was higher in males than females in winter for both conditions and in autumn for creatinine-adjusted levels. We quantified progesterone (P₄) in a subset of females. In nonpregnant females, urinary P₄ was constant across seasons except for unadjusted levels, which were highest in the summer. In pregnant females, urinary E₂ and P₄ levels peaked beginning ~20 d before parturition, with both steroids returning to baseline in the following weeks. Knowing how urinary steroid levels fluctuate with age and sex and across the annual season is key to understanding reproductive cycling in bats. Our research furthers the potential for bats as a model for medical reproductive research. Moreover, it complements previous studies on the potential role of steroids in primer pheromonal effects in bats.

Molecular Evolution of Tooth-Related Genes Provides New Insights into Dietary Adaptations of Mammals

Mammals have evolved different tooth phenotypes that are hypothesized to be associated with feeding habits. However, the genetic basis for the linkage has not been well explored. In this study, we investigated 13 tooth-related genes, including seven enamel-related genes (AMELX, AMBN, ENAM, AMTN, ODAM, KLK4 and MMP20) and six dentin-related genes (DSPP, COL1A1, DMP1, IBSP, MEPE and SPP1), from 63 mammals to determine their evolutionary history. Our results showed that different evolutionary histories have evolved among divergent feeding habits in mammals. There was stronger positive selection for eight genes (ENAM, AMTN, ODAM, KLK4, DSPP, DMP1, COL1A1, MEPE) in herbivore lineages. In addition, AMELX, AMBN, ENAM, AMTN, MMP20 and COL1A1 underwent accelerated evolution in herbivores. While relatively strong positive selection was detected in IBSP, SPP1, and DSPP, accelerated evolution was only detected for MEPE and SPP1 genes among the carnivorous lineages. We found positive selection on AMBN and ENAM genes for omnivorous primates in the catarrhini clade. Interestingly, a significantly positive association between the evolutionary rate of ENAM, ODAM, KLK4, MMP20 and the average enamel thickness was found in primates. Additionally, we found molecular convergence in some amino acid sites of tooth-related genes among the lineages whose feeding habit are similar. The positive selection of related genes might promote the formation and bio-mineralization of tooth enamel and dentin, which would make the tooth structure stronger. Our results revealed that mammalian tooth-related genes have experienced variable evolutionary histories, which provide some new insights into the molecular basis of dietary adaptation in mammals.

The earliest Asian bats (Mammalia: Chiroptera) address major gaps in bat evolution

Bats dispersed widely after evolving the capacity for powered flight, and fossil bats are known from the early Eocene of most continents. Until now, however, bats have been conspicuously absent from the early Eocene of mainland Asia. Here, we report two teeth from the Junggar Basin of northern Xinjiang, China belonging to the first known early Eocene bats from Asia, representing arguably the most plesiomorphic bat molars currently recognized. These teeth combine certain bat synapomorphies with primitive traits found in other placental mammals, thereby potentially illuminating dental evolution among stem bats. The Junggar Basin teeth suggest that the dentition of the stem chiropteran family Onychonycteridae is surprisingly derived, although their postcranial anatomy is more primitive than that of any other Eocene bats. Additional comparisons with stem bat families Icaronycteridae and Archaeonycteridae fail to identify unambiguous synapomorphies for the latter taxa, raising the possibility that neither is monophyletic as currently recognized. The presence of highly plesiomorphic bats in the early Eocene of central Asia suggests that this region was an important locus for the earliest, transitional phases of bat evolution, as has been demonstrated for other placental mammal orders including Lagomorpha and Rodentia.

ACPT gene is inactivated in mammalian lineages that lack enamel or teeth

Loss of tooth or enamel is widespread in multiple mammal lineages. Although several studies have been reported, the evolutionary mechanisms of tooth/enamel loss are still unclear. Most previous studies have found that some tooth-related genes have been inactivated in toothless and/or enamel-less mammals, such as ENAM, ODAM, C4orf26, AMBN, AMTN, DSPP, etc. Here, we conducted evolutionary analyses on ACPT playing a key role in amelogenesis, to interrogate the mechanisms. We obtained the ACPT sequences from 116 species, including edentulous and enamel-less mammals. The results shows that variant ORF-disrupting mutations were detected in ACPT coding region among nine edentulous baleen whales and three enamel-less taxa (pygmy sperm whale, aardvark, nine-banded armadillo). Furtherly, selective pressure uncovered that the selective constraints have been relaxed among all toothless and enamel-less lineages. Moreover, our results support the hypothesis that mineralized teeth were lost or degenerated in the common ancestor of crown Mysticeti through two shared single-base sites deletion in exon 4 and 5 of ACPT among all living baleen whales. D_N/d_S values on transitional branches were used to estimate ACPT inactivation records. In the case of aardvark, inactivation of ACPT was estimated at ~23.60–28.32 Ma, which is earlier than oldest aardvark fossil record (Orycteropus minutus, ~19 Ma), suggesting that ACPT inactivation may result in degeneration or loss of enamel. Conversely, the inactivation time of ACPT estimated in armadillo (~10.18–11.30 Ma) is later than oldest fossil record, suggesting that inactivation of ACPT may result from degeneration or loss of enamel in these mammals. Our findings suggested that different mechanisms of degeneration of tooth/enamel might exist among toothless and enamel-less lineages during evolution. Our study further considered that ACPT is a novel gene for studying tooth evolution.

Diverse phylogenomic datasets uncover a concordant scenario of laurasiatherian interordinal relationships

Resolving the interordinal relationships in the mammalian superorder Laurasiatheria has been among the most intractable problems in higher-level mammalian systematics, with many conflicting hypotheses having been proposed. The present study collected three different sources of genome-scale data with comprehensive taxon sampling of laurasiatherian species, including two protein-coding datasets (4,186 protein-coding genes for an amino acid dataset comprising 2,761,247 amino acid residues and a nucleotide dataset comprising 5,516,340 nucleotides from 1st and 2nd codon positions), an intronic dataset (1,210 introns comprising 1,162,723 nucleotides) and an ultraconserved elements (UCEs) dataset (1,246 UCEs comprising 1,946,472 nucleotides) from 40 species representing all six laurasiatherian orders and 7 non-laurasiatherian outgroups. Remarkably, phylogenetic trees reconstructed with the four datasets using different tree-building methods (RAxML, FastTree, ASTRAL and MP-EST) all supported the relationship (Eulipotyphla, (Chiroptera, ((Carnivora, Pholidota), (Cetartiodactyla, Perissodactyla)))). We find a resolution of interordinal relationships of Laurasiatheria among all types of markers used in the present study, and the likelihood ratio tests for tree comparisons confirmed that the present tree topology is the optimal hypothesis compared to other examined hypotheses. Jackknifing subsampling analyses demonstrate that the results of laurasiatherian tree reconstruction varied with the number of loci and ordinal representatives used, which are likely the two main contributors to phylogenetic disagreements of Laurasiatheria seen in previous studies. Our study provides significant insight into laurasiatherian evolution, and moreover, an important methodological strategy and reference for resolving phylogenies of adaptive radiation, which have been a long-standing challenge in the field of phylogenetics.

Comparative genomic analyses illuminate the distinct evolution of megabats within Chiroptera

The revision of the sub-order Microchiroptera is one of the most intriguing outcomes in recent mammalian molecular phylogeny. The unexpected sister–taxon relationship between rhinolophoid microbats and megabats, with the exclusion of other microbats, suggests that megabats arose in a relatively short period of time from a microbat-like ancestor. In order to understand the genetic mechanism underlying adaptive evolution in megabats, we determined the whole-genome sequences of two rousette megabats, Leschenault’s rousette (Rousettus leschenaultia) and the Egyptian fruit bat (R. aegyptiacus). The sequences were compared with those of 22 other mammals, including nine bats, available in the database. We identified that megabat genomes are distinct in that they have extremely low activity of SINE retrotranspositions, expansion of two chemosensory gene families, including the trace amine receptor (TAAR) and olfactory receptor (OR), and elevation of the dN/dS ratio in genes for immunity and protein catabolism. The adaptive signatures discovered in the genomes of megabats may provide crucial insight into their distinct evolution, including key processes such as virus resistance, loss of echolocation, and frugivorous feeding.

Insights into Dietary Switch in Cetaceans: Evidence from Molecular Evolution of Proteinases and Lipases

Fossil evidence suggests that cetaceans evolved from artiodactylans. Thus, there was a major dietary change from herbivorous to carnivorous during their transition from a terrestrial to an aquatic environment. However, the molecular evolutionary mechanisms underlying this dietary switch have not been well investigated. Evidence of positive selection of digestive proteinases and lipases of cetaceans was detected: (1) For the four pancreatic proteinase families (carboxypeptidase, trypsin, chymotrypsin, and elastase) examined in this study, each family included only a single intact gene (e.g., CPA1, PRSS1, CTRC, and CELA3B) that had no ORF-disrupted or premature stop codons, whereas other members of each family had become pseudogenized. Further selective pressure analysis showed that three genes (PRSS1, CTRC, and CELA3B) were subjected to significant positive selection in cetaceans. (2) For digestive proteinases from the stomach, PGA was identified to be under positive selection. (3) Intense positive selection was also detected for the lipase gene PLRP2 in cetaceans. In addition, parallel /convergent amino acid substitutions between cetaceans and carnivores, two groups of mammals that have evolved similar feeding habits, were identified in 10 of the 12 functional genes. Although pseudogenization resulted in each family of pancreatic proteinases only retaining one intact gene copy in cetacean genomes, positive selection might have driven pancreatic proteinases, stomach proteinases, and lipases to adaptively evolve a stronger ability to digest a relatively higher proportion of proteins and lipids from animal foods. This study can provide some novel insights into the molecular mechanism of cetacean dietary changes during their transition from land to sea.

The Prevalence and Impact of Model Violations in Phylogenetic Analysis

In phylogenetic inference, we commonly use models of substitution which assume that sequence evolution is stationary, reversible, and homogeneous (SRH). Although the use of such models is often criticized, the extent of SRH violations and their effects on phylogenetic inference of tree topologies and edge lengths are not well understood. Here, we introduce and apply the maximal matched-pairs tests of homogeneity to assess the scale and impact of SRH model violations on 3,572 partitions from 35 published phylogenetic data sets. We show that roughly one-quarter of all the partitions we analyzed (23.5%) reject the SRH assumptions, and that for 25% of data sets, tree topologies inferred from all partitions differ significantly from topologies inferred using the subset of partitions that do not reject the SRH assumptions. This proportion increases when comparing trees inferred using the subset of partitions that rejects the SRH assumptions, to those inferred from partitions that do not reject the SRH assumptions. These results suggest that the extent and effects of model violation in phylogenetics may be substantial. They highlight the importance of testing for model violations and possibly excluding partitions that violate models prior to tree reconstruction. Our results also suggest that further effort in developing models that do not require SRH assumptions could lead to large improvements in the accuracy of phylogenomic inference. The scripts necessary to perform the analysis are available in https://github.com/roblanf/SRHtests, and the new tests we describe are available as a new option in IQ-TREE (http://www.iqtree.org).

Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae)

BACKGROUND

Paris (Melanthiaceae) is an economically important but taxonomically difficult genus, which is unique in angiosperms because some species have extremely large nuclear genomes. Phylogenetic relationships within Paris have long been controversial. Based on complete plastomes and nuclear ribosomal DNA (nrDNA) sequences, this study aims to reconstruct a robust phylogenetic tree and explore historical biogeography and clade diversification in the genus.

RESULTS

All 29 species currently recognized in Paris were sampled. Whole plastomes and nrDNA sequences were generated by the genome skimming approach. Phylogenetic relationships were reconstructed using the maximum likelihood and Bayesian inference methods. Based on the phylogenetic framework and molecular dating, biogeographic scenarios and historical diversification of Paris were explored. Significant conflicts between plastid and nuclear datasets were identified, and the plastome tree is highly congruent with past interpretations of the morphology. Ancestral area reconstruction indicated that Paris may have originated in northeastern Asia and northern China, and has experienced multiple dispersal and vicariance events during its diversification. The rate of clade diversification has sharply accelerated since the Miocene/Pliocene boundary.

CONCLUSIONS

Our results provide important insights for clarifying some of the long-standing taxonomic debates in Paris. Cytonuclear discordance may have been caused by ancient and recent hybridizations in the genus. The climatic and geological changes since the late Miocene, such as the intensification of Asian monsoon and the rapid uplift of Qinghai-Tibet Plateau, as well as the climatic fluctuations during the Pleistocene, played essential roles in driving range expansion and radiative diversification in Paris. Our findings challenge the theoretical prediction that large genome sizes may limit speciation.

Divergent Evolution of TRC Genes in Mammalian Niche Adaptation

Mammals inhabit a wide variety of ecological niches, which in turn can be affected by various ecological factors, especially in relation to immunity. The canonical TRC repertoire (TRAC, TRBC, TRGC, and TRDC) codes C regions of T cell receptor chains that form the primary antigen receptors involved in the activation of cellular immunity. At present, little is known about the correlation between the evolution of mammalian TRC genes and ecological factors. In this study, four types canonical of TRC genes were identified from 37 mammalian species. Phylogenetic comparative methods (phyANOVA and PGLS) and selective pressure analyses among different groups of ecological factors (habitat, diet, and sociality) were carried out. The results showed that habitat was the major ecological factor shaping mammalian TRC repertoires. Specifically, trade-off between TRGC numbers and positive selection of TRAC and the balanced evolutionary rates between TRAC and TRDC genes were speculated as two main mechanisms in adaption to habitat and sociality. Overall, our study suggested divergent mechanisms for the evolution of TRCs, prompting mammalian immunity adaptions within diverse niches.

Elongation of the CA1 field of the septal hippocampus in ungulates

It is widely assumed that the hippocampal formation seen in laboratory rodents and in primates is typical of that seen in other mammals. We have tested this assumption by examining sections of brains of 56 mammals from 20 mammalian orders from images on the brainmuseum.org website. We found wide variation in the form of the hippocampal formation, the most extreme examples of which are seen in ungulates, which possess an unusual elongation of the distal CA1 of the septal hippocampus. This phenomenon has not previously been reported. In individual coronal sections of the brains of seven artiodactyl ungulates, the pyramidal layer of CA1 is four times as long as CA2 + CA3. In a perissodactyl ungulate (Burchell's zebra) the distal end of CA1 is so large that it forms a number of folds. A similar but less pronounced CA1 elongation was seen in the brains of 14 carnivores. A modest elongation of CA1 is also present in some other placental mammals, notably the elephant shrew, hyrax, capybara, beaver, and rabbit. The elongation was not present in brains of primates, marsupials, or monotremes. The distal part of CA1 has been shown to play a role in object integration into the spatial map. We hypothesize that the distal CA1 enlargement could serve to enhance the ability to integrate objects into spatial navigation, which would be an advantage for migrating herds of ungulates. We suggest that the remarkable elongation of Q5 CA1 represents a major evolutionary specialization in the ungulates.

Phylogenetic and mutational analyses of human LEUTX, a homeobox gene implicated in embryogenesis

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

Differences in Performance among Test Statistics for Assessing Phylogenomic Model Adequacy

Statistical phylogenetic analyses of genomic data depend on models of nucleotide or amino acid substitution. The adequacy of these substitution models can be assessed using a number of test statistics, allowing the model to be rejected when it is found to provide a poor description of the evolutionary process. A potentially valuable use of model-adequacy test statistics is to identify when data sets are likely to produce unreliable phylogenetic estimates, but their differences in performance are rarely explored. We performed a comprehensive simulation study to identify test statistics that are sensitive to some of the most commonly cited sources of phylogenetic estimation error. Our results show that, for many test statistics, traditional thresholds for assessing model adequacy can fail to reject the model when the phylogenetic inferences are inaccurate and imprecise. This is particularly problematic when analysing loci that have few informative sites. We propose new thresholds for assessing substitution model adequacy and demonstrate their effectiveness in analyses of three phylogenomic data sets. These thresholds lead to frequent rejection of the model for loci that yield topological inferences that are imprecise and are likely to be inaccurate. We also propose the use of a summary statistic that provides a practical assessment of overall model adequacy. Our approach offers a promising means of enhancing model choice in genome-scale data sets, potentially leading to improvements in the reliability of phylogenomic inference.

Using Phylogenomic Data to Explore the Effects of Relaxed Clocks and Calibration Strategies on Divergence Time Estimation: Primates as a Test Case

Primates have long been a test case for the development of phylogenetic methods for divergence time estimation. Despite a large number of studies, however, the timing of origination of crown Primates relative to the Cretaceous-Paleogene (K-Pg) boundary and the timing of diversification of the main crown groups remain controversial. Here, we analysed a data set of 372 taxa (367 Primates and 5 outgroups, 3.4 million aligned base pairs) that includes nine primate genomes. We systematically explore the effect of different interpretations of fossil calibrations and molecular clock models on primate divergence time estimates. We find that even small differences in the construction of fossil calibrations can have a noticeable impact on estimated divergence times, especially for the oldest nodes in the tree. Notably, choice of molecular rate model (autocorrelated or independently distributed rates) has an especially strong effect on estimated times, with the independent rates model producing considerably more ancient age estimates for the deeper nodes in the phylogeny. We implement thermodynamic integration, combined with Gaussian quadrature, in the program MCMCTree, and use it to calculate Bayes factors for clock models. Bayesian model selection indicates that the autocorrelated rates model fits the primate data substantially better, and we conclude that time estimates under this model should be preferred. We show that for eight core nodes in the phylogeny, uncertainty in time estimates is close to the theoretical limit imposed by fossil uncertainties. Thus, these estimates are unlikely to be improved by collecting additional molecular sequence data. All analyses place the origin of Primates close to the K-Pg boundary, either in the Cretaceous or straddling the boundary into the Palaeogene.

Genomic organization and adaptive evolution of IGHC genes in marine mammals

•

The genomic organization of IGHC genes in marine mammal is similar to that of terrestrial relatives.

•

The number of the immunoglobulin heavy chain constant region genes vary among different mammals.

•

Different levels of selective pressures were detected between marine and terrestrial mammalian lineages.

Immunoglobulins are important elements of the adaptive immune system that bind to an immense variety of microbial antigens to neutralize infectivity and specify effector functions. In the present study, the immunoglobulin heavy chain constant region (IGHC) genes from marine mammals were identified and compared with those of their terrestrial relatives to explore their genomic organization and evolutionary characteristics. The genomic organization of marine mammal IGHC genes was shown to be conservative with other eutherian mammals. Stronger signals of positive selection on IGHC were revealed in terrestrial mammals than that in marine mammals with the branch-site model, displaying different selective pressure, which might suggest their divergent adaptations to contrasted environments.

Phylogenomic Resolution of the Phylogeny of Laurasiatherian Mammals: Exploring Phylogenetic Signals within Coding and Noncoding Sequences

The interordinal relationships of Laurasiatherian mammals are currently one of the most controversial questions in mammalian phylogenetics. Previous studies mainly relied on coding sequences (CDS) and seldom used noncoding sequences. Here, by data mining public genome data, we compiled an intron data set of 3,638 genes (all introns from a protein-coding gene are considered as a gene) (19,055,073 bp) and a CDS data set of 10,259 genes (20,994,285 bp), covering all major lineages of Laurasiatheria (except Pholidota). We found that the intron data contained stronger and more congruent phylogenetic signals than the CDS data. In agreement with this observation, concatenation and species-tree analyses of the intron data set yielded well-resolved and identical phylogenies, whereas the CDS data set produced weakly supported and incongruent results. Further analyses showed that the phylogeny inferred from the intron data is highly robust to data subsampling and change in outgroup, but the CDS data produced unstable results under the same conditions. Interestingly, gene tree statistical results showed that the most frequently observed gene tree topologies for the CDS and intron data are identical, suggesting that the major phylogenetic signal within the CDS data is actually congruent with that within the intron data. Our final result of Laurasiatheria phylogeny is (Eulipotyphla,((Chiroptera, Perissodactyla),(Carnivora, Cetartiodactyla))), favoring a close relationship between Chiroptera and Perissodactyla. Our study 1) provides a well-supported phylogenetic framework for Laurasiatheria, representing a step towards ending the long-standing "hard" polytomy and 2) argues that intron within genome data is a promising data resource for resolving rapid radiation events across the tree of life.

Evolutionary changes of Hox genes and relevant regulatory factors provide novel insights into mammalian morphological modifications

The diversity of body plans of mammals accelerates the innovation of lifestyles and the extensive adaptation to different habitats, including terrestrial, aerial and aquatic habitats. However, the genetic basis of those phenotypic modifications, which have occurred during mammalian evolution, remains poorly explored. In the present study, we synthetically surveyed the evolutionary pattern of Hox clusters that played a powerful role in the morphogenesis along the head–tail axis of animal embryos and the main regulatory factors (Mll, Bmi1 and E2f6) that control the expression of Hox genes. A deflected density of repetitive elements and lineage‐specific radical mutations of Mll have been determined in marine mammals with morphological changes, suggesting that evolutionary changes may alter Hox gene expression in these lineages, leading to the morphological modification of these lineages. Although no positive selection was detected at certain ancestor nodes of lineages, the increased ω values of Hox genes implied the relaxation of functional constraints of these genes during the mammalian evolutionary process. More importantly, 49 positively‐selected sites were identified in mammalian lineages with phenotypic modifications, indicating adaptive evolution acting on Hox genes and regulatory factors. In addition, 3 parallel amino acid substitutions in some Hox genes were examined in marine mammals, which might be responsible for their streamlined body.

Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements

The phylogeny of eutherian mammals contains some of the most recalcitrant nodes in the tetrapod tree of life. We combined comprehensive taxon and character sampling to explore three of the most debated interordinal relationships among placental mammals. We performed in silico extraction of ultraconserved element loci from 72 published genomes and invitro enrichment and sequencing of ultraconserved elements from 28 additional mammals, resulting in alignments of 3,787 loci. We analyzed these data using concatenated and multispecies coalescent phylogenetic approaches, topological tests, and exploration of support among individual loci to identify the root of Eutheria and the sister groups of tree shrews (Scandentia) and horses (Perissodactyla). Individual loci provided weak, but often consistent support for topological hypotheses. Although many gene trees lacked accepted species-tree relationships, summary coalescent topologies were largely consistent with inferences from concatenation. At the root of Eutheria, we identified consistent support for a sister relationship between Xenarthra and Afrotheria (i.e., Atlantogenata). At the other nodes of interest, support was less consistent. We suggest Scandentia is the sister of Primatomorpha (Euarchonta), but we failed to reject a sister relationship between Scandentia and Glires. Similarly, we suggest Perissodactyla is sister to Cetartiodactyla (Euungulata), but a sister relationship between Perissodactyla and Chiroptera remains plausible.

Investigating Difficult Nodes in the Placental Mammal Tree with Expanded Taxon Sampling and Thousands of Ultraconserved Elements

The phylogeny of eutherian mammals contains some of the most recalcitrant nodes in the tetrapod tree of life. We combined comprehensive taxon and character sampling to explore three of the most debated interordinal relationships among placental mammals. We performed in silico extraction of ultraconserved element loci from 72 published genomes and invitro enrichment and sequencing of ultraconserved elements from 28 additional mammals, resulting in alignments of 3,787 loci. We analyzed these data using concatenated and multispecies coalescent phylogenetic approaches, topological tests, and exploration of support among individual loci to identify the root of Eutheria and the sister groups of tree shrews (Scandentia) and horses (Perissodactyla). Individual loci provided weak, but often consistent support for topological hypotheses. Although many gene trees lacked accepted species-tree relationships, summary coalescent topologies were largely consistent with inferences from concatenation. At the root of Eutheria, we identified consistent support for a sister relationship between Xenarthra and Afrotheria (i.e., Atlantogenata). At the other nodes of interest, support was less consistent. We suggest Scandentia is the sister of Primatomorpha (Euarchonta), but we failed to reject a sister relationship between Scandentia and Glires. Similarly, we suggest Perissodactyla is sister to Cetartiodactyla (Euungulata), but a sister relationship between Perissodactyla and Chiroptera remains plausible.

Speciation network in Laurasiatheria: retrophylogenomic signals

Rapid species radiation due to adaptive changes or occupation of new ecospaces challenges our understanding of ancestral speciation and the relationships of modern species. At the molecular level, rapid radiation with successive speciations over short time periods-too short to fix polymorphic alleles-is described as incomplete lineage sorting. Incomplete lineage sorting leads to random fixation of genetic markers and hence, random signals of relationships in phylogenetic reconstructions. The situation is further complicated when you consider that the genome is a mosaic of ancestral and modern incompletely sorted sequence blocks that leads to reconstructed affiliations to one or the other relative, depending on the fixation of their shared ancestral polymorphic alleles. The laurasiatherian relationships among Chiroptera, Perissodactyla, Cetartiodactyla, and Carnivora present a prime example for such enigmatic affiliations. We performed whole-genome screenings for phylogenetically diagnostic retrotransposon insertions involving the representatives bat (Chiroptera), horse (Perissodactyla), cow (Cetartiodactyla), and dog (Carnivora), and extracted among 162,000 preselected cases 102 virtually homoplasy-free, phylogenetically informative retroelements to draw a complete picture of the highly complex evolutionary relations within Laurasiatheria. All possible evolutionary scenarios received considerable retrotransposon support, leaving us with a network of affiliations. However, the Cetartiodactyla-Carnivora relationship as well as the basal position of Chiroptera and an ancestral laurasiatherian hybridization process did exhibit some very clear, distinct signals. The significant accordance of retrotransposon presence/absence patterns and flanking nucleotide changes suggest an important influence of mosaic genome structures in the reconstruction of species histories.

Macrosystematics of eutherian mammals combining HTS data to expand taxon coverage

In the last few years high-throughput sequencing technologies have permitted significant advances in mammalian phylogenetic studies from a genomic perspective. However, these studies have been restricted to a sparse number of species with available reference genomes. Thus, several issues inside the eutherian mammals phylogeny remain unresolved. This may be due in part to limited taxon sampling, as taxonomic density is known to affect phylogenetic resolution. In this context, we present a protocol to increase taxon coverage using high-throughput sequencing data (RNA or DNA) generated for other biological studies and available in public databases. Following this procedure we addressed pending or controversial issues concerning the phylogenetic position of Dermoptera, Pholidota and Chiroptera, considering multiple and independent loci. Also for Chiroptera and Arctoidea we evaluated the relationships of the lineages that compose it. Although the maximum number of genes used is moderate (95), in some cases taxon coverage doubles that of previous related studies. Globally, all coalescent-based (STAR, MP-EST and ASTRAL) and concatenated (IQ-TREE and BEAST2) methods used for species tree reconstruction were consistent to each other and most of interrogated nodes received high statistical support.

Diversification of Rosaceae since the Late Cretaceous based on plastid phylogenomics

Phylogenetic relationships in Rosaceae have long been problematic because of frequent hybridisation, apomixis and presumed rapid radiation, and their historical diversification has not been clarified. With 87 genera representing all subfamilies and tribes of Rosaceae and six of the other eight families of Rosales (outgroups), we analysed 130 newly sequenced plastomes together with 12 from GenBank in an attempt to reconstruct deep relationships and reveal temporal diversification of this family. Our results highlight the importance of improving sequence alignment and the use of appropriate substitution models in plastid phylogenomics. Three subfamilies and 16 tribes (as previously delimited) were strongly supported as monophyletic, and their relationships were fully resolved and strongly supported at most nodes. Rosaceae were estimated to have originated during the Late Cretaceous with evidence for rapid diversification events during several geological periods. The major lineages rapidly diversified in warm and wet habits during the Late Cretaceous, and the rapid diversification of genera from the early Oligocene onwards occurred in colder and drier environments. Plastid phylogenomics offers new and important insights into deep phylogenetic relationships and the diversification history of Rosaceae. The robust phylogenetic backbone and time estimates we provide establish a framework for future comparative studies on rosaceous evolution.

Resolving the relationships of Paleocene placental mammals

The 'Age of Mammals' began in the Paleocene epoch, the 10 million year interval immediately following the Cretaceous-Palaeogene mass extinction. The apparently rapid shift in mammalian ecomorphs from small, largely insectivorous forms to many small-to-large-bodied, diverse taxa has driven a hypothesis that the end-Cretaceous heralded an adaptive radiation in placental mammal evolution. However, the affinities of most Paleocene mammals have remained unresolved, despite significant advances in understanding the relationships of the extant orders, hindering efforts to reconstruct robustly the origin and early evolution of placental mammals. Here we present the largest cladistic analysis of Paleocene placentals to date, from a data matrix including 177 taxa (130 of which are Palaeogene) and 680 morphological characters. We improve the resolution of the relationships of several enigmatic Paleocene clades, including families of 'condylarths'. Protungulatum is resolved as a stem eutherian, meaning that no crown-placental mammal unambiguously pre-dates the Cretaceous-Palaeogene boundary. Our results support an Atlantogenata-Boreoeutheria split at the root of crown Placentalia, the presence of phenacodontids as closest relatives of Perissodactyla, the validity of Euungulata, and the placement of Arctocyonidae close to Carnivora. Periptychidae and Pantodonta are resolved as sister taxa, Leptictida and Cimolestidae are found to be stem eutherians, and Hyopsodontidae is highly polyphyletic. The inclusion of Paleocene taxa in a placental phylogeny alters interpretations of relationships and key events in mammalian evolutionary history. Paleocene mammals are an essential source of data for understanding fully the biotic dynamics associated with the end-Cretaceous mass extinction. The relationships presented here mark a critical first step towards accurate reconstruction of this important interval in the evolution of the modern fauna.

Identification of secreted and membrane-bound bat immunoglobulin using a Microchiropteran-specific mouse monoclonal antibody

Bat immunity has received increasing attention because some bat species are being decimated by the fungal disease, White Nose Syndrome, while other species are potential reservoirs of zoonotic viruses. Identifying specific immune processes requires new specific tools and reagents. In this study, we describe a new mouse monoclonal antibody (mAb) reactive with Eptesicus fuscus immunoglobulins. The epitope recognized by mAb BT1-4F10 was localized to immunoglobulin light (lambda) chains; hence, the mAb recognized serum immunoglobulins and B lymphocytes. The BT1-4F10 epitope appeared to be restricted to Microchiropteran immunoglobulins and absent from Megachiropteran immunoglobulins. Analyses of sera and other E. fuscus fluids showed that most, if not all, secreted immunoglobulins utilized lambda light chains. Finally, mAb BT1-4F10 permitted the identification of B cell follicles in splenic white pulp. This Microchiropteran-specific mAb has potential utility in seroassays; hence, this reagent may have both basic and practical applications for studying immune process.

Evolution of Digestive Enzymes and RNASE1 Provides Insights into Dietary Switch of Cetaceans

Although cetaceans (whales, porpoises, and dolphins) have multi-chambered stomachs, feeding habits of modern cetaceans have dramatically changed from herbivorous to carnivorous. However, the genetic basis underlying this dietary switch remains unexplored. Here, we present the first systematic investigation of 10 digestive enzymes genes (i.e., CYP7A1, CTRC, LIPC, LIPF, PNLIP, PGC, PRSS1, SI, SLC5A1, and TMPRSS15) of representative cetaceans, and the evolutionary trajectory of RNASE1 in cetartiodactylans. Positive selections were detected with proteinases (i.e., CTRC, PRSS1, and TMPRSS15) and lipases (i.e., CYP7A1, LIPF, and PNLIP) suggesting that cetaceans have evolved an enhanced digestion capacity for proteins and lipids, the major nutritional components of their prey (fishes and invertebrates). In addition, it was found that RNASE1 gene duplicated after the cetartiodactylan speciation and two independent gene duplication events took place in Camelidae and Ruminantia. Positive selection was detected with RNASE1 of Camelidae and Bovidae, suggesting enhanced digestive efficiency in the ruminants. Remarkably, even though the ancestors of cetaceans were terrestrial artiodactyls that are herbivorous, modern cetaceans lost the pancreatic RNASE1 copy with digestive function, which is in accordance with the dietary change from herbivorous to carnivorous. In sum, this is the first study that provides new insights into the evolutionary mechanism of dietary switch in cetaceans.

Dietary specialization drives multiple independent losses and gains in the bitter taste gene repertoire of Laurasiatherian Mammals

Background

Bitter taste perception is essential for species with selective food intake, enabling them to avoid unpalatable or toxic items. Previous studies noted a marked variation in the number of TAS2R genes among various vertebrate species, but the underlying causes are not well understood. Laurasiatherian mammals have highly diversified dietary niche, showing repeated evolution of specialized feeding preferences in multiple lineages and offering a unique chance to investigate how various feeding niches are associated with copy number variation for bitter taste receptor genes.

Results

Here we investigated the evolutionary trajectories of TAS2Rs and their implications on bitter taste perception in whole-genome assemblies of 41 Laurasiatherian species. The number of intact TAS2Rs copies varied considerably, ranging from 0 to 52. As an extreme example of a narrow dietary niche, the Chinese pangolin possessed the lowest number of intact TAS2Rs (n = 2) among studied terrestrial vertebrates. Marine mammals (cetacea and pinnipedia), which swallow prey whole, presented a reduced copy number of TAS2Rs (n = 0-5). In contrast, independent insectivorous lineages, such as the shrew and insectivorous bats possessed a higher TAS2R diversity (n = 52 and n = 20-32, respectively), exceeding that in herbivores (n = 9-22) and omnivores (n = 18-22).

Conclusions

Besides herbivores, insectivores in Laurasiatheria tend to have more functional TAS2Rs in comparison to carnivores and omnivores. Furthermore, animals swallowing food whole (cetacean, pinnipedia and pangolin) have lost most functional TAS2Rs. These findings provide the most comprehensive view of the bitter taste gene repertoire in Laurasiatherian mammals to date, casting new light on the relationship between losses and gains of TAS2Rs and dietary specialization in mammals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12983-016-0161-1) contains supplementary material, which is available to authorized users.

Molecular Evolution of the Nuclear Factor (Erythroid-Derived 2)-Like 2 Gene Nrf2 in Old World Fruit Bats (Chiroptera: Pteropodidae)

Mammals developed antioxidant systems to defend against oxidative damage in their daily life. Enzymatic antioxidants and low molecular weight antioxidants (LMWAs) constitute major parts of the antioxidant systems. Nuclear factor (erythroid-derived 2)-like 2 (Nrf2, encoded by the Nrf2 gene) is a central transcriptional regulator, regulating transcription, of many antioxidant enzymes. Frugivorous bats eat large amounts of fruits that contain high levels of LMWAs such as vitamin C, thus, a reliance on LMWAs might greatly reduce the need for antioxidant enzymes in comparison to insectivorous bats. Therefore, it is possible that frugivorous bats have a reduced need for Nrf2 function due to their substantial intake of diet-antioxidants. To test whether the Nrf2 gene has undergone relaxed evolution in fruit-eating bats, we obtained Nrf2 sequences from 16 species of bats, including four Old World fruit bats (Pteropodidae) and one New World fruit bat (Phyllostomidae). Our molecular evolutionary analyses revealed changes in the selection pressure acting on Nrf2 gene and identified seven specific amino acid substitutions that occurred on the ancestral lineage leading to Old World fruit bats. Biochemical experiments were conducted to examine Nrf2 in Old World fruit bats and showed that the amount of catalase, which is regulated by Nrf2, was significantly lower in the brain, heart and liver of Old World fruit bats despite higher levels of Nrf2 protein in Old World fruit bats. Computational predictions suggest that three of these seven amino acid replacements might be deleterious to Nrf2 function. Therefore, the results suggest that Nrf2 gene might have experienced relaxed constraint in Old World fruit bats, however, we cannot rule out the possibility of positive selection. Our study provides the first data on the molecular adaptation of Nrf2 gene in frugivorous bats in compensation to the increased levels of LWMAs from their fruit-diet.

Comparative Transcriptome and Chloroplast Genome Analyses of Two Related Dipteronia Species

Dipteronia (order Sapindales) is an endangered genus endemic to China and has two living species, D.sinensis and D. dyeriana. The plants are closely related to the genus Acer, which is also classified in the order Sapindales. Evolutionary studies on Dipteronia have been hindered by the paucity of information on their genomes and plastids. Here, we used next generation sequencing to characterize the transcriptomes and complete chloroplast genomes of both Dipteronia species. A comparison of the transcriptomes of both species identified a total of 7814 orthologs. Estimation of selection pressures using Ka/Ks ratios showed that only 30 of 5435 orthologous pairs had a ratio significantly >1, i.e., showing positive selection. However, 4041 orthologs had a Ka/Ks < 0.5 (p < 0.05), suggesting that most genes had likely undergone purifying selection. Based on orthologous unigenes, 314 single copy nuclear genes (SCNGs) were identified. Through a combination of de novo and reference guided assembly, plastid genomes were obtained; that of D. sinensis was 157,080 bp and that of D. dyeriana was 157,071 bp. Both plastid genomes encoded 87 protein coding genes, 40 tRNAs, and 8 rRNAs; no significant differences were detected in the size, gene content, and organization of the two plastomes. We used the whole chloroplast genomes to determine the phylogeny of D. sinensis and D. dyeriana and confirmed that the two species were highly divergent. Overall, our study provides comprehensive transcriptomic and chloroplast genomic resources, which will be valuable for future evolutionary studies of Dipteronia.

Gene structure and evolution of transthyretin in the order Chiroptera

Bats are mammals in the order Chiroptera. Although many extensive morphologic and molecular genetics analyses have been attempted, phylogenetic relationships of bats has not been completely resolved. The paraphyly of microbats is of particular controversy that needs to be confirmed. In this study, we attempted to use the nucleotide sequence of transthyretin (TTR) intron 1 to resolve the relationship among bats. To explore its utility, the complete sequences of TTR gene and intron 1 region of bats in Vespertilionidae: genus Eptesicus (Eptesicus fuscus) and genus Myotis (Myotis brandtii, Myotis davidii, and Myotis lucifugus), and Pteropodidae (Pteropus alecto and Pteropus vampyrus) were extracted from the retrieved sequences, whereas those of Rhinoluphus affinis and Scotophilus kuhlii were amplified and sequenced. The derived overall amino sequences of bat TTRs were found to be very similar to those in other eutherians but differed from those in other classes of vertebrates. However, missing of amino acids from N-terminal or C-terminal region was observed. The phylogenetic analysis of amino acid sequences suggested bat and other eutherian TTRs lineal descent from a single most recent common ancestor which differed from those of non-placental mammals and the other classes of vertebrates. The splicing of bat TTR precursor mRNAs was similar to those of other eutherian but different from those of marsupial, bird, reptile and amphibian. Based on TTR intron 1 sequence, the inferred evolutionary relationship within Chiroptera revealed more closely relatedness of R. affinis to megabats than to microbats. Accordingly, the paraphyly of microbats was suggested.

Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi

We studied the evolutionary history of the Parmeliaceae (Lecanoromycetes, Ascomycota), one of the largest families of lichen-forming fungi with complex and variable morphologies, also including several lichenicolous fungi. We assembled a six-locus data set including nuclear, mitochondrial and low-copy protein-coding genes from 293 operational taxonomic units (OTUs). The lichenicolous lifestyle originated independently three times in lichenized ancestors within Parmeliaceae, and a new generic name is introduced for one of these fungi. In all cases, the independent origins occurred c. 24 million yr ago. Further, we show that the Paleocene, Eocene and Oligocene were key periods when diversification of major lineages within Parmeliaceae occurred, with subsequent radiations occurring primarily during the Oligocene and Miocene. Our phylogenetic hypothesis supports the independent origin of lichenicolous fungi associated with climatic shifts at the Oligocene-Miocene boundary. Moreover, diversification bursts at different times may be crucial factors driving the diversification of Parmeliaceae. Additionally, our study provides novel insight into evolutionary relationships in this large and diverse family of lichen-forming ascomycetes.

The evolution of bat nucleic acid-sensing Toll-like receptors

We characterized the nucleic acid-sensing Toll-like receptors (TLR) of a New World bat species, the common vampire bat (Desmodus rotundus), and through a comparative molecular evolutionary approach searched for general adaptation patterns among the nucleic acid-sensing TLRs of eight different bats species belonging to three families (Pteropodidae, Vespertilionidae and Phyllostomidae). We found that the bat TLRs are evolving slowly and mostly under purifying selection and that the divergence pattern of such receptors is overall congruent with the species tree, consistent with the evolution of many other mammalian nuclear genes. However, the chiropteran TLRs exhibited unique mutations fixed in ligand-binding sites, some of which involved nonconservative amino acid changes and/or targets of positive selection. Such changes could potentially modify protein function and ligand-binding properties, as some changes were predicted to alter nucleic acid binding motifs in TLR 9. Moreover, evidence for episodic diversifying selection acting specifically upon the bat lineage and sublineages was detected. Thus, the long-term adaptation of chiropterans to a wide variety of environments and ecological niches with different pathogen profiles is likely to have shaped the evolution of the bat TLRs in an order-specific manner. The observed evolutionary patterns provide evidence for potential functional differences between bat and other mammalian TLRs in terms of resistance to specific pathogens or recognition of nucleic acids in general.