Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of Microtus voles

Phylogenetic reconstruction of the evolutionary history of closely related organisms may be difficult because of the presence of unsorted lineages and of a relatively high proportion of heterozygous sites that are usually not handled well by phylogenetic programs. Genomic data may provide enough fixed polymorphisms to resolve phylogenetic trees, but the diploid nature of sequence data remains analytically challenging. Here, we performed a phylogenomic reconstruction of the evolutionary history of the common vole (Microtus arvalis) with a focus on the influence of heterozygosity on the estimation of intraspecific divergence times. We used genome-wide sequence information from 15 voles distributed across the European range. We provide a novel approach to integrate heterozygous information in existing phylogenetic programs by repeated random haplotype sampling from sequences with multiple unphased heterozygous sites. We evaluated the impact of the use of full, partial, or no heterozygous information for tree reconstructions on divergence time estimates. All results consistently showed four deep and strongly supported evolutionary lineages in the vole data. These lineages undergoing divergence processes split only at the end or after the last glacial maximum based on calibration with radiocarbon-dated paleontological material. However, the incorporation of information from heterozygous sites had a significant impact on absolute and relative branch length estimations. Ignoring heterozygous information led to an overestimation of divergence times between the evolutionary lineages of M. arvalis. We conclude that the exclusion of heterozygous sites from evolutionary analyses may cause biased and misleading divergence time estimates in closely related taxa.

Bayesian mixed models and divergence time estimation of Chinese cavefishes (Cyprinidae : Sinocyclocheilus)

The genus Sinocyclocheilus is distributed in Yun-Gui Plateau and its surrounding region only, within more than 10 cave species showing different degrees of degeneration of eyes and pigmentation with wonderful adaptations. To present, published morphological and molecular phylogenetic hypotheses of Sinocyclocheilus from prior works are very different and the relationships within the genus are still far from clear. We obtained the sequences of cytochrome b (cyt b) and NADH dehydrogenase subunit 4 (ND4) of 34 species within Sinocyclocheilus, which represent the most dense taxon sampling to date. We performed Bayesian mixed models analyses with this data set. Under this phylogenetic framework, we estimated the divergence times of recovered clades using different methods under relaxed molecular clock. Our phyloegentic results supported the monophyly of Sinocyclocheilus and showed that this genus could be subdivided into 6 major clades. In addition, an earlier finding demonstrating the polyphyletic of cave species and the most basal position of S. jii was corroborated. Relaxed divergence-time estimation suggested that Sinocyclocheilus originated at the late Miocene, about 11 million years ago (Ma), which is older than what have been assumed.

An estimate of divergence time of parazoa and eumetazoa and that of cephalochordata and vertebrata by aldolase and triose phosphate isomerase clocks

Previously we suggested that four proteins including aldolase and triose phosphate isomerase (TPI) evolved with approximately constant rates over long periods covering the whole animal phyla. The constant rates of aldolase and TPI evolution were reexamined based on three different models for estimating evolutionary distances, It was shown that the evolutionary rates remain essentially unchanged in comparisons not only between different classes of vertebrates but also between vertebrates and arthropods and even between animals and plants, irrespective of the models used, Thus these enzymes might be useful molecular clocks for inferring divergence times of animal phyla, To know the divergence time of Parazoa and Eumetazoa and that of Cephalochordata and Vertebrata, the aldolase cDNAs from Ephydatia fluviatilis, a freshwater sponge, and the TPI cDNAs from Ephydatia fluviatilis and Branchiostoma belcheri an amphioxus, have been cloned and sequenced, Comparisons of the deduced amino acid sequences of aldolase and TPI from the freshwater sponge with known sequences revealed that the Parazoa-Eumetazoa split occurred about 940 million years ago (Ma) as determined by the average of two proteins and three models, Similarly, the aldolase and TPI clocks suggest that vertebrates and amphioxus last shared a common ancestor around 700 Ma and they possibly diverged shortly after the divergence of deuterostomes and protostomes.

Improving real-time estimates of output and inflation gaps with multiple-vintage models

Real-time estimates of output gaps and inflation gaps differ from the values that are obtained using data available long after the event. Part of the problem is that the data on which the real-time estimates are based is subsequently revised. We show that vector-autoregressive models of data vintages provide forecasts of post-revision values of future observations and of already-released observations capable of improving estimates of output and inflation gaps in real time. Our findings indicate that annual revisions to output and inflation data are in part predictable based on their past vintages.

Quantification of uncertainties associated with space-time estimates of short-term soil CO2 emissions in a sugar cane area

The characterization of soil CO2 emissions (FCO2) is important for the study of the global carbon cycle. This phenomenon presents great variability in space and time, a characteristic that makes attempts at modeling and forecasting FCO2 challenging. Although spatial estimates have been performed in several studies, the association of these estimates with the uncertainties inherent in the estimation procedures is not considered. This study aimed to evaluate the local, spatial, local-temporal and spatial-temporal uncertainties of short-term FCO2 after harvest period in a sugar cane area. The FCO2 was featured in a sampling grid of 60m×60m containing 127 points with minimum separation distances from 0.5 to 10m between points. The FCO2 was evaluated 7 times within a total period of 10 days. The variability of FCO2 was described by descriptive statistics and variogram modeling. To calculate the uncertainties, 300 realizations made by sequential Gaussian simulation were considered. Local uncertainties were evaluated using the probability values exceeding certain critical thresholds, while the spatial uncertainties considering the probability of regions with high probability values together exceed the adopted limits. Using the daily uncertainties, the local-spatial and spatial-temporal uncertainty (Ftemp) was obtained. The daily and mean emissions showed a variability structure that was described by spherical and Gaussian models. The differences between the daily maps were related to variations in the magnitude of FCO2, covering mean values ranging from 1.28±0.11μmolm-2s-1 (F197) to 1.82±0.07μmolm-2s-1 (F195). The Ftemp showed low spatial uncertainty coupled with high local uncertainty estimates. The average emission showed great spatial uncertainty of the simulated values. The evaluation of uncertainties associated with the knowledge of temporal and spatial variability is an important tool for understanding many phenomena over time, such as the quantification of greenhouse gases or the identification of areas with high crop productivity. © 2013 Elsevier B.V.

Systematics and biogeography of the Gondwanan Orthocladiinae (Diptera: Chironomidae)

Restrictions to effective dispersal and gene flow caused by the fragmentation of ancient supercontinents are considered to have driven diversification and speciation on disjunct landmasses globally. Investigating the role that these processes have played in the development of diversity within and among taxa is crucial to understanding the origins and evolution of regional biotas. Within the chironomid (non-biting midge) subfamily Orthocladiinae (Diptera: Chironomidae), a group of genera that are distributed across the austral continents (Australia, New Zealand, South America) have been proposed to represent a relict Gondwanan clade. We used a molecular approach to resolve relationships among taxa with the aim to determine the relative roles that vicariance and dispersal may have played in the evolution of this group. Continental biotas did not form monophyletic groups, in accordance with expectations given existing morphological evidence. Patterns of phylogenetic relationships among taxa did not accord with expected patterns based on the geological sequence of break-up of the Gondwanan supercontinent. Likewise, divergence time estimates, particularly for New Zealand taxa, largely post-dated continental fragmentation and implied instead that several transoceanic dispersal events may have occurred post-vicariance. Passive dispersal of gravid female chironomid adults is the most likely mechanism for transoceanic movement, potentially facilitated by West Wind Drift or anti-cyclone fronts. Estimated timings of divergence among Australian and South American Botryocladius, on the other hand, were congruent with the proposed ages of separation of the two continents from Antarctica. Taken together, these data suggest that a complex relationship between both vicariance and dispersal may explain the evolution of this group. The sampling regime we implemented here was the most intensive yet performed for austral members of the Orthocladiinae and unsurprisingly revealed several novel taxa that will require formal description.

Not drowning, (hand)waving? Molecular phylogenetics, biogeography and evolutionary tempo of the ‘Gondwanan’ midge Stictocladius Edwards (Diptera: Chironomidae)

Many insect clades, especially within the Diptera (true flies), have been considered classically ‘Gondwanan’, with an inference that distributions derive from vicariance of the southern continents. Assessing the role that vicariance has played in the evolution of austral taxa requires testing the location and tempo of diversification and speciation against the well-established predictions of fragmentation of the ancient super-continent. Several early (anecdotal) hypotheses that current austral distributions originate from the breakup of Gondwana derive from studies of taxa within the family Chironomidae (non-biting midges). With the advent of molecular phylogenetics and biogeographic analytical software, these studies have been revisited and expanded to test such conclusions better. Here we studied the midge genus Stictocladius Edwards, from the subfamily Orthocladiinae, which contains austral-distributed clades that match vicariance-based expectations. We resolve several issues of systematic relationships among morphological species and reveal cryptic diversity within many taxa. Time-calibrated phylogenetic relationships among taxa accorded partially with the predicted tempo from geology. For these apparently vagile insects, vicariance-dated patterns persist for South America and Australia. However, as often found, divergence time estimates for New Zealand at c. 50 mya post-date separation of Zealandia from Antarctica and the remainder of Gondwana, but predate the proposed Oligocene ‘drowning’ of these islands. We detail other such ‘anomalous’ dates and suggest a single common explanation rather than stochastic processes. This could involve synchronous establishment following recovery from ‘drowning’ and/or deleteriously warming associated with the mid-Eocene climatic optimum (hence ‘waving’, which refers to cycles of drowning events) plus new availability of topography providing of cool running waters, or all these factors in combination. Alternatively a vicariance explanation remains available, given the uncertain duration of connectivity of Zealandia to Australia–Antarctic–South America via the Lord Howe and Norfolk ridges into the Eocene.

An inordinate fondness for Fusarium: Phylogenetic diversity of fusaria cultivated by ambrosia beetles in the genus Euwallacea on avocado and other plant hosts

Ambrosia beetle fungiculture represents one of the most ecologically and evolutionarily successful symbioses, as evidenced by the 11 independent origins and 3500 species of ambrosia beetles. Here we document the evolution of a clade within Fusarium associated with ambrosia beetles in the genus Euwallacea (Coleoptera: Scolytinae). Ambrosia Fusarium Clade (AFC) symbionts are unusual in that some are plant pathogens that cause significant damage in naive natural and cultivated ecosystems, and currently threaten avocado production in the United States, Israel and Australia. Most AFC fusaria produce unusual clavate macroconidia that serve as a putative food source for their insect mutualists. AFC symbionts were abundant in the heads of four Euwallacea spp., which suggests that they are transported within and from the natal gallery in mandibular mycangia. In a four-locus phylogenetic analysis, the AFC was resolved in a strongly supported monophyletic group within the previously described Cade 3 of the Fusarium solani species complex (FSSC). Divergence-time estimates place the origin of the AFC in the early Miocene similar to 21.2 Mya, which coincides with the hypothesized adaptive radiation of the Xyleborini. Two strongly supported clades within the AFC (Clades A and B) were identified that include nine species lineages associated with ambrosia beetles, eight with Euwallacea spp. and one reportedly with Xyleborus ferrugineus, and two lineages with no known beetle association. More derived lineages within the AFC showed fixation of the clavate (club-shaped) macroconidial trait, while basal lineages showed a mix of clavate and more typical fusiform macroconidia. AFC lineages consisted mostly of genetically identical individuals associated with specific insect hosts in defined geographic locations, with at least three interspecific hybridization events inferred based on discordant placement in individual gene genealogies and detection of recombinant loci. Overall, these data are consistent with a strong evolutionary trend toward obligate symbiosis coupled with secondary contact and interspecific hybridization. (C) 2013 Elsevier Inc. All rights reserved.

Sun skink diversification across the Indian-Southeast Asian biogeographical interface

Aim Widespread, transcontinental vertebrate groups represent ideal systems for biogeographical studies, because they can shed light on a wide range of questions relating to species diversification across the geographical template. We combined extensive geographical and genetic sampling from across multiple biogeographical realms to examine the timing and location of diversification in Asian sun skinks, a clade characterized by problematic species boundaries and a particularly enigmatic evolutionary history. Location Indian subcontinent, the Philippines, Southeast Asia and Sundaland. Methods We sequenced one mitochondrial and nine nuclear genes for most species in the genus Eutropis, and estimated phylogenetic relationships and divergence times using coalescent methods. To investigate the location of diversification events, we also estimated ancestral geographical ranges using several methods. Finally, we explored patterns of genetic diversity within several poorly understood, but widely distributed species. Results Divergence-time estimates indicate that Eutropis began to diversify during the Eocene. Biogeographical reconstructions show that species diversification was associated with dispersal into three biogeographical realms: India, Sundaland and the Philippines. Main conclusions The results of this study clarify several questions related to the evolutionary history of Eutropis, and place them in the context of classic Southeast Asian biogeography. Our study represents one of the first to compile a heavily sampled multilocus dataset ranging across international boundaries in southern Asia that have historically prevented a unified understanding of biogeographical and evolutionary processes involving the Indian subcontinent, mainland southern Asia and the island archipelagos of Southeast Asia.

Independent evolutionary origins of landlocked alewife populations and rapid parallel evolution of phenotypic traits.

Alewife, Alosa pseudoharengus, populations occur in two discrete life-history variants, an anadromous form and a landlocked (freshwater resident) form. Landlocked populations display a consistent pattern of life-history divergence from anadromous populations, including earlier age at maturity, smaller adult body size, and reduced fecundity. In Connecticut (USA), dams constructed on coastal streams separate anadromous spawning runs from lake-resident landlocked populations. Here, we used sequence data from the mtDNA control region and allele frequency data from five microsatellite loci to ask whether coastal Connecticut landlocked alewife populations are independently evolved from anadromous populations or whether they share a common freshwater ancestor. We then used microsatellite data to estimate the timing of the divergence between anadromous and landlocked populations. Finally, we examined anadromous and landlocked populations for divergence in foraging morphology and used divergence time estimates to calculate the rate of evolution for foraging traits. Our results indicate that landlocked populations have evolved multiple times independently. Tests of population divergence and estimates of gene flow show that landlocked populations are genetically isolated, whereas anadromous populations exchange genes. These results support a 'phylogenetic raceme' model of landlocked alewife divergence, with anadromous populations forming an ancestral core from which landlocked populations independently diverged. Divergence time estimates suggest that landlocked populations diverged from a common anadromous ancestor no longer than 5000 years ago and perhaps as recently as 300 years ago, depending on the microsatellite mutation rate assumed. Examination of foraging traits reveals landlocked populations to have significantly narrower gapes and smaller gill raker spacings than anadromous populations, suggesting that they are adapted to foraging on smaller prey items. Estimates of evolutionary rates (in haldanes) indicate rapid evolution of foraging traits, possibly in response to changes in available resources.

Response to Comment on "Whole-genome analyses resolve early branches in the tree of life of modern birds".

Mitchell et al. argue that divergence-time estimates for our avian phylogeny were too young because of an "inappropriate" maximum age constraint for the most recent common ancestor of modern birds and that, as a result, most modern bird orders diverged before the Cretaceous-Paleogene mass extinction event 66 million years ago instead of after. However, their interpretations of the fossil record and timetrees are incorrect.

Molecular dating and biogeography of fig-pollinating wasps

Figs and fig-pollinating wasps are obligate mutualists that have coevolved for over 60 million years. But when and where did pollinating fig wasps (Agaonidae) originate? Some studies suggest that agaonids arose in the Late Cretaceous and the current distribution of fig-wasp faunas can be explained by the break-up of the Gondwanan landmass. However, recent molecular-dating studies suggest divergence time estimates that are inconsistent with the Gondwanan vicariance hypothesis and imply that long distance oceanic dispersal could have been an important process for explaining the current distribution of both figs and fig wasps. Here, we use a combination of phylogenetic and biogeographical data to infer the age, the major period of diversification, and the geographic origin of pollinating fig wasps. Age estimates ranged widely depending on the molecular-dating method used and even when using the same method but with slightly different constraints, making it difficult to assess with certainty a Gondwanan origin of agaonids. The reconstruction of ancestral areas suggests that the most recent common ancestor of all extant fig-pollinating wasps was most likely Asian, although a southern Gondwana origin cannot be rejected. Our analysis also suggests that dispersal has played a more important role in the development of the fig-wasp biota than previously assumed. (C) 2009 Elsevier Inc. All rights reserved.