Dates from the molecular clock: how wrong can we be?

Comparative Expression Profiling and Sequence Characterization of ATP1A1 Gene Associated with Heat Tolerance in Tropically Adapted Cattle

Climate change is an imminent threat to livestock production. One adaptation strategy is selection for heat tolerance. While it is established that the ATP1A1 gene and its product play an important role in the response to many stressors, there has been no attempt to characterize the sequence or to perform expression profiling of the gene in production animals. We undertook a field experiment to compare the expression profiles of ATP1A1 in heat-tolerant Vechur and Kasaragod cattle (Bos taurus indicus) with the profile of a heat-susceptible crossbreed (B. t. taurus × B. t. indicus). The cattle were exposed to heat stress while on pasture in the hot summer season. The environmental stress was quantified using the temperature humidity index (THI), while the heat tolerance of each breed was assessed using a heat tolerance coefficient (HTC). The ATP1A1 mRNA of Vechur cattle was amplified from cDNA and sequenced. The HTC varied significantly between the breeds and with time-of-day (p < 0.01). The breed-time-of-day interaction was also significant (p < 0.01). The relative expression of ATP1A1 differed between heat-tolerant and heat-susceptible breeds (p = 0.02). The expression of ATP1A1 at 08:00, 10:00 and 12:00, and the breed-time-of-day interaction, were not significant. The nucleotide sequence of Vechur ATP1A1 showed 99% homology with the B. t. taurus sequence. The protein sequence showed 98% homology with B. t. taurus cattle and with B. grunniens (yak) and 97.7% homology with Ovis aries (sheep). A molecular clock analysis revealed evidence of divergent adaptive evolution of the ATP1A1 gene favoring climate resilience in Vechur cattle. These findings further our knowledge of the relationship between the ATP1A1 gene and heat tolerance in phenotypically incongruent animals. We propose that ATP1A1 could be used in marker assisted selection (MAS) for heat tolerance.

The germline mutational process in rhesus macaque and its implications for phylogenetic dating

BACKGROUND

Understanding the rate and pattern of germline mutations is of fundamental importance for understanding evolutionary processes.

RESULTS

Here we analyzed 19 parent-offspring trios of rhesus macaques (Macaca mulatta) at high sequencing coverage of ∼76× per individual and estimated a mean rate of 0.77 × 10-8de novo mutations per site per generation (95% CI: 0.69 × 10-8 to 0.85 × 10-8). By phasing 50% of the mutations to parental origins, we found that the mutation rate is positively correlated with the paternal age. The paternal lineage contributed a mean of 81% of the de novo mutations, with a trend of an increasing male contribution for older fathers. Approximately 3.5% of de novo mutations were shared between siblings, with no parental bias, suggesting that they arose from early development (postzygotic) stages. Finally, the divergence times between closely related primates calculated on the basis of the yearly mutation rate of rhesus macaque generally reconcile with divergence estimated with molecular clock methods, except for the Cercopithecoidea/Hominoidea molecular divergence dated at 58 Mya using our new estimate of the yearly mutation rate.

CONCLUSIONS

When compared to the traditional molecular clock methods, new estimated rates from pedigree samples can provide insights into the evolution of well-studied groups such as primates.

Log Transformation Improves Dating of Phylogenies

Phylogenetic trees inferred from sequence data often have branch lengths measured in the expected number of substitutions and therefore, do not have divergence times estimated. These trees give an incomplete view of evolutionary histories since many applications of phylogenies require time trees. Many methods have been developed to convert the inferred branch lengths from substitution unit to time unit using calibration points, but none is universally accepted as they are challenged in both scalability and accuracy under complex models. Here, we introduce a new method that formulates dating as a nonconvex optimization problem where the variance of log-transformed rate multipliers is minimized across the tree. On simulated and real data, we show that our method, wLogDate, is often more accurate than alternatives and is more robust to various model assumptions.

Enrichment in conservative amino acid changes among fixed and standing missense variations in slowly evolving proteins

The process of molecular evolution has many elements that are not yet fully understood. Evolutionary rates are known to vary among protein coding and noncoding DNAs, and most of the observed changes in amino acid or nucleotide sequences are assumed to be non-adaptive by the neutral theory of molecular evolution. However, it remains unclear whether fixed and standing missense changes in slowly evolving proteins are more or less neutral compared to those in fast evolving genes. Here, based on the evolutionary rates as inferred from identity scores between orthologs in human and Rhesus Macaques (Macaca mulatta), we found that the fraction of conservative substitutions between species was significantly higher in their slowly evolving proteins. Similar results were obtained by using four different methods of scoring conservative substitutions, including three that remove the impact of substitution probability, where conservative changes require fewer mutations. We also examined the single nucleotide polymorphisms (SNPs) by using the 1000 Genomes Project data and found that missense SNPs in slowly evolving proteins also had a higher fraction of conservative changes, especially for common SNPs, consistent with more non-conservative substitutions and hence stronger natural selection for SNPs, particularly rare ones, in fast evolving proteins. These results suggest that fixed and standing missense variants in slowly evolving proteins are more likely to be neutral.

Multilocus phylogeny of Paratelmatobiinae (Anura: Leptodactylidae) reveals strong spatial structure and previously unknown diversity in the Atlantic Forest hotspot

The Brazilian Atlantic Forest harbors high levels of anuran diversity and endemism, including several taxa restricted to small geographic ranges. Here, we provide a multilocus phylogeny for Paratelmatobiinae, a leptodactylid subfamily composed of small-ranged species distributed in the Brazilian Atlantic Forest and in the campo rupestre ecosystem. We performed Bayesian inference and maximum likelihood analyses using three mitochondrial and five nuclear markers, and a matrix comprising a broad taxonomic sampling. We then delimitated independently evolving lineages within the group. We recovered Paratelmatobiinae and each of its four genera as monophyletic and robustly supported. Five putatively new species included in our analyses were unambiguously supported in the phylogenetic trees and delimitation analyses. We also recovered other deeply divergent and geographically structured lineages within the four genera of Paratelmatobiinae. Our estimation of divergence times indicates that diversification in the subfamily began in the Eocene and continued until the Pleistocene. We discuss possible scenarios of diversification for the four genera of Paratelmatobiinae, and outline the implications of our findings for taxonomy and conservation.

The coincidence of ecological opportunity with hybridization explains rapid adaptive radiation in Lake Mweru cichlid fishes

The process of adaptive radiation was classically hypothesized to require isolation of a lineage from its source (no gene flow) and from related species (no competition). Alternatively, hybridization between species may generate genetic variation that facilitates adaptive radiation. Here we study haplochromine cichlid assemblages in two African Great Lakes to test these hypotheses. Greater biotic isolation (fewer lineages) predicts fewer constraints by competition and hence more ecological opportunity in Lake Bangweulu, whereas opportunity for hybridization predicts increased genetic potential in Lake Mweru. In Lake Bangweulu, we find no evidence for hybridization but also no adaptive radiation. We show that the Bangweulu lineages also colonized Lake Mweru, where they hybridized with Congolese lineages and then underwent multiple adaptive radiations that are strikingly complementary in ecology and morphology. Our data suggest that the presence of several related lineages does not necessarily prevent adaptive radiation, although it constrains the trajectories of morphological diversification. It might instead facilitate adaptive radiation when hybridization generates genetic variation, without which radiation may start much later, progress more slowly or never occur.

Population genetic structure and geographical variation in Neotricula aperta (Gastropoda: Pomatiopsidae), the snail intermediate host of Schistosoma mekongi (Digenea: Schistosomatidae)

Background

Neotricula aperta is the snail-intermediate host of the parasitic blood-fluke Schistosoma mekongi which causes Mekong schistosomiasis in Cambodia and the Lao PDR. Despite numerous phylogenetic studies only one DNA-sequence based population-genetic study of N. aperta had been published, and the origin, structure and persistence of N. aperta were poorly understood. Consequently, a phylogenetic and population genetic study was performed, with addition of new data to pre-existing DNA-sequences for N. aperta from remote and inaccessible habitats, including one new taxon from Laos and 505 bp of additional DNA-sequence for all sampled taxa,.

Principal findings

Spatial Principal Component Analysis revealed the presence of significant spatial-genetic clustering. Genetic-distance-based clustering indicated four populations with near perfect match to a priori defined ecogeographical regions. Spring-dwelling taxa were found to form an ecological isolate relative to other N. aperta. The poor dispersal capabilities suggested by spatial-genetic analyses were confirmed by Bayesian inference of migration rates. Population divergence time estimation implied a mid-Miocene colonisation of the present range, with immediate and rapid radiation in each ecogeographical region. Estimated effective population sizes were large (120–310 thousand).

Conclusions

The strong spatial-genetic structure confirmed the poor dispersal capabilities of N. aperta—suggesting human-mediated reintroduction of disease to controlled areas as the primary reason for control failure. The isolation of the spring-dwelling taxa and ecogeographical structure suggests adaptation of sub-populations to different habitats; the epidemiological significance of this needs investigation. The large effective population sizes indicate that the high population densities observed in surveyed habitats are also present in inaccessible areas; affording great potential for recrudescence driven by animal-reservoir transmission in remote streams. Mid-Miocene colonisation implies heterochronous evolution of these snails and associated schistosomes and suggests against coevolution of snail and parasite. Heterochronicity favours ecological factors as shapers of host-parasite specificity and greater potential for escape from schistosomiasis control through host-switching.

The disease Mekong schistosomiasis poses a threat to the health of about 1.5 million people living near the Mekong river and its tributaries in Cambodia and Laos. It is a water-borne parasite transmitted by direct contact with water in which freshwater snails of the species Neotricula aperta live. Control of the snails is an effective approach to control of the parasite; however, because many suitable habitats for N. aperta occur in remote and inaccessible areas, knowledge of N. aperta population sizes and interconnectivity is insufficient for the design of effective snail control interventions. Although much of the region is difficult to survey by conventional means, population genetics can be used to estimate population structure and total size from small samples of accessible populations. The study added to existing data-sets, to give more population samples and longer DNA-sequences, together with improved analytical approaches to provide a better overview of N. aperta. The findings suggest that N. aperta in different kinds of habitats are also genetically different, with very low levels of migration between them; this genetic clustering is greater than expected from spatial distance alone. Further work is needed to determine if these different clusters vary in ability to transmit the parasite. The overall population size estimates were very large; thus suggesting that high snail population densities observed in accessible habitats are also characteristic of inaccessible populations—parasites are therefore more likely to return after disease control by immigration from remote areas. Finally, the timing of evolutionary events for snails and parasites was found to differ; this implies that the parasite may not be as strongly restricted to one species of snail as originally thought, which has implications for avoidance of parasite control by host-switching.

Persistence of phylogeographic footprints helps to understand cryptic diversity detected in two marine amphipods widespread in the Mediterranean basin

Amphipods of the genus Gammarus are a vital component of macrozoobenthic communities in European inland and coastal, marine and brackish waters of the Mediterranean and the Black Sea. Exceptional levels of cryptic diversity have been revealed for several widespread freshwater Gammarus species in Europe. No comprehensive assessment has yet been made for brackishwater counterparts, such as Gammarus aequicauda and G. insensibilis, which are among the most widely dispersed members of the so-called "G. locusta group" in the Mediterranean and in the Black Sea. Here we probe the diversity of these morphospecies examining the partitioning of mtDNA and nDNA across multiple populations along their distribution range and discuss it within the regional paleogeographic framework. We gathered molecular data from a collection of 166 individuals of G. aequicauda and G. insensibilis from 47 locations along their distribution range in the Mediterranean including the Black Sea. They were amplified for both mitochondrial COI and 16S rRNA as well as the nuclear 28S rRNA. All five MOTU delimitation methods (ABGD, BIN, bPTP, GMYC single and multiple threshold models) applied revealed deep divergence between Black Sea and Mediterranean populations in both G. aequicauda and G. insensibilis. There were eight distinct MOTUs delimited for G. aequicauda (6-18% K2P) and 4 MOTUs for G. insensibilis (4-14% K2P). No sympatric MOTUs were detected throughout their distribution range. Multimarker time-calibrated phylogeny indicated that divergence of both G. aequicauda and G. insensibilis species complexes started already in the late Oligocene/early Miocene with the split between clades inhabiting eastern and western part of the Mediterranean occurring in both species at the similar time. Our results indicate a high cryptic diversity within Mediterranean brackishwater Gammarus, similar to that observed for freshwater counterparts. Moreover, the phylogenetic history combined with the current geographic distribution indicate that the evolution of both studied Gammarus morphogroups has been strongly connected with the geological events in the Mediterranean Basin and it reflect the turbulent history of the area.

Molecular Phylogeny and Dating of Forsythieae (Oleaceae) Provide Insight into the Miocene History of Eurasian Temperate Shrubs

Tribe Forsythieae (Oleaceae), containing two genera (Abeliophyllum and Forsythia) and 13 species, is economically important plants used as ornamentals and in traditional medicine. This tribe species occur primarily in mountainous regions of Eurasia with the highest species diversity in East Asia. Here, we examine 11 complete chloroplast genome and nuclear cycloidea2 (cyc2) DNA sequences of 10 Forsythia species and Abeliophyllum distichum using Illumina platform to provide the phylogeny and biogeographic history of the tribe. The chloroplast genomes of the 11 Forsythieae species are highly conserved, except for a deletion of about 400 bp in the accD-psaI region detected only in Abeliophyllum. Within Forsythieae species, analysis of repetitive sequences revealed a total of 51 repeats comprising 26 forward repeats, 22 palindromic repeats, and 3 reverse repeats. Of those, 19 repeats were common and 32 were unique to one or more Forsythieae species. Our phylogenetic analyses supported the monophyly of Forsythia and its sister group is Abeliophyllum using the concatenated dataset of 78 chloroplast genes. Within Forsythia, Forsythia likiangensis and F. giraldiana were basal lineages followed by F. europaea; the three species are characterized by minutely serrate or entire leaf margins. The remaining species, which are distributed in East Asia, formed two major clades. One clade included F. ovata, F. velutina, and F. japonica; they are morphologically supported by broadly ovate leaves. Another clade of F. suspensa, F. saxatilis, F. viridissima, and F. koreana characterized by lanceolate leaves (except F. suspensa which have broad ovate leaves). Although cyc2 phylogeny is largely congruent to chloroplast genome phylogeny, we find the discordance between two phylogenies in the position of F. ovata suggesting that introgression of the chloroplast genome from one species into the nuclear background of another by interspecific hybridization in East Asian Forsythia species. Molecular dating and biogeographic reconstructions suggest an origin of the Forsythieae species in East China in the Miocene. Distribution patterns in Forsythia indicated that the species were radially differentiated from East China, and the speciation of the European F. europaea was the result of both vicariance and dispersal in the late Miocene to Pliocene.

Genetic equidistance at nucleotide level

The genetic equidistance phenomenon shows complex taxa to be all approximately equidistant to a less complex species in amino acid percentage identity. The overlooked mystery was re-interpreted by the maximum genetic diversity hypothesis (MGD). Here, we studied 14 proteomes and their coding DNA sequences (CDS) to see if the equidistance phenomenon also holds at the CDS level. We found that the outgroup taxon was equidistant to the two more complex taxa species. When two sister taxa were compared to human as the outgroup, the more complex taxon was closer to human, confirming species complexity to be the primary determinant of MGD. Finally, we found the fraction of overlap sites to be inversely correlated with CDS conservation, indicating saturation to be more common in less conserved DNAs. These results establish the genetic equidistance phenomenon to be universal at the DNA level and provide additional evidence for the MGD theory.

A Possible Role for Stochastic Astrophysical Ionizing Radiation Events in the Systematic Disparity between Molecular and Fossil Dates

Major discrepancies have been noted for some time between fossil ages and molecular divergence dates for a variety of taxa. Recently, systematic trends within avian clades have been uncovered. The trends show that the disparity is much larger for mitochondrial DNA than for nuclear DNA, also that it is larger for crown fossil dates than stem fossil dates. It has been argued that this pattern is largely inconsistent with incompleteness of the fossil record as the principal driver of the disparity. A case is presented that, given the expected mutations from a fluctuating background of astrophysical radiation from such sources as supernovae, the rate of molecular clocks is variable and should increase back to a few million years, before returning to the long-term average rate. This is a possible explanation for the disparity. One test of this hypothesis is to look for an acceleration of molecular clocks at 2 to 2.5 Ma due to one or more moderately nearby supernovae known to have happened at that time. Another is to look for reduced disparity in benthic organisms of the deep ocean. In addition, due to the importance of highly penetrating muon irradiation, the disparity should be magnified for megafauna. Key Words: Extreme events in Earth history-Molecular clock-Radiation physics-Evolution. Astrobiology 17, 87-90.

Phytochrome Interacting Factors (PIFs) in Solanum lycopersicum: Diversity, Evolutionary History and Expression Profiling during Different Developmental Processes

Although the importance of light for tomato plant yield and edible fruit quality is well known, the PHYTOCHROME INTERACTING FACTORS (PIFs), main components of phytochrome-mediated light signal transduction, have been studied almost exclusively in Arabidopsis thaliana. Here, the diversity, evolution and expression profile of PIF gene subfamily in Solanum lycopersicum was characterized. Eight tomato PIF loci were identified, named SlPIF1a, SlPIF1b, SlPIF3, SlPIF4, SlPIF7a, SlPIF7b, SlPIF8a and SlPIF8b. The duplication of SlPIF1, SlPIF7 and SlPIF8 genes were dated and temporally coincided with the whole-genome triplication event that preceded tomato and potato divergence. Different patterns of mRNA accumulation in response to light treatments were observed during seedling deetiolation, dark-induced senescence, diel cycle and fruit ripening. SlPIF4 showed similar expression profile as that reported for A. thaliana homologs, indicating an evolutionary conserved function of PIF4 clade. A comprehensive analysis of the evolutionary and transcriptional data allowed proposing that duplicated SlPIFs have undergone sub- and neofunctionalization at mRNA level, pinpointing the importance of transcriptional regulation for the maintenance of duplicated genes. Altogether, the results indicate that genome polyploidization and functional divergence have played a major role in diversification of the Solanum PIF gene subfamily.

Human Germline Mutation and the Erratic Evolutionary Clock

Our understanding of the chronology of human evolution relies on the “molecular clock” provided by the steady accumulation of substitutions on an evolutionary lineage. Recent analyses of human pedigrees have called this understanding into question by revealing unexpectedly low germline mutation rates, which imply that substitutions accrue more slowly than previously believed. Translating mutation rates estimated from pedigrees into substitution rates is not as straightforward as it may seem, however. We dissect the steps involved, emphasizing that dating evolutionary events requires not “a mutation rate” but a precise characterization of how mutations accumulate in development in males and females—knowledge that remains elusive.

Comparative Microbial Genomics and Forensics

Forensic science concerns the application of scientific techniques to questions of a legal nature and may also be used to address questions of historical importance. Forensic techniques are often used in legal cases that involve crimes against persons or property, and they increasingly may involve cases of bioterrorism, crimes against nature, medical negligence, or tracing the origin of food- and crop-borne disease. Given the rapid advance of genome sequencing and comparative genomics techniques, we ask how these might be used to address cases of a forensic nature, focusing on the use of microbial genome sequence analysis. Such analyses rely on the increasingly large numbers of microbial genomes present in public databases, the ability of individual investigators to rapidly sequence whole microbial genomes, and an increasing depth of understanding of their evolution and function. Suggestions are made as to how comparative microbial genomics might be applied forensically and may represent possibilities for the future development of forensic techniques. A particular emphasis is on the nascent field of genomic epidemiology, which utilizes rapid whole-genome sequencing to identify the source and spread of infectious outbreaks. Also discussed is the application of comparative microbial genomics to the study of historical epidemics and deaths and how the approaches developed may also be applicable to more recent and actionable cases.

A Phylogeographic Assessment of the Malagasy Giant Chameleons (Furcifer verrucosus and Furcifer oustaleti)

The Malagasy giant chameleons (Furcifer oustaleti and Furcifer verrucosus) are sister species that are both broadly distributed in Madagascar, and also endemic to the island. These species are also morphologically similar and, because of this, have been frequently misidentified in the field. Previous studies have suggested that cryptic species are nested within this chameleon group, and two subspecies have been described in F. verrucosus. In this study, we utilized a phylogeographic approach to assess genetic diversification within these chameleons. This was accomplished by (1) identifying clades within each species supported by both mitochondrial and nuclear DNA, (2) assessing divergence times between clades, and (3) testing for niche divergence or conservatism. We found that both F. oustaleti and F. verrucosus could be readily identified based on genetic data, and within each species, there are two well-supported clades. However, divergence times are not contemporary and spatial patterns are not congruent. Diversification within F. verrucosus occurred during the Plio-Pleistocene, and there is evidence for niche divergence between a southwestern and southeastern clade, in a region of Madagascar that shows no obvious landscape barriers to dispersal. Diversification in F. oustaleti occurred earlier in the Pliocene or Miocene, and niche conservatism is supported with two genetically distinct clades separated at the Sofia River in northwestern Madagascar. Divergence within F. verrucosus is most consistent with patterns expected from ecologically mediated speciation, whereas divergence in F. oustaleti most strongly matches the patterns expected from the riverine barrier hypothesis.

If Dung Beetles (Scarabaeidae: Scarabaeinae) Arose in Association with Dinosaurs, Did They Also Suffer a Mass Co-Extinction at the K-Pg Boundary?

The evolutionary success of beetles and numerous other terrestrial insects is generally attributed to co-radiation with flowering plants but most studies have focused on herbivorous or pollinating insects. Non-herbivores represent a significant proportion of beetle diversity yet potential factors that influence their diversification have been largely unexamined. In the present study, we examine the factors driving diversification within the Scarabaeidae, a speciose beetle family with a range of both herbivorous and non-herbivorous ecologies. In particular, it has been long debated whether the key event in the evolution of dung beetles (Scarabaeidae: Scarabaeinae) was an adaptation to feeding on dinosaur or mammalian dung. Here we present molecular evidence to show that the origin of dung beetles occurred in the middle of the Cretaceous, likely in association with dinosaur dung, but more surprisingly the timing is consistent with the rise of the angiosperms. We hypothesize that the switch in dinosaur diet to incorporate more nutritious and less fibrous angiosperm foliage provided a palatable dung source that ultimately created a new niche for diversification. Given the well-accepted mass extinction of non-avian dinosaurs at the Cretaceous-Paleogene boundary, we examine a potential co-extinction of dung beetles due to the loss of an important evolutionary resource, i.e., dinosaur dung. The biogeography of dung beetles is also examined to explore the previously proposed “out of Africa” hypothesis. Given the inferred age of Scarabaeinae as originating in the Lower Cretaceous, the major radiation of dung feeders prior to the Cenomanian, and the early divergence of both African and Gondwanan lineages, we hypothesise that that faunal exchange between Africa and Gondwanaland occurred during the earliest evolution of the Scarabaeinae. Therefore we propose that both Gondwanan vicariance and dispersal of African lineages is responsible for present day distribution of scarabaeine dung beetles and provide examples.

New thoughts on an old riddle: What determines genetic diversity within and between species?

The question of what determines genetic diversity has long remained unsolved by the modern evolutionary theory (MET). However, it has not deterred researchers from producing interpretations of genetic diversity by using MET. We examine the two observations of genetic diversity made in the 1960s that contributed to the development of MET. The interpretations of these observations by MET are widely known to be inadequate. We review the recent progress of an alternative framework, the maximum genetic diversity (MGD) hypothesis, that uses axioms and natural selection to explain the vast majority of genetic diversity as being at equilibrium that is largely determined by organismal complexity. The MGD hypothesis absorbs the proven virtues of MET and considers its assumptions relevant only to a much more limited scope. This new synthesis has accounted for the overlooked phenomenon of progression towards higher complexity, and more importantly, been instrumental in directing productive research.

Out-of-the tropics or trans-tropical dispersal? The origins of the disjunct distribution of the gooseneck barnacle Pollicipes elegans

Background

Studying species with disjunct distributions allows biogeographers to evaluate factors controlling species ranges, limits on gene flow, and allopatric speciation. Here, we use phylogeographic and population genetic studies of the barnacle Pollicipes elegans to discriminate between two primary hypotheses about the origin of disjunct distributions of extra-tropical populations: trans-tropical stepping-stone colonization versus an out-of-the tropics origin.

Results

Nucleotide diversity peaked in the centre of the species’ range in samples from El Salvador and was lower in samples from higher latitudes at Mexico and Peru. Haplotypes from El Salvador samples also had a deeper coalescent, or an older time to a most recent common ancestor. A deep phylogeographical break exists between Mexico and all samples taken to the south (El Salvador and Peru). Isolation-with-migration analyses showed no significant gene flow between any of the three regions indicating that the difference in genetic differentiation among all three regions is explained primarily by differences in population separation times. Approximate Bayesian Computation model testing found strong support for an out-of-the tropics origin of extra-tropical populations in P. elegans.

Conclusions

We found little evidence consistent with a stepping-stone history of trans-tropical colonization, but instead found strong evidence for a tropical origin model for the largely disjunct distribution of P. elegans. Sea surface temperature and habitat suitability are likely mechanisms driving decline of populations in tropical regions, causing the disjunct distribution.

Electronic supplementary material

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

Phylogeny of Dictyoptera: Dating the Origin of Cockroaches, Praying Mantises and Termites with Molecular Data and Controlled Fossil Evidence

Understanding the origin and diversification of organisms requires a good phylogenetic estimate of their age and diversification rates. This estimate can be difficult to obtain when samples are limited and fossil records are disputed, as in Dictyoptera. To choose among competing hypotheses of origin for dictyopteran suborders, we root a phylogenetic analysis (~800 taxa, 10 kbp) within a large selection of outgroups and calibrate datings with fossils attributed to lineages with clear synapomorphies. We find the following topology: (mantises, (other cockroaches, (Cryptocercidae, termites)). Our datings suggest that crown-Dictyoptera-and stem-mantises-would date back to the Late Carboniferous (~ 300 Mya), a result compatible with the oldest putative fossil of stem-dictyoptera. Crown-mantises, however, would be much more recent (~ 200 Mya; Triassic/Jurassic boundary). This pattern (i.e., old origin and more recent diversification) suggests a scenario of replacement in carnivory among polyneopterous insects. The most recent common ancestor of (cockroaches + termites) would date back to the Permian (~275 Mya), which contradicts the hypothesis of a Devonian origin of cockroaches. Stem-termites would date back to the Triassic/Jurassic boundary, which refutes a Triassic origin. We suggest directions in extant and extinct species sampling to sharpen this chronological framework and dictyopteran evolutionary studies.

Benthic perspective on Earth's oldest evidence for oxygenic photosynthesis

The Great Oxidation Event (GOE) is currently viewed as a protracted process during which atmospheric oxygen increased above ∼10(-5) times the present atmospheric level (PAL). This threshold represents an estimated upper limit for sulfur isotope mass-independent fractionation (S-MIF), an Archean signature of atmospheric anoxia that begins to disappear from the rock record at 2.45 Ga. However, an increasing number of papers have suggested that the timing for oxidative continental weathering, and by conventional thinking the onset of atmospheric oxygenation, was hundreds of million years earlier than previously thought despite the presence of S-MIF. We suggest that this apparent discrepancy can be resolved by the earliest oxidative-weathering reactions occurring in benthic and soil environments at profound redox disequilibrium with the atmosphere, such as biological soil crusts and freshwater microbial mats covering riverbed, lacustrine, and estuarine sediments. We calculate that oxygenic photosynthesis in these millimeter-thick ecosystems provides sufficient oxidizing equivalents to mobilize sulfate and redox-sensitive trace metals from land to the oceans while the atmosphere itself remained anoxic with its attendant S-MIF signature. As continental freeboard increased significantly between 3.0 and 2.5 Ga, the chemical and isotopic signatures of benthic oxidative weathering would have become more globally significant from a mass-balance perspective. These observations help reconcile evidence for pre-GOE oxidative weathering with the history of atmospheric chemistry, and support the plausible antiquity of a terrestrial biosphere populated by cyanobacteria well before the GOE.

Transcriptome analysis reveals positive selection on the divergent between topmouth culter and zebrafish

The topmouth culter (Erythroculter ilishaeformis) is a predatory cyprinid fish that distributes widely in the East Asia. Here we report the liver transcriptome in this organism as a model of predatory fish. Sequencing of 5 Gb raw reads led to 27,741 unigenes and produced 11,131 annotatable genes. A total of 7093 (63.7%) genes were found to have putative functions by gene ontology analysis. Importantly, a blast search revealed 4033 culter genes that were orthologous to the zebrafish. Extracted from 38 candidate positive selection genes, 4 genes exhibit strong positive selection based on the ratio of nonsynonymous (Ka) to synonymous substitutions (Ks). In addition, the four genes also indicated the strong positive selection by comparing them between blunt snout bream (Megalobrama amblycephala) and zebrafish. These genes were involved in activator of gene expression, metabolic processes and development. The transcriptome variation may be reflective of natural selection in the early life history of Cyprinidae. Based on Ks ratios, date of the separation between topmouth culter and zebrafish is approximately 64 million years ago. We conclude that natural selection acts in diversifying the genomes between topmouth culter and zebrafish.

Phylogenetic informativeness reconciles ray-finned fish molecular divergence times

Background

Discordance among individual molecular age estimates, or between molecular age estimates and the fossil record, is observed in many clades across the Tree of Life. This discordance is attributed to a variety of variables including calibration age uncertainty, calibration placement, nucleotide substitution rate heterogeneity, or the specified molecular clock model. However, the impact of changes in phylogenetic informativeness of individual genes over time on phylogenetic inferences is rarely analyzed. Using nuclear and mitochondrial sequence data for ray-finned fishes (Actinopterygii) as an example, we extend the utility of phylogenetic informativeness profiles to predict the time intervals when nucleotide substitution saturation results in discordance among molecular ages estimated.

Results

We demonstrate that even with identical calibration regimes and molecular clock methods, mitochondrial based molecular age estimates are systematically older than those estimated from nuclear sequences. This discordance is most severe for highly nested nodes corresponding to more recent (i.e., Jurassic-Recent) divergences. By removing data deemed saturated, we reconcile the competing age estimates and highlight that the older mtDNA based ages were driven by nucleotide saturation.

Conclusions

Homoplasious site patterns in a DNA sequence alignment can systematically bias molecular divergence time estimates. Our study demonstrates that PI profiles can provide a non-arbitrary criterion for data exclusion to mitigate the influence of homoplasy on time calibrated branch length estimates. Analyses of actinopterygian molecular clocks demonstrate that scrutiny of the time scale on which sequence data is informative is a fundamental, but generally overlooked, step in molecular divergence time estimation.

The changing face of the molecular evolutionary clock

The molecular clock has played an important role in biological research, both as a description of the evolutionary process and as a tool for inferring evolutionary timescales. Genomic data have provided valuable insights into the molecular clock, allowing the patterns and causes of evolutionary rate variation to be characterized in increasing detail. I explain how genome sequences offer exciting opportunities for estimating the timescale of the Tree of Life. I describe the different approaches that have been used to deal with the computational and statistical challenges encountered in molecular clock analyses of genomic data. Finally, I offer a perspective on the future of molecular clocks, highlighting some of the key limitations and the most promising research directions.

Speciation and demographic history of Atlantic eels (Anguilla anguilla and A. rostrata) revealed by mitogenome sequencing

Processes leading to speciation in oceanic environments without obvious physical barriers remain poorly known. European and American eel (Anguilla anguilla and A. rostrata) spawn in partial sympatry in the Sargasso Sea. Larvae are advected by the Gulf Stream and other currents towards the European/North African and North American coasts, respectively. We analyzed 104 mitogenomes from the two species along with mitogenomes of other Anguilla and outgroup species. We estimated divergence time between the two species to identify major events involved in speciation. We also considered two previously stated hypotheses: one where the ancestral species was present in only one continent but was advected across the Atlantic by ocean current changes and another where population declines during Pleistocene glaciations led to increasing vicariance, facilitating speciation. Divergence time was estimated to ∼3.38 Mya, coinciding with the closure of the Panama Gateway that led to reinforcement of the Gulf Stream. This could have advected larvae towards European/North African coasts, in which case American eel would be expected to be the ancestral species. This scenario could, however, not be unequivocally confirmed by analyses of dN/dS, nucleotide diversity and effective population size estimates. Extended bayesian skyline plots showed fluctuations of effective population sizes and declines during glaciations, and thus also lending support to the importance of vicariance during speciation. There was evidence for positive selection at the ATP6 and possibly ND5 genes, indicating a role in speciation. The findings suggest an important role of ocean current changes in speciation of marine organisms.

A phylogeny for the pomatiopsidae (Gastropoda: Rissooidea): a resource for taxonomic, parasitological and biodiversity studies

BACKGROUND

The Pomatiopsidae are reported from northern India into southern China and Southeast Asia, with two sub-families, the Pomatiopsinae (which include freshwater, amphibious, terrestrial and marine species) and the freshwater Triculinae. Both include species acting as intermediate host for species of the blood-fluke Schistosoma which cause a public health problem in East Asia. Also, with around 120 species, triculine biodiversity exceeds that of any other endemic freshwater molluscan fauna. Nevertheless, the origins of the Pomatiopsidae, the factors driving such a diverse radiation and aspects of their co-evolution with Schistosoma are not fully understood. Many taxonomic questions remain; there are problems identifying medically relevant species. The predicted range is mostly unsurveyed and the true biodiversity of the family is underestimated. Consequently, the aim of the study was to collect DNA-sequence data for as many pomatiopsid taxa as possible, as a first step in providing a resource for identification of epidemiologically significant species (by non-malacologists), for use in resolving taxonomic confusion and for testing phylogeographical hypotheses.

RESULTS

The evolutionary radiation of the Triculinae was shown to have been rapid and mostly post late Miocene. Molecular dating indicated that the radiation of these snails was driven first by the uplift of the Himalaya and onset of a monsoon system, and then by late-Pliocene global warming. The status of Erhaia as Anmicolidae is supported. The genera Tricula and Neotricula are shown to be non-monophyletic and the tribe Jullieniini may be polyphyletic (based on convergent characters). Triculinae from northern Vietnam could be derived from Gammatricula of Fujian/Yunnan, China.

CONCLUSIONS

The molecular dates and phylogenetic estimates in this study are consistent with an Australasian origin for the Pomatiopsidae and an East to West radiation via Oligocene Borneo-Philippines island hopping to Japan and then China (Triculinae arising mid-Miocene in Southeast China), and less so with a triculine origin in Tibet. The lack of monophyly in the medically important genera and indications of taxonomic inaccuracies, call for further work to identify epidemiologically significant taxa (e.g., Halewisia may be potential hosts for Schistosoma mekongi) and highlight the need for surveys to determine the true biodiversity of the Triculinae.

Molecular phylogeny and SNP variation of polar bears (Ursus maritimus), brown bears (U. arctos), and black bears (U. americanus) derived from genome sequences

We assessed the relationships of polar bears (Ursus maritimus), brown bears (U. arctos), and black bears (U. americanus) with high throughput genomic sequencing data with an average coverage of 25× for each species. A total of 1.4 billion 100-bp paired-end reads were assembled using the polar bear and annotated giant panda (Ailuropoda melanoleuca) genome sequences as references. We identified 13.8 million single nucleotide polymorphisms (SNP) in the 3 species aligned to the polar bear genome. These data indicate that polar bears and brown bears share more SNP with each other than either does with black bears. Concatenation and coalescence-based analysis of consensus sequences of approximately 1 million base pairs of ultraconserved elements in the nuclear genome resulted in a phylogeny with black bears as the sister group to brown and polar bears, and all brown bears are in a separate clade from polar bears. Genotypes for 162 SNP loci of 336 bears from Alaska and Montana showed that the species are genetically differentiated and there is geographic population structure of brown and black bears but not polar bears.

Cryptic genetic diversity is paramount in small-bodied amphibians of the genus Euparkerella (Anura: Craugastoridae) endemic to the Brazilian Atlantic forest

Morphological similarity associated to restricted distributions and low dispersal abilities make the direct developing “Terrarana” frogs of the genus Euparkerella a good model for examining diversification processes. We here infer phylogenetic relationships within the genus Euparkerella, using DNA sequence data from one mitochondrial and four nuclear genes coupled with traditional Bayesian phylogenetic reconstruction approaches and more recent coalescent methods of species tree inference. We also used Bayesian clustering analysis and a recent Bayesian coalescent-based approach specifically to infer species delimitation. The analysis of 39 individuals from the four known Euparkerella species uncovered high levels of genetic diversity, especially within the two previously morphologically-defined E. cochranae and E. brasiliensis. Within these species, the gene trees at five independent loci and trees from combined data (concatenated dataset and the species tree) uncovered six deeply diverged and geographically coherent evolutionary units, which may have diverged between the Miocene and the Pleistocene. These six units were also uncovered in the Bayesian clustering analysis, and supported by the Bayesian coalescent-based species delimitation (BPP), and Genealogical Sorting Index (GSI), providing thus strong evidence for underestimation of the current levels of diversity within Euparkerella. The cryptic diversity now uncovered opens new opportunities to examine the origins and maintenance of microendemism in the context of spatial heterogeneity and/or human induced fragmentation of the highly threatened Brazilian Atlantic forest hotspot.

A multi-calibrated mitochondrial phylogeny of extant Bovidae (Artiodactyla, Ruminantia) and the importance of the fossil record to systematics

BACKGROUND

Molecular phylogenetics has provided unprecedented resolution in the ruminant evolutionary tree. However, molecular age estimates using only one or a few (often misapplied) fossil calibration points have produced a diversity of conflicting ages for important evolutionary events within this clade. I here identify 16 fossil calibration points of relevance to the phylogeny of Bovidae and Ruminantia and use these, individually and together, to construct a dated molecular phylogeny through a reanalysis of the full mitochondrial genome of over 100 ruminant species.

RESULTS

The new multi-calibrated tree provides ages that are younger overall than found in previous studies. Among these are young ages for the origin of crown Ruminantia (39.3-28.8 Ma), and crown Bovidae (17.3-15.1 Ma). These are argued to be reasonable hypotheses given that many basal fossils assigned to these taxa may in fact lie on the stem groups leading to the crown clades, thus inflating previous age estimates. Areas of conflict between molecular and fossil dates do persist, however, especially with regard to the base of the rapid Pecoran radiation and the sister relationship of Moschidae to Bovidae. Results of the single-calibrated analyses also show that a very wide range of molecular age estimates are obtainable using different calibration points, and that the choice of calibration point can influence the topology of the resulting tree. Compared to the single-calibrated trees, the multi-calibrated tree exhibits smaller variance in estimated ages and better reflects the fossil record.

CONCLUSIONS

The use of a large number of vetted fossil calibration points with soft bounds is promoted as a better approach than using just one or a few calibrations, or relying on internal-congruency metrics to discard good fossil data. This study also highlights the importance of considering morphological and ecological characteristics of clades when delimiting higher taxa. I also illustrate how phylogeographic and paleoenvironmental hypotheses inferred from a tree containing only extant taxa can be problematic without consideration of the fossil record. Incorporating the fossil record of Ruminantia is a necessary step for future analyses aiming to reconstruct the evolutionary history of this clade.

The genetic equidistance result: misreading by the molecular clock and neutral theory and reinterpretation nearly half of a century later

In 1963, Margoliash discovered the unexpected genetic equidistance result after comparing cytochrome c sequences from different species. This finding, together with the hemoglobin analyses of Zuckerkandl and Pauling in 1962, directly inspired the ad hoc molecular clock hypothesis. Unfortunately, however, many biologists have since mistakenly viewed the molecular clock as a genuine reality, which in turn inspired Kimura, King, and Jukes to propose the neutral theory of molecular evolution. Many years of studies have found numerous contradictions to the theory, and few today believe in a universal constant clock. What is being neglected, however, is that the failure of the molecular clock hypothesis has left the original equidistance result an unsolved mystery. In recent years, we fortuitously rediscovered the equidistance result, which remains unknown to nearly all researchers. Incorporating the proven virtues of existing evolutionary theories and introducing the novel concept of maximum genetic diversity, we proposed a more complete hypothesis of evolutionary genetics and reinterpreted the equidistance result and other major evolutionary phenomena. The hypothesis may rewrite molecular phylogeny and population genetics and solve major biomedical problems that challenge the existing framework of evolutionary biology.