Molecular phylogeny and evolution of the freshwater eel, genus Anguilla

A molecular phylogenetic analysis was conducted on all of the known catadromous eel species of the genus Anguilla to assess their relationships and evolutionary history. The analyses of a total of 1427 bp of the mitochondrial 16S ribosomal RNA and 1140 bp of the complete cytochrome b gene sequences suggested that the genus Anguilla was monophyletic in origin, with A. borneensis as the most basal species. Four clades/species groups that correspond to their geographical ranges were indicated, Indo-Atlantic (three species), Oceania (two species), tropical Pacific (two species), and Indo-Pacific (five species), with ambiguous positions for A. japonica and A. reinhardti. This grouping conflicts with that of a previous morphological study, since the broad undivided maxillary and short-fin type, which were thought to be phylogenetically important, were paraphyletic in the molecular analysis. However, the molecular phylogeny and the present geographic distribution of species suggested historical dispersion of the genus Anguilla according to the Tethys corridor hypothesis, which proposed that anguillid eels originated near present-day Indonesia and dispersed westward along paleo-circumglobal equatorial currents. The westward-moving strain entered the paleo-Atlantic through the Tethys Sea and was ancestral to present-day European and American species.     

Directional selection in the evolution of elongated upper canines in clouded leopards and sabre‐toothed cats

Extremely developed or specialized traits such as the elongated upper canines of extinct sabre‐toothed cats are often not analogous to those of any extant species, which limits our understanding of their evolutionary cause. However, an extant species may have undergone directional selection for a similar extreme phenotype. Among living felids, the clouded leopard, Neofelis nebulosa, has exceptionally long upper canines for its body size. We hypothesized that directional selection generated the elongated upper canines of clouded leopards in a manner similar to the process in extinct sabre‐toothed cats. To test this, we developed an approach that compared the effect of directional selection among lineages in a phylogeny using a simulation of trait evolution and approximate Bayesian computation. This approach was applied to analyse the evolution of upper canine length in the Felidae phylogeny. Our analyses consistently showed directional selection favouring longer upper canines in the clouded leopard lineage and a lineage leading to the sabre‐toothed cat with the longest upper canines, Smilodon. Most of our analyses detected an effect of directional selection for longer upper canines in the lineage leading to another sabre‐toothed cat, Homotherium, although this selection may have occurred exclusively in the primitive species. In all the analyses, the clouded leopard and Smilodon lineages showed comparable directional selection. This implies that clouded leopards share a selection advantage with sabre‐toothed cats in having elongated upper canines.     

Morphological evolution of silica scales in the freshwater genus Synura (Stramenopiles)

A high degree of morphological variability is expressed between the ornately sculptured siliceous scales formed by species in the chrysophycean genus, Synura. In this study, we aimed to uncover the general principles and trends underlying the evolution of scale morphology in this genus. We assessed the relationships among thirty extant Synura species using a robust molecular analysis that included six genes, coupled with morphological characterization of the species‐specific scales. The analysis was further enriched with addition of morphological information from fossil specimens and by including the unique modern species, Synura punctulosa. We inferred the phylogenetic position of the morphologically unique S. punctulosa, to be an ancient Synura lineage related to S. splendida in the section Curtispinae. Some morphological traits, including development of a keel or a labyrinth ribbing pattern on the scale, appeared once in evolution, whereas other structures, such as a hexagonal meshwork pattern, originated independently several times over geologic time. We further uncovered numerous construction principles governing scale morphology and evolution, as follows: (i) scale roundness and pore diameter decreased during evolution; (ii) elongated scales became strengthened by a higher number of struts or ribs; (iii) as a consequence of scale biogenesis, scales with spines possessed smaller basal holes than scales with a keel and; and (iv) the keel area was proportional to scale area, indicating its potential value in strengthening the scale against breakage.     

Molecular phylogeny of the diatom genera Amphora and Halamphora (Bacillariophyta) with a focus on morphological and ecological evolution

The taxonomic history of the diatom genus Amphora is one of a broad early morphological concept resulting in the inclusion of a diversity of taxa, followed by an extended period of revision and refinement. The introduction of molecular systematics has increased the pace of revision and has largely resolved the relationships between the major lineages, indicating homoplasy in the evolution of amphoroid symmetry. Within the two largest monophyletic lineages, the genus Halamphora and the now taxonomically refined genus Amphora, the intrageneric morphological and ecological relationships have yet to be explored within a phylogenetic framework. Critical among this is whether the range of morphological features exhibited within these diverse genera are reflective of evolutionary groupings or, as with many previously studied amphoroid features, are nonhomologous when examined phylogenetically. Presented here is a four‐marker molecular phylogeny that includes 31 taxa from the genus Amphora and 77 taxa from the genus Halamphora collected from fresh, brackish, and salt waters from coastal and inland habitats of the United States and Japan. These phylogenies illustrate complex patterns in the evolution of frustule morphology and ecology within the genera and the implications of this on the taxonomy, classification, and organization of the genera are discussed.     

Colonisation of toxic environments drives predictable life‐history evolution in livebearing fishes (Poeciliidae)

New World livebearing fishes (family Poeciliidae) have repeatedly colonised toxic, hydrogen sulphide‐rich waters across their natural distribution. Physiological considerations and life‐history theory predict that these adverse conditions should favour the evolution of larger offspring. Here, we examined nine poeciliid species that independently colonised toxic environments, and show that these fishes have indeed repeatedly evolved much larger offspring size at birth in sulphidic waters, thus uncovering a widespread pattern of predictable evolution. However, a second pattern, only indirectly predicted by theory, proved additionally common: a reduction in the number of offspring carried per clutch (i.e. lower fecundity). Our analyses reveal that this secondary pattern represents a mere consequence of a classic life‐history trade‐off combined with strong selection on offspring size alone. With such strong natural selection in extreme environments, extremophile organisms may commonly exhibit multivariate phenotypic shifts even though not all diverging traits necessarily represent adaptations to the extreme conditions.     

Phylogeny and species delineation in the marine diatom Pseudo‐nitzschia (Bacillariophyta) using cox1, LSU,and ITS2 rRNA genes: A perspective in character evolution

Analyses of the mitochondrial cox1, the nuclear‐encoded large subunit (LSU), and the internal transcribed spacer 2 (ITS2) RNA coding region of Pseudo‐nitzschia revealed that the P. pseudodelicatissima complex can be phylogenetically grouped into three distinct clades (Groups I–III), while the P. delicatissima complex forms another distinct clade (Group IV) in both the LSU and ITS2 phylogenetic trees. It was elucidated that comprehensive taxon sampling (sampling of sequences), selection of appropriate target genes and outgroup, and alignment strategies influenced the phylogenetic accuracy. Based on the genetic divergence, ITS2 resulted in the most resolved trees, followed by cox1 and LSU. The morphological characters available for Pseudo‐nitzschia, although limited in number, were overall in agreement with the phylogenies when mapped onto the ITS2 tree. Information on the presence/absence of a central nodule, number of rows of poroids in each stria, and of sectors dividing the poroids mapped onto the ITS2 tree revealed the evolution of the recently diverged species. The morphologically based species complexes showed evolutionary relevance in agreement with molecular phylogeny inferred from ITS2 sequence–structure data. The data set of the hypervariable region of ITS2 improved the phylogenetic inference compared to the cox1 and LSU data sets. The taxonomic status of P. cuspidata and P. pseudodelicatissima requires further elucidation.     

Molecular analyses identify hybridization‐mediated nuclear evolution in newly discovered fungal hybrids

Hybridization may be a major driver in the evolution of plant pathogens. In a high elevation Alpine larch stand in Montana, a novel hybrid fungal pathogen of trees originating from the mating of Heterobasidion irregulare with H. occidentale has been recently discovered. In this study, sequence analyses of one mitochondrial and four nuclear loci from 11 Heterobasidion genotypes collected in the same Alpine larch stand indicated that hybridization has increased allelic diversity by generating novel polymorphisms unreported in either parental species. Sequence data and ploidy analysis through flow cytometry confirmed that heterokaryotic (n + n) genotypes were not first‐generation hybrids, but were the result of multiple backcrosses, indicating hybrids are fertile. Additionally, all admixed genotypes possessed the H. occidentale mitochondrion, indicating that the hybrid progeny may have been backcrossing mostly with H. occidentale. Based on reticulate phylogenetic network analysis by PhyloNet, Bayesian assignment, and ordination tests, alleles can be defined as H. irregulare‐like or H. occidentale‐like. H. irregulare‐like alleles are clearly distinct from all known H. irregulare alleles and are derived from the admixing of both Heterobasidion species. Instead, all but one H. occidentale alleles found in hybrids, although novel, were not clearly distinct from alleles found in the parental H. occidentale population. This discovery demonstrates that Alpine larch can be a universal host favouring the interspecific hybridization between H. irregulare and H. occidentale and the hybridization‐mediated evolution of a nucleus, derived from H. irregulare parental species but clearly distinct from it.     

Multi‐locus phylogeny of African pipits and longclaws (Aves: Motacillidae) highlights taxonomic inconsistencies

The globally distributed avian family Motacillidae consists of five to seven genera (Anthus, Dendronanthus, Tmetothylacus, Macronyx and Motacilla, and depending on the taxonomy followed, Amaurocichla and Madanga) and 66–68 recognized species, of which 32 species in four genera occur in sub‐Saharan Africa. The taxonomy of the Motacillidae has been contentious, with variable numbers of genera, species and subspecies proposed and some studies suggesting greater taxonomic diversity than currently recognized (five genera and 67 species). Using one nuclear (Mb) and two mitochondrial (cyt b and CO1) gene regions amplified from DNA extracted from contemporary and museum specimens, we investigated the taxonomic status of 56 of the currently recognized motacillid species and present the most taxonomically complete and expanded phylogeny of this family to date. Our results suggest that the family comprises six clades broadly reflecting continental distributions: sub‐Saharan Africa (two clades), the New World (one clade), Palaearctic (one clade), a widespread large‐bodied Anthus clade, and a sixth widespread genus, Motacilla. Within the Afrotropical region, our phylogeny further supports recognition of Wood Pipit Anthus nyassae as a valid species, and the treatment of Long‐tailed Pipit Anthus longicaudatus and Kimberley Pipit Anthus pseudosimilis as junior subjective synonyms of Buffy Pipit Anthus vaalensis and African Pipit Anthus cinnamomeus, respectively. As the disjunct populations of Long‐billed Pipit Anthus similis in southern and East Africa are genetically distinct and geographically separated, we propose a specific status for the southern African population under the earliest available name, Nicholson's Pipit Anthus nicholsoni. Further, as our analyses indicate that Yellow‐breasted Pipit Anthus chloris and Golden Pipit Tmetothylacus tenellus are both nested within the Macronyx longclaws, we propose transferring these species to the latter genus.     

Temperature‐dependent phagotrophy and phototrophy in a mixotrophic chrysophyte

The roles of temperature and light on grazing and photosynthesis were examined for Dinobryon sociale, a common freshwater mixotrophic alga. Photosynthetic rate was determined for D. sociale adapted to temperatures of 8, 12, 16, and 20°C under photosynthetically active radiation light irradiances of 25, 66, and 130 μmol photons · m^?2 · s^?1, with concurrent measurement of bacterial ingestion at all temperatures under medium and high light (66 and 130 μmol photons · m^?2 · s^?1). Rates of ingestion and photosynthesis increased with temperature to a maximum at 16°C under the two higher light regimes, and declined at 20°C. Although both light and temperature had a marked effect on photosynthesis, there was no significant difference in bacterivory at medium and high irradiances at any given temperature. At the lowest light condition (25 μmol photons · m^?2 · s^?1), photosynthesis remained low and relatively stable at all temperatures. D. sociale acquired the majority of carbon from photosynthesis, although the low photosynthetic rate without a concurrent decline in feeding rate at 8°C suggested 20%–30% of the carbon budget could be attributed to bacterivory at low temperatures. Grazing experiments in nutrient‐modified media revealed that this mixotroph had increased ingestion rates when either dissolved nitrogen or phosphorus was decreased. This work increases our understanding of environmental effects on mixotrophic nutrition. Although the influence of abiotic factors on phagotrophy and phototrophy in pure heterotrophs and phototrophs has been well studied, much less is known for mixotrophic organisms.     

DNA‐based phylogeny of the marine genus Heterodrilus (Annelida,Clitellata, Naididae)

Heterodrilus is a group of marine Naididae, common worldwide in subtropical and tropical areas, and unique among the oligochaetes by their tridentate chaetae. The phylogenetic relationships within the group are assessed from the nuclear 18S rDNA gene, and the mitochondrial cytochrome c oxidase subunit I (COI) and 16S rDNA genes. Sequence data were obtained from 16 Heterodrilus species and 13 out‐group taxa; 48 sequences are new for this study. The data were analysed by Bayesian inference. Monophyly of the genus is corroborated by the resulting tree, with Heterodrilus ersei (a taxon representing a small group of species with aberrant male genitalia) proposed to be outside all other sampled species. Although earlier regarded as a member of the subfamily Rhyacodrilinae, both molecular and morphological data seem to support that Heterodrilus is closely related to Phallodrilinae. However, the results are not conclusive as to whether the genus is the sister group of, or a group nested inside, or separate from this latter subfamily. The studied sample of species suggests at least two major clades in Heterodrilus with different geographical distributions, in one of the clades, most species are from the Indo‐West Pacific Ocean, while in the other, the majority are from the Western Atlantic Ocean. Morphological characters traditionally used in Heterodrilus taxonomy are optimized on the phylogenetic tree, revealing a high degree of homoplasy.     

SIMULATION‐BASED LIKELIHOOD APPROACH FOR EVOLUTIONARY MODELS OF PHENOTYPIC TRAITS ON PHYLOGENY

Phylogenetic comparative methods (PCMs) have been used to test evolutionary hypotheses at phenotypic levels. The evolutionary modes commonly included in PCMs are Brownian motion (genetic drift) and the Ornstein–Uhlenbeck process (stabilizing selection), whose likelihood functions are mathematically tractable. More complicated models of evolutionary modes, such as branch‐specific directional selection, have not been used because calculations of likelihood and parameter estimates in the maximum‐likelihood framework are not straightforward. To solve this problem, we introduced a population genetics framework into a PCM, and here, we present a flexible and comprehensive framework for estimating evolutionary parameters through simulation‐based likelihood computations. The method does not require analytic likelihood computations, and evolutionary models can be used as long as simulation is possible. Our approach has many advantages: it incorporates different evolutionary modes for phenotypes into phylogeny, it takes intraspecific variation into account, it evaluates full likelihood instead of using summary statistics, and it can be used to estimate ancestral traits. We present a successful application of the method to the evolution of brain size in primates. Our method can be easily implemented in more computationally effective frameworks such as approximate Bayesian computation (ABC), which will enhance the use of computationally intensive methods in the study of phenotypic evolution.     

Establishing the phylogeny of Prochlorococcus with a new alignment‐free method

Prochlorococcus marinus, one of the most abundant marine cyanobacteria in the global ocean, is classified into low‐light (LL) and high‐light (HL) adapted ecotypes. These two adapted ecotypes differ in their ecophysiological characteristics, especially whether adapted for growth at high‐light or low‐light intensities. However, some evolutionary relationships of Prochlorococcus phylogeny remain to be resolved, such as whether the strains SS120 and MIT9211 form a monophyletic group. We use the Natural Vector (NV) method to represent the sequence in order to identify the phylogeny of the Prochlorococcus. The natural vector method is alignment free without any model assumptions. This study added the covariances of amino acids in protein sequence to the natural vector method. Based on these new natural vectors, we can compute the Hausdorff distance between the two clades which represents the dissimilarity. This method enables us to systematically analyze both the dataset of ribosomal proteomes and the dataset of 16s‐23s rRNA sequences in order to reconstruct the phylogeny of Prochlorococcus. Furthermore, we apply classification to inspect the relationship of SS120 and MIT9211. From the reconstructed phylogenetic trees and classification results, we may conclude that the SS120 does not cluster with MIT9211. This study demonstrates a new method for performing phylogenetic analysis. The results confirm that these two strains do not form a monophyletic clade in the phylogeny of Prochlorococcus.     

Deeply conserved susceptibility in a multi‐host,multi‐parasite system

Variation in susceptibility is ubiquitous in multi‐host, multi‐parasite assemblages, and can have profound implications for ecology and evolution in these systems. The extent to which susceptibility to parasites is phylogenetically conserved among hosts can be revealed by analysing diverse regional communities. We screened for haemosporidian parasites in 3983 birds representing 40 families and 523 species, spanning ~ 4500 m elevation in the tropical Andes. To quantify the influence of host phylogeny on infection status, we applied Bayesian phylogenetic multilevel models that included a suite of environmental, spatial, temporal, life history and ecological predictors. We found evidence of deeply conserved susceptibility across the avian tree; host phylogeny explained substantial variation in infection status, and results were robust to phylogenetic uncertainty. Our study suggests that susceptibility is governed, in part, by conserved, latent aspects of anti‐parasite defence. This demonstrates the importance of deep phylogeny for understanding present‐day ecological interactions.     

Phylogenomic analysis of gene co‐expression networks reveals the evolution of functional modules

Molecular evolutionary studies correlate genomic and phylogenetic information with the emergence of new traits of organisms. These traits are, however, the consequence of dynamic gene networks composed of functional modules, which might not be captured by genomic analyses. Here, we established a method that combines large‐scale genomic and phylogenetic data with gene co‐expression networks to extensively study the evolutionary make‐up of modules in the moss Physcomitrella patens, and in the angiosperms Arabidopsis thaliana and Oryza sativa (rice). We first show that younger genes are less annotated than older genes. By mapping genomic data onto the co‐expression networks, we found that genes from the same evolutionary period tend to be connected, whereas old and young genes tend to be disconnected. Consequently, the analysis revealed modules that emerged at a specific time in plant evolution. To uncover the evolutionary relationships of the modules that are conserved across the plant kingdom, we added phylogenetic information that revealed duplication and speciation events on the module level. This combined analysis revealed an independent duplication of cell wall modules in bryophytes and angiosperms, suggesting a parallel evolution of cell wall pathways in land plants. We provide an online tool allowing plant researchers to perform these analyses at http://www.gene2function.de .     

Trade‐offs and evolution of thermal adaptation in the Irish potato famine pathogen Phytophthora infestans

Temperature is one of the most important environmental parameters with crucial impacts on nearly all biological processes. Due to anthropogenic activity, average air temperatures are expected to increase by a few degrees in coming decades, accompanied by an increased occurrence of extreme temperature events. Such global trends are likely to have various major impacts on human society through their influence on natural ecosystems, food production and biotic interactions, including diseases. In this study, we used a combination of statistical genetics, experimental evolution and common garden experiments to investigate the evolutionary potential for thermal adaptation in the potato late blight pathogen, Phytophthora infestans, and infer its likely response to changing temperatures. We found a trade‐off associated with thermal adaptation to heterogeneous environments in P. infestans, with the degree of the trade‐off peaking approximately at the pathogen's optimum growth temperature. A genetic trade‐off in thermal adaptation was also evidenced by the negative association between a strain's growth rate and its thermal range for growth, and warm climates selecting for a low pathogen growth rate. We also found a mirror effect of phenotypic plasticity and genetic adaptation on growth rate. At below the optimum, phenotypic plasticity enhances pathogen's growth rate but nature selects for slower growing genotypes when temperature increases. At above the optimum, phenotypic plasticity reduces pathogen's growth rate but natural selection favours for faster growing genotypes when temperature increases further. We conclude from these findings that the growth rate of P. infestans will only be marginally affected by global warming.     

Optimizing the widely used nuclear protein‐coding gene primers in beetle phylogenies and their application in the genus Sasajiscymnus Vandenberg (Coleoptera: Coccinellidae)

Advances in genomic biology and the increasing availability of genomic resources allow developing hundreds of nuclear protein‐coding (NPC) markers, which can be used in phylogenetic research. However, for low taxonomic levels, it may be more practical to select a handful of suitable molecular loci for phylogenetic inference. Unfortunately, the presence of degenerate primers of NPC markers can be a major impediment, as the amplification success rate is low and they tend to amplify nontargeted regions. In this study, we optimized five NPC fragments widely used in beetle phylogenetics (i.e., two parts of carbamoyl‐phosphate synthetase: CADXM and CADMC, Topoisomerase, Wingless and Pepck) by reducing the degenerate site of primers and the length of target genes slightly. These five NPC fragments and 6 other molecular loci were amplified to test the monophyly of the coccinellid genus Sasajiscymnus Vandenberg. The analysis of our molecular data set clearly supported the genus Sasajiscymnus may be monophyletic but confirmation with an extended sampling is required. A fossil‐calibrated chronogram was generated by BEAST, indicating an origin of the genus at the end of the Cretaceous (77.87 Myr). Furthermore, a phylogenetic informativeness profile was generated to compare the phylogenetic properties of each gene more explicitly. The results showed that COI provides the strongest phylogenetic signal among all the genes, but Pepck, Topoisomerase, CADXM and CADMC are also relatively informative. Our results provide insight into the evolution of the genus Sasajiscymnus, and also enrich the molecular data resources for further study.     

Parallel evolution of larval morphology and habitat in the snail-killing fly genus Tetanocera

In this study, we sequenced one nuclear and three mitochondrial DNA loci to construct a robust estimate of phylogeny for all available species of Tetanocera. Character optimizations suggested that aquatic habitat was the ancestral condition for Tetanocera larvae, and that there were at least three parallel transitions to terrestrial habitat, with one reversal. Maximum likelihood analyses of character state transformations showed significant correlations between habitat transitions and changes in four larval morphological characteristics (cuticular pigmentation and three characters associated with the posterior spiracular disc). We provide evidence that phylogenetic niche conservatism has been responsible for the maintenance of aquatic-associated larval morphological character states, and that concerted convergence and/or gene linkage was responsible for parallel morphological changes that were derived in conjunction with habitat transitions. These habitat-morphology associations were consistent with the action of natural selection in facilitating the morphological changes that occurred during parallel aquatic to terrestrial habitat transitions in Tetanocera.     

A pre‐Miocene Irano‐Turanian cradle: Origin and diversification of the species‐rich monocot genus Gagea (Liliaceae)

The Irano‐Turanian (IT) floristic region is considered an important center of origin for many taxa. However, there is a lack of studies dealing with typical IT genera that also occur in neighboring areas. The species‐rich monocot genus Gagea Salisb. shows a center of diversity in IT region and a distribution in adjacent regions, therefore representing a good study object to investigate spatial and temporal relationships among IT region and its neighboring areas (East Asia, Euro‐Siberia, Himalaya, and Mediterranean). We aimed at (a) testing the origin of the genus and of its major lineages in the IT region, (b) reconstructing divergence times, and (c) reconstructing colonization events. To address these problems, sequences of the ribosomal DNA internal transcribed spacer (ITS) region of 418 individuals and chloroplast intergenic spacers sequences (psbA‐trnH, trnL‐trnF) of 497 individuals, representing 116 species from all sections of the genus and nearly its entire distribution area were analyzed. Divergence times were estimated under a random molecular clock based on nrITS phylogeny, which was the most complete data set regarding the representation of species and distribution areas. Ancestral distribution ranges were estimated for the nrITS data set as well as for a combined data set, revealing that Gagea most likely originated in southwestern Asia. This genus first diversified there starting in the Early Miocene. In the Middle Miocene, Gagea migrated to the Mediterranean and to East Asia, while migration into Euro‐Siberia took place in the Late Miocene. During the Pleistocene, the Arctic was colonized and Gagea serotina, the most widespread species, reached North America. The Mediterranean basin was colonized multiple times from southwestern Asia or Euro‐Siberia. Most of the currently existing species originated during the last 3 Ma.