BMC-EVOL.BIOL.

Background: We have analyzed the complete mitochondrial genomes of 22 Pan paniscus (bonobo, pygmy chimpanzee) individuals to assess the detailed mitochondrial DNA (mtDNA) phylogeny of this close relative of Homo sapiens. Results: We identified three major clades among bonobos that separated approximately 540,000 years ago, as suggested by Bayesian analysis. Incidentally, we discovered that the current reference sequence for bonobo likely is a hybrid of the mitochondrial genomes of two distant individuals. When comparing spectra of polymorphic mtDNA sites in bonobos and humans, we observed two major differences: (i) Of all 31 bonobo mtDNA homoplasies, i.e. nucleotide changes that occurred independently on separate branches of the phylogenetic tree, 13 were not homoplasic in humans. This indicates that at least a part of the unstable sites of the mitochondrial genome is species-specific and difficult to be explained on the basis of a mutational hotspot concept. (ii) A comparison of the ratios of non-synonymous to synonymous changes (dN/dS) among polymorphic positions in bonobos and in 4902 Homo sapiens mitochondrial genomes revealed a remarkable difference in the strength of purifying selection in the mitochondrial genes of the F0F1-ATPase complex. While in bonobos this complex showed a similar low value as complexes I and IV, human haplogroups displayed 2.2 to 7.6 times increased dN/dS ratios when compared to bonobos. Conclusions: Some variants of mitochondrially encoded subunits of the ATPase complex in humans very likely decrease the efficiency of energy conversion leading to production of extra heat. Thus, we hypothesize that the species-specific release of evolutionary constraints for the mitochondrial genes of the proton-translocating ATPase is a consequence of altered heat homeostasis in modern humans.

Background: Patterns of spatial variation in discrete phenotypic traits can be used to draw inferences about the adaptive significance of traits and evolutionary processes, especially when compared to patterns of neutral genetic variation. Population divergence in adaptive traits such as color morphs can be influenced by both local ecology and stochastic factors such as genetic drift or founder events. Here, we use quantitative color measurements of males and females of Skyros wall lizard, Podarcis gaigeae, to demonstrate that this species is polymorphic with respect to throat color, and the morphs form discrete phenotypic clusters with limited overlap between categories. We use divergence in throat color morph frequencies and compare that to neutral genetic variation to infer the evolutionary processes acting on islet- and mainland populations. Results: Geographically close islet- and mainland populations of the Skyros wall lizard exhibit strong divergence in throat color morph frequencies. Population variation in throat color morph frequencies between islets was higher than that between mainland populations, and the effective population sizes on the islets were small (Ne:s < 100). Population divergence (FST) for throat color morph frequencies fell within the neutral FST-distribution estimated from microsatellite markers, and genetic drift could thus not be rejected as an explanation for the pattern. Moreover, for both comparisons among mainland-mainland population pairs and between mainland-islet population pairs, morph frequency divergence was significantly correlated with neutral divergence, further pointing to some role for genetic drift in divergence also at the phenotypic level of throat color morphs. Conclusions: Genetic drift could not be rejected as an explanation for the pattern of population divergence in morph frequencies. In spite of an expected stabilising selection, throat color frequencies diverged in the islet populations. These results suggest that there is an interaction between selection and genetic drift causing divergence even at a phenotypic level in these small, subdivided populations.

Background: The grape powdery mildew fungus, Erysiphe necator, was introduced into Europe more than 160 years ago and is now distributed everywhere that grapes are grown. To understand the invasion history of this pathogen we investigated the evolutionary relationships between introduced populations of Europe, Australia and the western United States (US) and populations in the eastern US, where E. necator is thought to be native. Additionally, we tested the hypothesis that populations of E. necator in the eastern US are structured based on geography and Vitis host species. Results: We sequenced three nuclear gene regions covering 1803 nucleotides from 146 isolates of E. necator collected from the eastern US, Europe, Australia, and the western US. Phylogeographic analyses show that the two genetic groups in Europe represent two separate introductions and that the genetic groups may be derived from eastern US ancestors. Populations from the western US and Europe share haplotypes, suggesting that the western US population was introduced from Europe. Populations in Australia are derived from European populations. Haplotype richness and nucleotide diversity were significantly greater in the eastern US populations than in the introduced populations. Populations within the eastern US are geographically differentiated; however, no structure was detected with respect to host habitat (i.e., wild or cultivated). Populations from muscadine grapes, V. rotundifolia, are genetically distinct from populations from other Vitis host species, yet no differentiation was detected among populations from other Vitis species. Conclusions: Multilocus sequencing analysis of the grape powdery mildew fungus is consistent with the hypothesis that populations in Europe, Australia and the western US are derived from two separate introductions and their ancestors were likely from native populations in the eastern US. The invasion history of E. necator follows a pattern consistent with plant-mediated dispersal, however, more exhaustive sampling is required to make more precise conclusions as to origin. E. necator shows no genetic structure across Vitis host species, except with respect to V. rotundifolia.

Background: The human FOXI1 gene codes for a transcription factor involved in the physiology of the inner ear, testis, and kidney. Using three interspecies comparisons, it has been suggested that this may be a gene under human-specific selection. We sought to confirm this finding by using an extended set of orthologous sequences. Additionally, we explored for signals of natural selection within humans by sequencing the gene in 20 Europeans, 20 East Asians and 20 Yorubas and by analysing SNP variation in a 2 Mb region centered on FOXI1 in 39 worldwide human populations from the HGDP-CEPH diversity panel. Results: The genome sequences recently available from other primate and non-primate species showed that FOXI1 divergence patterns are compatible with neutral evolution. Sequence-based neutrality tests were not significant in Europeans, East Asians or Yorubas. However, the Long Range Haplotype (LRH) test, as well as the iHS and XP-Rsb statistics revealed significantly extended tracks of homozygosity around FOXI1 in Africa, suggesting a recent episode of positive selection acting on this gene. A functionally relevant SNP, as well as several SNPs either on the putatively selected core haplotypes or with significant iHS or XP-Rsb values, displayed allele frequencies strongly correlated with the absolute geographical latitude of the populations sampled. Conclusions: We present evidence for recent positive selection in the FOXI1 gene region in Africa. Climate might be related to this recent adaptive event in humans. Of the multiple functions of FOXI1, its role in kidney-mediated water-electrolyte homeostasis is the most obvious candidate for explaining a climate-related adaptation.

Background: Ecological speciation is a process in which a transiently resource-polymorphic species divides into two specialized sister lineages as a result of divergent selection pressures caused by the use of multiple niches or environments. Ecology-based speciation has been studied intensively in plant-feeding insects, in which both sympatric and allopatric shifts onto novel host plants could speed up diversification. However, while numerous examples of species pairs likely to have originated by resource shifts have been found, the overall importance of ecological speciation in relation to other, non-ecological speciation modes remains unknown. Here, we apply phylogenetic information on sawflies belonging to the 'Higher' Nematinae (Hymenoptera: Tenthredinidae) to infer the frequency of niche shifts in relation to speciation events. Results: Phylogenetic trees reconstructed on the basis of DNA sequence data show that the diversification of higher nematines has involved frequent shifts in larval feeding habits and in the use of plant taxa. However, the inferred number of resource shifts is considerably lower than the number of past speciation events, indicating that the majority of divergences have occurred by non-ecological allopatric speciation; based on a time-corrected analysis of sister species, we estimate that a maximum of c. 20% of lineage splits have been triggered by a change in resource use. In addition, we find that postspeciational changes in geographic distributions have led to broad sympatry in many species having identical host-plant ranges. Conclusion: Our analysis indicates that the importance of niche shifts for the diversification of herbivorous insects is at present implicitly and explicitly overestimated. In the case of the Higher Nematinae, employing a time correction for sister-species comparisons lowered the proportion of apparent ecology-based speciation events from c. 50-60% to around 20%, but such corrections are still lacking in other herbivore groups. The observed convergent but asynchronous shifting among dominant northern plant taxa in many higher-nematine clades, in combination with the broad overlaps in the geographic distributions of numerous nematine species occupying near-identical niches, indicates that host-plant shifts and herbivore community assembly are largely unconstrained by direct or indirect competition among species. More phylogeny-based studies on connections between niche diversification and speciation are needed across many insect taxa, especially in groups that exhibit few host shifts in relation to speciation.

Background: Leuciscinae is a subfamily belonging to the Cyprinidae fish family that is widely distributed in Circum-Mediterranean region. Many efforts have been carried out to deciphering the evolutionary history of this group. Thus, different biogeographical scenarios have tried to explain the colonization of Europe and Mediterranean area by cyprinids, such as the "north dispersal" or the "Lago Mare dispersal" models. Most recently, Pleistocene glaciations influenced the distribution of leuciscins, especially in North and Central Europe. Weighing up these biogeographical scenarios, this paper constitutes not only the first attempt at deciphering the mitochondrial and nuclear relationships of Mediterranean leuciscins but also a test of biogeographical hypotheses that could have determined the current distribution of Circum-Mediterranean leuciscins. Results: A total of 4439 characters (mitochondrial + nuclear) from 321 individuals of 176 leuciscine species rendered a well-supported phylogeny, showing fourteen main lineages. Analyses of independent mitochondrial and nuclear markers supported the same main lineages, but basal relationships were not concordant. Moreover, some incongruence was found among independent mitochondrial and nuclear phylogenies. The monophyly of some poorly known genera such as Pseudophoxinus and Petroleuciscus was rejected. Representatives of both genera belong to different evolutionary lineages. Timing of cladogenetic events among the main leuciscine lineages was gained using mitochondrial and all genes data set. Conclusions: Adaptations to a predatory lifestyle or miniaturization have superimposed the morphology of some species. These species have been separated into different genera, which are not supported by a phylogenetic framework. Such is the case of the genera Pseudophoxinus and Petroleuciscus, which real taxonomy is not well known. The diversification of leuciscine lineages has been determined by intense vicariant events following the paleoclimatological and hydrogeological history of Mediterranean region. We propose different colonization models of Mediterranean region during the early Oligocene. Later vicariance events promoted Leuciscinae diversification during Oligocene and Miocene periods. Our data corroborate the presence of leuciscins in North Africa before the Messinian salinity crisis. Indeed, Messinian period appears as a stage of gradually Leuciscinae diversification. The rise of humidity at the beginning of the Pliocene promoted the colonization and posterior isolation of newly established freshwater populations. Finally, Pleistocene glaciations determined the current European distribution of some leuciscine species.

Background: More than 2 billion individuals worldwide suffer from helminth infections. The highest parasite burdens occur in children and helminth infection during pregnancy is a risk factor for preterm delivery and reduced birth weight. Therefore, helminth infections can be regarded as a strong selective pressure. Results: Here we propose that candidate susceptibility genes for parasitic worm infections can be identified by searching for SNPs that display a strong correlation with the diversity of helminth species/genera transmitted in different geographic areas. By a genome-wide search we identified 3478 variants that correlate with helminth diversity. These SNPs map to 810 distinct human genes including loci involved in regulatory T cell function and in macrophage activation, as well as leukocyte integrins and co-inhibitory molecules. Analysis of functional relationships among these genes identified complex interaction networks centred around Th2 cytokines. Finally, several genes carrying candidate targets for helminth-driven selective pressure also harbour susceptibility alleles for asthma/allergy or are involved in airway hyper-responsiveness, therefore expanding the known parallelism between these conditions and parasitic infections. Conclusions: Our data provide a landscape of human genes that modulate susceptibility to helminths and indicate parasitic worms as one of the major selective forces in humans.

Background: Among bacteria and archaea, amino acid usage is correlated with habitat temperatures. In particular, protein surfaces in species thriving at higher temperatures appear to be enriched in amino acids that stabilize protein structure and depleted in amino acids that decrease thermostability. Does this observation reflect a causal relationship, or could the apparent trend be caused by phylogenetic relatedness among sampled organisms living at different temperatures? And do proteins from endothermic and exothermic vertebrates show similar differences? Results: We find that the observed correlations between the frequencies of individual amino acids and prokaryotic habitat temperature are strongly influenced by evolutionary relatedness between the species analysed; however, a proteome-wide bias towards increased thermostability remains after controlling for phylogeny. Do eukaryotes show similar effects of thermal adaptation? A small shift of amino acid usage in the expected direction is observed in endothermic ('warm-blooded') mammals and chicken compared to ectothermic ('cold-blooded') vertebrates with lower body temperatures; this shift is not simply explained by nucleotide usage biases. Conclusion: Protein homologs operating at different temperatures have different amino acid composition, both in prokaryotes and in vertebrates. Thus, during the transition from ectothermic to endothermic life styles, the ancestors of mammals and of birds may have experienced weak genome-wide positive selection to increase the thermostability of their proteins.

Background: Harpalinae is a species rich clade of carabid beetles with many unusual morphological forms and ecological interactions. How this diversity evolved has been difficult to reconstruct, perhaps because harpalines underwent a rapid burst of diversification early in their evolutionary history. Here we investigate the tempo of evolution in harpalines using molecular divergence dating techniques and explore the rates of lineage accumulation in harpalines and their sister group. Results: According to molecular divergence date estimates, harpalines originated in the mid Cretaceous but did not diversify extensively until the late Cretaceous or early Paleogene about 32 years after their origin. In a relatively small window of time, harpalines underwent rapid speciation. Harpalines have a relative high net diversification rate and increased cladogenesis in some regions of the clade. We did not see a significant decrease in diversification rate through time in the MCCR test, but a model of diversification with two shift points to lower diversification rates fit the harpaline lineage accumulation through time the best. Conclusions: Our results indicate harpalines are significantly more diverse and have higher diversification than their sistergroup. Instead of an immediate burst of explosive diversification, harpalines may have had a long "fuse" before major lineages diversified during the early Paleogene when other taxa such as mammals, birds, and some flowering plants were also rapidly diversifying.

Background: Female condition-dependent variation in mate preference may have important evolutionary implications, not only within the same population but also among populations. There are few experiments, however, on how condition and/or genotype influences female mate preferences. The black throat patch of the male house sparrow, Passer domesticus, is an intensively studied plumage trait. It is often referred to as a 'badge of status' and seems to be involved in female mate choice, but differences exist among populations. Between-population variation in mate preference may occur for condition-dependent mate preferences. We tested the hypothesis that female preference may vary with female quality (body condition). Therefore, we measured female preference for badge size using an aviary two-choice test in which females were presented with two males that had different sizes of badges (enlarged or averaged). Results: Overall we did not find a female preference for enlarged or average badges, but low-quality females spent more time near average badge males. Conversely, high-quality females did not show a clear preference. Conclusions: Collectively, these results indicate that female preference varies with female quality. Differences in female condition are causes of within-population variation in mating preferences. To our knowledge, our results provide one of the first experimental evidences that variation in preference for a male ornament is associated with female condition. In our study, however, only females of low condition displayed a clear mate preference. Differences observed among populations could be partly explained by differences in female condition.

Background: Hatching enzyme, belonging to the astacin metallo-protease family, digests egg envelope at embryo hatching. Orthologous genes of the enzyme are found in all vertebrate genomes. Recently, we found that exon-intron structures of the genes were conserved among tetrapods, while the genes of teleosts frequently lost their introns. Occurrence of such intron losses in teleostean hatching enzyme genes is an uncommon evolutionary event, as most eukaryotic genes are generally known to be interrupted by introns and the intron insertion sites are conserved from species to species. Here, we report on extensive studies of the exon-intron structures of teleostean hatching enzyme genes for insight into how and why introns were lost during evolution. Results: We investigated the evolutionary pathway of intron-losses in hatching enzyme genes of 27 species of Teleostei. Hatching enzyme genes of basal teleosts are of only one type, which conserves the 9-exon-8-intron structure of an assumed ancestor. On the other hand, otocephalans and euteleosts possess two types of hatching enzyme genes, suggesting a gene duplication event in the common ancestor of otocephalans and euteleosts. The duplicated genes were classified into two clades, clades I and II, based on phylogenetic analysis. In otocephalans and euteleosts, clade I genes developed a phylogeny-specific structure, such as an 8-exon-7-intron, 5-exon-4-intron, 4-exon-3-intron or intron-less structure. In contrast to the clade I genes, the structures of clade II genes were relatively stable in their configuration, and were similar to that of the ancestral genes. Expression analyses revealed that hatching enzyme genes were high-expression genes, when compared to that of housekeeping genes. When expression levels were compared between clade I and II genes, clade I genes tends to be expressed more highly than clade II genes. Conclusions: Hatching enzyme genes evolved to lose their introns, and the intron-loss events occurred at the specific points of teleostean phylogeny. We propose that the high-expression hatching enzyme genes frequently lost their introns during the evolution of teleosts, while the low-expression genes maintained the exon-intron structure of the ancestral gene.

Background: The phenotype of an organism is an outcome of both its genotype, encoding the primary sequence of proteins, and the developmental orchestration of gene expression. The substrate of gene expression in eukaryotes is the chromatin, whose fundamental units are nucleosomes composed of DNA wrapped around each two of the core histone types H2A, H2B, H3 and H4. Key regulatory steps involved in the determination of chromatin conformations are posttranslational modifications (PTM) at histone tails as well as the assembly of histone variants into nucleosomal arrays. Although the mechanistic background is fragmentary understood, it appears that the chromatin signature of metazoan cell types is inheritable over generations. Even less understood is the conservation of epigenetic mechanisms among eukaryotes and their origins. Results: In the light of recent progress in understanding the tree of eukaryotic life we discovered the origin of histone H3 by phylogenetic analyses of variants from all supergroups, which allowed the reconstruction of ancestral states. We found that H3 variants evolved frequently but independently within related species of almost all eukaryotic supergroups. Interestingly, we found all core histone types encoded in the genome of a basal dinoflagellate and H3 variants in two other species, although is was reported that dinoflagellate chromatin is not organized into nucleosomes.Most probably one or more animal/nuclearid H3.3-like variants gave rise to H3 variants of all opisthokonts (animals, choanozoa, fungi, nuclearids, Amoebozoa). H3.2 and H3.1 as well as H3.1t are derivatives of H3.3, whereas H3.2 evolved already in early branching animals, such as Trichoplax. H3.1 and H3.1t are probably restricted to mammals.We deduced a model for protoH3 of the last eukaryotic common ancestor (LECA) confirming a remarkable degree of sequence conservation in comparison to canonical human H3.1. We found evidence that multiple PTMs are conserved even in putatively early branching eukaryotic taxa (Euglenozoa/Excavata). Conclusions: At least a basal repertoire of chromatin modifying mechanisms appears to share old common ancestry and may thus be inherent to all eukaryotes. We speculate that epigenetic principles responsive to environmental triggers may have had influenced phenotypic variation and concomitantly may potentially have had impact on eukaryotic diversification.

Background: Phasmids show noteworthy abilities to overcome species-specific reproductive isolation mechanisms, including hybridization, polyploidy, parthenogenesis, hybridogenesis and androgenesis. From an evolutionary standpoint, such tangled reproductive interactions lead to the complex phyletic relationships known as "reticulate evolution". Moroccan stick insects of the genus Clonopsis include one bisexual (C. felicitatis) and two closely related parthenogenetic forms (C. gallica, C. soumiae), which represent a polyploid series in chromosome number, but with apparent diploid karyotypes. Moreover, two Clonopsis strains of ameiotic males have been described, C. androgenes-35 and C. androgenes-53. As a consequence, Clonopsis stick insects may have experienced complex micro-evolutionary events, which we try to disentangle in this study. Results: Mitochondrial cox2 analysis supports a recent divergence of Clonopsis, while AFLPs evidence genetic differentiation not linked to karyotypes, so that parthenogenetic C. gallica and C. soumiae appear to be a mix of strains of polyphyletic origin rather than single parthenogenetic species. Moreover, an admixed hybrid origin seems to be confirmed for C. androgenes. Conclusion: On the whole, Clonopsis is an intriguing case of reticulate evolution. Actually, complex cladogenetic events should be taken into account to explain the observed genetic structure, including diploidization of polyploid karyotypes, possibly coupled with hybridization and androgenesis. We also proposed a "working hypothesis" to account for the observed data, which deserves further studies, but fits the observed data very well.

Background: The tree of life is significantly asymmetrical - a result of differential speciation and extinction - but general causes of such asymmetry are unclear. Differences in niche partitioning are thought to be one possible general explanation. Ecological specialization might lead to increases in diversification rate or, alternatively, specialization might limit the evolutionary potential of specialist lineages and increase their extinction risk. Here we compare the diversification rates of gall-inducing and non-galling insect lineages. Compared with other insect herbivores feeding on the same host plant, gall-inducing insects feed on plant tissue that is more nutritious and less defended, and they do so in a favorable microhabitat that also provides some protection from natural enemies. We use sister-taxon comparisons to test whether gall-inducing lineages are more host-specific than non-galling lineages, and more or less diverse than non-gallers. We evaluate the significance of diversity bipartitions under equal rates Markov models, and use maximum likelihood model-fitting to test for shifts in diversification rates. Results: We find that, although gall-inducing insect groups are more host-specific than their non-galling relatives, there is no general significant increase in diversification rate in gallers. However, gallers are found at both extremes - two gall-inducing lineages are exceptionally diverse (Euurina sawflies on Salicaceae and Apiomorpha scale insects on Eucalytpus), and one gall-inducing lineage is exceptionally species-poor (Maskellia armored scales on Eucalyptus). Conclusions: The effect of ecological specialization on diversification rates is complex in the case of gall-inducing insects, but host-breadth may be an important factor. When a gall-inducing lineage has a host range approximate to that of its non-galling sister, the gallers are more diverse. When the non-galler clade has a much wider host range than the galler, the non-galler is also much more diverse. There are also lineage-specific effects, with gallers on the same host group exhibiting very different diversities. No single general model explains the observed pattern.

Background: Local adaptation to divergent environmental conditions can promote population genetic differentiation even in the absence of geographic barriers and hence, lead to speciation. Perturbations by catastrophic events, however, can distort such parapatric ecological speciation processes. Here, we asked whether an exceptionally strong flood led to homogenization of gene pools among locally adapted populations of the Atlantic molly (Poecilia mexicana, Poeciliidae) in the Cueva del Azufre system in southern Mexico, where two strong environmental selection factors (darkness within caves and/or presence of toxic H2S in sulfidic springs) drive the diversification of P. mexicana. Nine nuclear microsatellites as well as heritable female life history traits (both as a proxy for quantitative genetics and for trait divergence) were used as markers to compare genetic differentiation, genetic diversity, and especially population mixing (immigration and emigration) before and after the flood. Results: Habitat type (i.e., non-sulfidic surface, sulfidic surface, or sulfidic cave), but not geographic distance was the major predictor of genetic differentiation. Before and after the flood, each habitat type harbored a genetically distinct population. Only a weak signal of individual dislocation among ecologically divergent habitat types was uncovered (with the exception of slightly increased dislocation from the Cueva del Azufre into the sulfidic creek, El Azufre). By contrast, several lines of evidence are indicative of increased flood-induced dislocation within the same habitat type, e.g., between different cave chambers of the Cueva del Azufre. Conclusions: The virtual absence of individual dislocation among ecologically different habitat types indicates strong natural selection against migrants. Thus, our current study exemplifies that ecological speciation in this and other systems, in which extreme environmental factors drive speciation, may be little affected by temporary perturbations, as adaptations to physico-chemical stressors may directly affect the survival probability in divergent habitat types.

Background: The composition of the arthropod head is one of the most contentious issues in animal evolution. In particular, controversy surrounds the homology and innervation of segmental cephalic appendages by the brain. Onychophora (velvet worms) play a crucial role in understanding the evolution of the arthropod brain, because they are close relatives of arthropods and have apparently changed little since the Early Cambrian. However, the segmental origins of their brain neuropils and the number of cephalic appendages innervated by the brain - key issues in clarifying brain composition in the last common ancestor of Onychophora and Arthropoda - remain unclear. Results: Using immunolabelling and neuronal tracing techniques in the developing and adult onychophoran brain, we found that the major brain neuropils arise from only the anterior-most body segment, and that two pairs of segmental appendages are innervated by the brain. The region of the central nervous system corresponding to the arthropod tritocerebrum is not differentiated as part of the onychophoran brain but instead belongs to the ventral nerve cords. Conclusions: Our results contradict the assumptions of a tripartite (three-segmented) brain in Onychophora and instead confirm the hypothesis of bipartite (two-segmented) brain composition. They suggest that the last common ancestor of Onychophora and Arthropoda possessed a brain consisting of protocerebrum and deutocerebrum whereas the tritocerebrum evolved in arthropods.

Background: A phylogenetic network is a generalization of phylogenetic trees that allows the representation of conflicting signals or alternative evolutionary histories in a single diagram. There are several methods for constructing these networks. Some of these methods are based on distances among taxa. In practice, themethods which are based on distance perform faster in comparison with other methods. The Neighbor-Net(N-Net) is a distance-based method. The N-Net produces a circular ordering from a distance matrix, thenconstructs a collection of weighted splits using circular ordering. The SplitsTree which is a program using theseweighted splits makes a phylogenetic network. In general, finding an optimal circular ordering is an NP-hardproblem. The N-Net is a heuristic algorithm to find the optimal circular ordering which is based onneighbor-joining algorithm. Results: In this paper, we present a heuristic algorithm to find an optimal circular ordering based on theMonte-Carlo method, called MC-Net algorithm. In order to show that MC-Net performs better than N-Net, weapply both algorithms on different data sets. Then we draw phylogenetic networks corresponding to outputs ofthese algorithms using SplitsTree and compare the results. Conclusions: We find that the circular ordering produced by the MC-Net is closer to optimal circular orderingthan the N-Net. Furthermore, the networks corresponding to outputs of MC-Net made by SplitsTree are simplerthan N-Net.

Background: The influenza A virus is an important infectious cause of morbidity and mortality in humans and was responsible for 3 pandemics in the 20th century. As the replication of the influenza virus is based on its host's machinery, codon usage of its viral genes might be subject to host selection pressures, especially after interspecies transmission. A better understanding of viral evolution and host adaptive responses might help control this disease. Results: Relative Synonymous Codon Usage (RSCU) values of the genes from segment 1 to segment 6 of avian and human influenza viruses, including pandemic H1N1, were studied via Correspondence Analysis (CA). The codon usage patterns of seasonal human influenza viruses were distinct among their subtypes and different from those of avian viruses. Newly isolated viruses could be added to the CA results, creating a tool to investigate the host origin and evolution of viral genes. It was found that the 1918 pandemic H1N1 virus contained genes with mammalian-like viral codon usage patterns, indicating that the introduction of this virus to humans was not through in toto transfer of an avian influenza virus.Many human viral genes had directional changes in codon usage over time of viral isolation, indicating the effect of host selection pressures. These changes reduced the overall GC content and the usage of G at the third codon position in the viral genome. Limited evidence of translational selection pressure was found in a few viral genes. Conclusions: Codon usage patterns from CA allowed identification of host origin and evolutionary trends in influenza viruses, providing an alternative method and a tool to understand the evolution of influenza viruses. Human influenza viruses are subject to selection pressure on codon usage which might assist in understanding the characteristics of newly emerging viruses.

Background: Many butterfly species have been experiencing the northward range expansion and physiological adaptation, probably due to climate warming. Here, we document an extraordinary field case of a species of lycaenid butterfly, Zizeeria maha, for which plastic phenotypes of wing color-patterns were revealed at the population level in the course of range expansion. Furthermore, we examined whether this outbreak of phenotypic changes was able to be reproduced in a laboratory. Results: In the recently expanded northern range margins of this species, more than 10% of the Z. maha population exhibited characteristic color-pattern modifications on the ventral wings for three years. We physiologically reproduced similar phenotypes by an artificial cold-shock treatment of a normal southern population, and furthermore, we genetically reproduced a similar phenotype after selective breeding of a normal population for ten generations, demonstrating that the cold-shock-induced phenotype was heritable and partially assimilated genetically in the breeding line. Similar genetic process might have occurred in the previous and recent range-margin populations as well. Relatively minor modifications expressed in the tenth generation of the breeding line together with other data suggest a role of founder effect in this field case. Conclusions: Our results support the notion that the outbreak of the modified phenotypes in the recent range-margin population was primed by the revelation of plastic phenotypes in response to temperature stress and by the subsequent genetic process in the previous range-margin population, followed by migration and temporal establishment of genetically unstable founders in the recent range margins. This case presents not only an evolutionary role of phenotypic plasticity in the field but also a novel evolutionary aspect of range expansion at the species level.

Background: Aphids are striking in their prodigious reproductive capacity and reliance on microbial endosymbionts, which provision their hosts with necessary amino acids and provide protection against parasites and heat stress. Perhaps as a result of this bacterial dependence, aphids have limited immune function that may leave them vulnerable to bacterial pathogens. An alternative, non-immunological response that may be available to infected aphids is to increase reproduction, thereby ameliorating fitness loss from infection. Such a response would reduce the need to mount a potentially energetically costly immune response, and would parallel that of other hosts that alter life-history traits when there is a risk of infection. Here we examined whether pea aphids (Acyrthosiphon pisum) respond to immunological challenges by increasing reproduction. As a comparison to the response to the internal cue of risk elicited by immunological challenge, we also exposed pea aphids to an external cue of risk - the aphid alarm pheromone (E)-beta-farnesene (EBF), which is released in the presence of predators. For each challenge, we also examined whether the presence of symbionts modified the host response, as maintaining host fitness in the face of challenge would benefit both the host and its dependent bacteria. Results: We found that aphids stabbed abdominally with a sterile needle had reduced fecundity relative to control aphids but that aphids stabbed with a needle bearing heat-killed bacteria had reproduction intermediate, and statistically indistinguishable, to the aphids stabbed with a sterile needle and the controls. Aphids with different species of facultative symbiotic bacteria had different reproductive patterns overall, but symbionts in general did not alter aphid reproduction in response to bacterial exposure. However, in response to exposure to alarm pheromone, aphids with Hamiltonella defensa or Serratia symbiotica symbiotic infections increased reproduction but those without a facultative symbiont or with Regiella insecticola did not. Conclusions: Overall, our results suggest that pea aphids are able to increase their reproduction in response to specific cues and that symbiont presence sometimes moderates this response. Such increased reproduction in response to risk of death increases the fitness of both aphids and their vertically transmitted symbionts, and since these organisms have high reproductive capacity, slight increases in reproduction could lead to a very large numerical advantage later in the season. Thus both symbiotic partners can benefit by increasing host fecundity under dangerous conditions.