Can 15-Locus Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeat Analysis Provide Insight into the Evolution of Mycobacterium tuberculosis?

The phylogeny and evolution of the bacterium Mycobacterium tuberculosis is still poorly understood despite the application of a variety of molecular techniques. We analyzed 469 M. tuberculosis and 49 Mycobacterium bovis isolates to evaluate if the mycobacterial interspersed repetitive units-variable-number tandem repeats (MIRU-VNTR) commonly used for epidemiological studies can define the phylogeny of the M. tuberculosis complex. This population was characterized by previously identified silent single-nucleotide polymorphisms (sSNPs) or by a macroarray based on these sSNPs that was developed in this study. MIRU-VNTR phylogenetic codes capable of differentiating between phylogenetic lineages were identified. Overall, there was 90.9% concordance between the lineages of isolates as defined by the MIRU-VNTR and sSNP analyses. The MIRU-VNTR phylogenetic code was unique to M. bovis and was not observed in any M. tuberculosis isolates. The codes were able to differentiate between different M. tuberculosis strain families such as Beijing, Delhi, and East African-Indian. Discrepant isolates with similar but not identical MIRU-VNTR codes often displayed a stepwise trend suggestive of bidirectional evolution. A lineage-specific panel of MIRU-VNTR can be used to subdivide each lineage for epidemiological purposes. MIRU-VNTR is a valuable tool for phylogenetic studies and could define an evolutionarily uncharacterized population of M. tuberculosis complex organisms.     

Phylogenetic incongruence in Cymbidium orchids

Cymbidium, which includes approximately 80 species, is one of the most ornamental and cultivated orchid genera. However, a lack of markers and sparse sampling have posed great challenges to resolving the phylogenetic relationships within the genus. In the present study, we reconstructed the phylogenetic relationships by utilizing one nuclear DNA (nrITS) and seven plastid genes (rbcL, trnS, trnG, matK, trnL, psbA, and atpI) from 70 species (varieties) in Cymbidium. We also examined the occurrence of phylogenetic conflict between nuclear (nrITS) and plastid loci and investigated how phylogenetic conflict bears on taxonomic classification within the genus. We found that phylogenetic conflict and low support values may be explained by hybridization and a lack of informative characteristics. Our results do not support previous classification of the subgenera and sections within Cymbidium. Discordance between gene trees and network analysis indicate that reticulate evolution occurred in the genus Cymbidium. Overall, our study indicates that Cymbidium has undergone a complex evolution.     

Phylogeny Reconstruction and Hybrid Analysis of Populus (Salicaceae) Based on Nucleotide Sequences of Multiple Single-Copy Nuclear Genes and Plastid Fragments

Populus (Salicaceae) is one of the most economically and ecologically important genera of forest trees. The complex reticulate evolution and lack of highly variable orthologous single-copy DNA markers have posed difficulties in resolving the phylogeny of this genus. Based on a large data set of nuclear and plastid DNA sequences, we reconstructed robust phylogeny of Populus using parsimony, maximum likelihood and Bayesian inference methods. The resulting phylogenetic trees showed better resolution at both inter- and intra-sectional level than previous studies. The results revealed that (1) the plastid-based phylogenetic tree resulted in two main clades, suggesting an early divergence of the maternal progenitors of Populus; (2) three advanced sections (Populus, Aigeiros and Tacamahaca) are of hybrid origin; (3) species of the section Tacamahaca could be divided into two major groups based on plastid and nuclear DNA data, suggesting a polyphyletic nature of the section; and (4) many species proved to be of hybrid origin based on the incongruence between plastid and nuclear DNA trees. Reticulate evolution may have played a significant role in the evolution history of Populus by facilitating rapid adaptive radiations into different environments.     

Large-scale patterns of signal evolution: an interspecific study of Liolaemus lizard headbob displays

Although many factors have been shown to influence the evolution of species recognition signals in a wide variety of taxa, it is difficult to draw general conclusions because of fundamental differences in the morphologies and ecologies of the animals considered. In this study, two morphologically and ecologically similar lizard genera (Sceloporus and Liolaemus) are used to provide replicate examples of the evolution of a complex visual display. New data on the headbob displays of 16 Liolaemus species are presented. As in other taxa, phylogenetic analyses show that evolutionary changes in display structure have been rapid, leaving little, if any, phylogenetic information in the display structure. Evolutionary changes in display structure also do not appear to be closely associated with any major habitat characteristics. Despite this rapid evolution, Liolaemus lizards produce headbob displays that are remarkably simple in structure in comparison to those produced by Sceloporus, perhaps compensating for lower complexity by frequent use of other visual displays such as forelimb and tail waves.     

Species replacement along a linear coastal habitat: phylogeography and speciation in the red alga Mazzaella laminarioides along the south east pacific

Background

Conserved domains are recognized as the building blocks of eukaryotic proteins. Domains showing a tendency to occur in diverse combinations (??promiscuous?? domains) are involved in versatile architectures in proteins with different functions. Current models, based on global-level analyses of domain combinations in multiple genomes, have suggested that the propensity of some domains to associate with other domains in high-level architectures increases with organismal complexity. Alternative models using domain-based phylogenetic trees propose that domains have become promiscuous independently in different lineages through convergent evolution and are, thus, random with no functional or structural preferences. Here we test whether complex protein architectures have occurred by accretion from simpler systems and whether the appearance of multidomain combinations parallels organismal complexity. As a model, we analyze the modular evolution of the PWWP domain and ask whether its appearance in combinations with other domains into multidomain architectures is linked with the occurrence of more complex life-forms. Whether high-level combinations of domains are conserved and transmitted as stable units (cassettes) through evolution is examined in the genomes of plant or metazoan species selected for their established position in the evolution of the respective lineages.

Results

Using the domain-tree approach, we analyze the evolutionary origins and distribution patterns of the promiscuous PWWP domain to understand the principles of its modular evolution and its existence in combination with other domains in higher-level protein architectures. We found that as a single module the PWWP domain occurs only in proteins with a limited, mainly, species-specific distribution. Earlier, it was suggested that domain promiscuity is a fast-changing (volatile) feature shaped by natural selection and that only a few domains retain their promiscuity status throughout evolution. In contrast, our data show that most of the multidomain PWWP combinations in extant multicellular organisms (humans or land plants) are present in their unicellular ancestral relatives suggesting they have been transmitted through evolution as conserved linear arrangements (??cassettes??). Among the most interesting biologically relevant results is the finding that the genes of the two plant Trithorax family subgroups (ATX1/2 and ATX3/4/5) have different phylogenetic origins. The two subgroups occur together in the earliest land plants Physcomitrella patens and Selaginella moellendorffii.

Conclusion

Gain/loss of a single PWWP domain is observed throughout evolution reflecting dynamic lineage- or species-specific events. In contrast, higher-level protein architectures involving the PWWP domain have survived as stable arrangements driven by evolutionary descent. The association of PWWP domains with the DNA methyltransferases in O. tauri and in the metazoan lineage seems to have occurred independently consistent with convergent evolution. Our results do not support models wherein more complex protein architectures involving the PWWP domain occur with the appearance of more evolutionarily advanced life forms.     

A phylogenomic approach to microbial evolution

To study the origin and evolution of biochemical pathways in microorganisms, we have developed methods and software for automatic, large-scale reconstructions of phylogenetic relationships. We define the complete set of phylogenetic trees derived from the proteome of an organism as the phylome and introduce the term phylogenetic connection as a concept that describes the relative relationships between taxa in a tree. A query system has been incorporated into the system so as to allow searches for defined categories of trees within the phylome. As a complement, we have developed the pyphy system for visualising the results of complex queries on phylogenetic connections, genomic locations and functional assignments in a graphical format. Our phylogenomics approach, which links phylogenetic information to the flow of biochemical pathways within and among microbial species, has been used to examine more than 8000 phylogenetic trees from seven microbial genomes. The results have revealed a rich web of phylogenetic connections. However, the separation of Bacteria and Archaea into two separate domains remains robust.     

Modular Evolution of the Drosophila Metabolome

Comparative phylogenetic studies offer a powerful approach to study the evolution of complex traits. Although much effort has been devoted to the evolution of the genome and to organismal phenotypes, until now relatively little work has been done on the evolution of the metabolome, despite the fact that it is composed of the basic structural and functional building blocks of all organisms. Here we explore variation in metabolite levels across 50 My of evolution in the genus Drosophila, employing a common garden design to measure the metabolome within and among 11 species of Drosophila. We find that both sex and age have dramatic and evolutionarily conserved effects on the metabolome. We also find substantial evidence that many metabolite pairs covary after phylogenetic correction, and that such metabolome coevolution is modular. Some of these modules are enriched for specific biochemical pathways and show different evolutionary trajectories, with some showing signs of stabilizing selection. Both observations suggest that functional relationships may ultimately cause such modularity. These coevolutionary patterns also differ between sexes and are affected by age. We explore the relevance of modular evolution to fitness by associating modules with lifespan variation measured in the same common garden. We find several modules associated with lifespan, particularly in the metabolome of older flies. Oxaloacetate levels in older females appear to coevolve with lifespan, and a lifespan-associated module in older females suggests that metabolic associations could underlie 50 My of lifespan evolution.     

Evolution and phylogenetic relationships in subterranean rodents of the Ctenomys mendocinus species complex: Effects of Late Quaternary landscape changes of Central Argentina

We conducted a thorough survey of populations of Ctenomys within the distributional ranges of all species that currently represent the Ctenomys mendocinus species complex in Central Argentina. From the Atlantic Southeastern coast to the heights of the Los Andes mountain range we sampled individuals for all the previously recognized species within this species complex (C. australis, C. azarae, C. porteousi, C. “chasiquensis” and C. mendocinus). Sequences of different fragments of mitochondrial DNA were used to assess the phylogenetic relationships between individuals and to understand the processes responsible for the observed geographic distribution of genetic variation. The genetic differentiation among most of the sequences was moderated, although highly divergent haplotypes were recorded in some localities in the West of our study area, suggestive of the occurrence of a new so far undescribed species. The phylogenetic reconstructions suggested the presence of three major clades into C. mendocinus complex species that do not agree with the currently proposed taxonomy. Phylogeographic and demographic estimates indicated that the C. mendocinus species complex experienced a strong and recent process of population expansion. We propose that the pattern of distribution of haplotype diversity in the C. mendocinus species complex corresponds to a scenario of mild local differentiation with subsequent expansion to the current distribution ranges. Demographic changes through time suggest that glacial-interglacial cycles of Late Quaternary, and the consequent variation in the aridity of Central plains in Argentina, had strong effects on the evolution and cladogenesis of the C. mendocinus species complex.     

Pythium intermedium,a species complex consisting of three phylogenetic species found in cool-temperate forest ecosystems

Pythium intermedium plays a vital role in the carbon cycle of cool-temperate forests and is widely distributed in Japan's forest soils. In this study, we performed a phylogenetic analysis of the P. intermedium species complex using DNA sequences from multiple loci. The study included 35 isolates from cool-temperate forest soils, seven known P. intermedium isolates, and six known Pythium attrantheridium isolates. We also performed morphological observations and mating tests. Our results showed that all the isolates formed one large clade but were divided into three subclades. Furthermore, we observed many mating reactions between isolates from different subclades, including between P. attrantheridium and P. intermedium. Therefore, we suggest that P. intermedium, P. attrantheridium, and another phylogenetic species belong to one species complex. This is the first report of a species complex within P. intermedium and will be helpful in understanding the evolution of Pythium species in natural ecosystems.     

Origin and Evolution of the Sodium -Pumping NADH: Ubiquinone Oxidoreductase

The sodium -pumping NADH: ubiquinone oxidoreductase (Na⁺-NQR) is the main ion pump and the primary entry site for electrons into the respiratory chain of many different types of pathogenic bacteria. This enzymatic complex creates a transmembrane gradient of sodium that is used by the cell to sustain ionic homeostasis, nutrient transport, ATP synthesis, flagellum rotation and other essential processes. Comparative genomics data demonstrate that the nqr operon, which encodes all Na⁺-NQR subunits, is found in a large variety of bacterial lineages with different habitats and metabolic strategies. Here we studied the distribution, origin and evolution of this enzymatic complex. The molecular phylogenetic analyses and the organizations of the nqr operon indicate that Na⁺-NQR evolved within the Chlorobi/Bacteroidetes group, after the duplication and subsequent neofunctionalization of the operon that encodes the homolog RNF complex. Subsequently, the nqr operon dispersed through multiple horizontal transfer events to other bacterial lineages such as Chlamydiae, Planctomyces and α, β, γ and δ -proteobacteria. Considering the biochemical properties of the Na⁺-NQR complex and its physiological role in different bacteria, we propose a detailed scenario to explain the molecular mechanisms that gave rise to its novel redox- dependent sodium -pumping activity. Our model postulates that the evolution of the Na⁺-NQR complex involved a functional divergence from its RNF homolog, following the duplication of the rnf operon, the loss of the rnfB gene and the recruitment of the reductase subunit of an aromatic monooxygenase.     

Episodes of adaptive evolution of mitochondrial genome in asiatic salamanders (Amphibia,Caudata, Hynobiidae)

To elucidate the effect of natural selection on the evolution of mitochondrial DNA (mtDNA) in Asiatic salamanders of the family Hynobiidae, nucleotide sequences of 12 protein-coding genes were analyzed. Using a mixed effects model of evolution, it was found that, in spite of the pronounced effect of negative selection on the mtDNA evolution in Hynobiidae (which is typical for the animals in general), two phylogenetic clusters, the West Asian one, represented by the genera Ranodon and Paradactylodon, and North Eurasian one, represented by the genus Salamandrella, were formed under the influence of episodic positive selection. Analysis of protein sequences encoded by the mitochondrial genome also supported the influence of positive selection on the evolution of Hynobiidae at some stages of their cladogenesis. It is suggested that the signatures of adaptive evolution detected in the mtDNA of Hynobiidae were determined by the complex and long-lasting history of their formation, accompanied by adaptation to the changing environment.     

The role of pollinator diversity in the evolution of corolla-shape integration in a pollination-generalist plant clade

Flowers of animal-pollinated plants are integrated structures shaped by the action of pollinator-mediated selection. It is widely assumed that pollination specialization increases the magnitude of floral integration. However, empirical evidence is still inconclusive. In this study, we explored the role of pollinator diversity in shaping the evolution of corolla-shape integration in Erysimum, a plant genus with generalized pollination systems. We quantified floral integration in Erysimum using geometric morphometrics and explored its evolution using phylogenetic comparative methods. Corolla-shape integration was low but significantly different from zero in all study species. Spatial autocorrelation and phylogenetic signal in corolla-shape integration were not detected. In addition, integration in Erysimum seems to have evolved in a way that is consistent with Brownian motion, but with frequent convergent evolution. Corolla-shape integration was negatively associated with the number of pollinators visiting the flowers of each Erysimum species. That is, it was lower in those species having a more generalized pollination system. This negative association may occur because the co-occurrence of many pollinators imposes conflicting selection and cancels out any consistent selection on specific floral traits, preventing the evolution of highly integrated flowers.     

Environmental correlates of phenotypic evolution in ecologically diverse Liolaemus lizards

Evolutionary correlations between phenotypic and environmental traits characterize adaptive radiations. However, the lizard genus Liolaemus, one of the most ecologically diverse terrestrial vertebrate radiations on earth, has so far shown limited or mixed evidence of adaptive diversification in phenotype. Restricted use of comprehensive environmental data, incomplete taxonomic representation and not considering phylogenetic uncertainty may have led to contradictory evidence. We compiled a 26‐taxon dataset for the Liolaemus gracilis species group, representing much of the ecological diversity represented within Liolaemus and used environmental data to characterize how environments occupied by species'' relate to phenotypic evolution. Our analyses, explicitly accounting for phylogenetic uncertainty, suggest diversification in phenotypic traits toward the present, with body shape evolution rapidly evolving in this group. Body shape evolution correlates with the occupation of different structural habitats indicated by vegetation axes suggesting species have adapted for maximal locomotory performance in these habitats. Our results also imply that the effects of phylogenetic uncertainty and model misspecification may be more extensive on univariate, relative to multivariate analyses of evolutionary correlations, which is an important consideration in analyzing data from rapidly radiating adaptive radiations.     

Accommodating natural and sexual selection in butterfly wing pattern evolution

Visual patterns in animals may serve different functions, such as attracting mates and deceiving predators. If a signal is used for multiple functions, the opportunity arises for conflict among the different functions, preventing optimization for any one visual signal. Here we investigate the hypothesis that spatial separation of different visual signal functions has occurred in Bicyclus butterflies. Using phylogenetic reconstructions of character evolution and comparisons of evolutionary rates, we found dorsal surface characters to evolve at higher rates than ventral characters. Dorsal characters also displayed sex-based differences in evolutionary rates more often than did ventral characters. Thus, dorsal characters corresponded to our predictions of mate signalling while ventral characters appear to play an important role in predator avoidance. Forewing characters also fit a model of mate signalling, and displayed higher rates of evolution than hindwing characters. Our results, as well as the behavioural and developmental data from previous studies of Bicyclus species, support the hypothesis that spatial separation of visual signal functions has occurred in Bicyclus butterflies. This study is the first to demonstrate, in a phylogenetic framework, that spatial separation of signals used for mate signalling and those used for predator avoidance is a viable strategy to accommodate multiple signal functions. This signalling strategy has important ramifications on the developmental evolution of wing pattern elements and diversification of butterfly species.     

Evolution of the tRNALeu (UAA) Intron and Congruence of Genetic Markers in Lichen-Symbiotic Nostoc

The group I intron interrupting the tRNALeu UAA gene (trnL) is present in most cyanobacterial genomes as well as in the plastids of many eukaryotic algae and all green plants. In lichen symbiotic Nostoc, the P6b stem-loop of trnL intron always involves one of two different repeat motifs, either Class I or Class II, both with unresolved evolutionary histories. Here we attempt to resolve the complex evolution of the two different trnL P6b region types. Our analysis indicates that the Class II repeat motif most likely appeared first and that independent and unidirectional shifts to the Class I motif have since taken place repeatedly. In addition, we compare our results with those obtained with other genetic markers and find strong evidence of recombination in the 16S rRNA gene, a marker widely used in phylogenetic studies on Bacteria. The congruence of the different genetic markers is successfully evaluated with the recently published software Saguaro, which has not previously been utilized in comparable studies.     

BIOLOGICAL AND HISTORICAL FACTORS INFLUENCING GENETIC DIVERSITY IN THE SCUTELLARIA ANGUSTIFOLIA COMPLEX (LABIATAE)

Studies of genetic diversity at isozyme loci were used to examine the phylogenetic distribution of several frequently reported population-genetic parameters in a putatively monophyletic group of plant species, the Scutellaria angustifolia complex. The influence of taxon-specific differences in habitat preference, breeding system, degree of endemism, and phylogenetic relatedness was examined. Many characters of reproductive morphology traditionally used in phylogenetic inference vary with breeding system. To the extent that reproductive systems are conservative markers of phylogenetic relationships, one would expect the distribution of genetic variation to be similar in closely related taxa. Results showed that closely related taxa may exhibit very different genetic-diversity statistics and that distantly related taxa may exhibit very similar genetic-diversity statistics. This suggests that complex patterns of evolution of breeding systems and morphological characters have occurred in the ten taxa included in the Scutellaria angustifolia complex. Similarity in habitat is not associated with similarity in genetic diversity in this group of species.     

Congruence between comparative morphology and molecular phylogenies: Evolution of the male genital skeletal/muscular system in the subtribe Polyommatina (Lepidoptera,Lycaenidae)

The skeleton and musculature of male genitalia were studied in species of a model butterfly group (subtribe Polyommatina, Lycaenidae). In total, we analyzed 45 species of the tribe Polyommatini most of which were previously used in the molecular phylogenetic study (Talavera et al., 2013). The studied morphological characters were mapped on the molecular trees, which allowed us to reveal trends of morphological changes and to estimate the age of their origin. As a result, chronology of evolution of skeleton and musculature traits was established. It was shown that periods of slow morphological evolution alternated in the subtribe Polyommatina with those of a high rate of origin of new traits. For example, topography of the intravalvar muscles has not changed for 26 MY preserving their initial fan-shaped attachment. The evolution of intravalvar muscles started 10 MYA, proceeded slowly during the first 5 MY, and then accelerated during the last 5 MY resulting in the extensive splitting of the musculature in most monophyletic lineages. Mapping the morphological characters on the phylogeny demonstrated that the rates of skeleton and muscle evolution within the skeleton/musculature apparatus were different. In most cases the intravalvar musculature evolved much faster than the skeleton. The cladistic interpretation of states of morphological traits was found to be consistent with phylogenetic reconstructions based on analysis of multiple molecular markers. Moreover, morphological synapomorphies were found for the lineages Alpherakya + Glabroculus and Aricia + (Alpherakya + Glabroculus), which had low statistical support in molecular phylogenetic analysis. Additionally, in some cases molecular studies helped to reveal trends in the evolution of morphological traits. For example, the unpaired uncus and the compact juxta are not plesiomorphic for Cupidina as previously thought; instead, they were shown to have evolved secondarily within this subtribe.     

Pythium insidiosum complex hides a cryptic novel species: Pythium periculosum

Early phylogenetic analysis of Pythium insidiosum, the etiologic agent of pythiosis in mammals, showed the presence of a complex comprising three monophyletic clusters. Two included isolates recovered from cases of pythiosis in the Americas (Cluster I) and Asia (Cluster II), whereas the third cluster included four diverged isolates three from humans in Thailand and the USA, and one isolate from a USA spectacled bear (Cluster III). Thereafter, several phylogenetic analyses confirmed the presence of at least three monophyletic clusters, with most isolates placed in clusters I and II. Recent phylogenetic analyses using isolates from environmental sources and from human cases in India, Spain, Thailand, and dogs in the USA, however, showed the presence of two monophyletic groups each holding two sub-clusters. These studies revealed that P. insidiosum possesses different phylogenetic patterns to that described by early investigators. In this study, phylogenetic, population genetic and protein MALDI-TOF analyses of the P. insidiosum isolates in our culture collection, as well as those available in the database, showed members in the proposed cluster III and IV are phylogenetically different from that in clusters I and II. Our analyses of the complex showed a novel group holding two sub-clusters the USA (Cluster III) and the other from different world regions (Cluster IV). The data showed the original P. insidiosum cluster III is a cryptic novel species, now identified as P. periculosum. The finding of a novel species within P. insidiosum complex has direct implications in the epidemiology, diagnosis, and management of pythiosis in mammalian hosts.     

Genetic Structure and Evolution of the Leishmania Genus in Africa and Eurasia: What Does MLSA Tell Us

Leishmaniasis is a complex parasitic disease from a taxonomic, clinical and epidemiological point of view. The role of genetic exchanges has been questioned for over twenty years and their recent experimental demonstration along with the identification of interspecific hybrids in natura has revived this debate. After arguing that genetic exchanges were exceptional and did not contribute to Leishmania evolution, it is currently proposed that interspecific exchanges could be a major driving force for rapid adaptation to new reservoirs and vectors, expansion into new parasitic cycles and adaptation to new life conditions.To assess the existence of gene flows between species during evolution we used MLSA-based (MultiLocus Sequence Analysis) approach to analyze 222 Leishmania strains from Africa and Eurasia to accurately represent the genetic diversity of this genus. We observed a remarkable congruence of the phylogenetic signal and identified seven genetic clusters that include mainly independent lineages which are accumulating divergences without any sign of recent interspecific recombination. From a taxonomic point of view, the strong genetic structuration of the different species does not question the current classification, except for species that cause visceral forms of leishmaniasis (L. donovani, L. infantum and L. archibaldi). Although these taxa cause specific clinical forms of the disease and are maintained through different parasitic cycles, they are not clearly distinct and form a continuum, in line with the concept of species complex already suggested for this group thirty years ago. These results should have practical consequences concerning the molecular identification of parasites and the subsequent therapeutic management of the disease.