Gene trees and species trees are not the same

Archaeological and genetic insights into the origins of domesticated rice

Rice (Oryza sativa) is one of the most important cereal grains in the world today and serves as a staple food source for more than half of the world's population. Research into when, where, and how rice was brought into cultivation and eventually domesticated, along with its development into a staple food source, is thus essential. These questions have been a point of nearly continuous research in both archaeology and genetics, and new information has continually come to light as theory, data acquisition, and analytical techniques have advanced over time. Here, we review the broad history of our scientific understanding of the rice domestication process from both an archaeological and genetic perspective and examine in detail the information that has come to light in both of these fields in the last 10 y. Current findings from genetics and archaeology are consistent with the domestication of O. sativa japonica in the Yangtze River valley of southern China. Interestingly, although it appears rice was cultivated in the area by as early 8000 BP, the key domestication trait of nonshattering was not fixed for another 1,000 y or perhaps longer. Rice was also cultivated in India as early as 5000 BP, but the domesticated indica subspecies currently appears to be a product of the introgression of favorable alleles from japonica. These findings are reshaping our understanding of rice domestication and also have implications for understanding the complex evolutionary process of plant domestication.

Identification of putative orthologous genes for the phylogenetic reconstruction of temperate woody bamboos (Poaceae: Bambusoideae)

The temperate woody bamboos (Arundinarieae) are highly diverse in morphology but lack a substantial amount of genetic variation. The taxonomy of this lineage is intractable, and the relationships within the tribe have not been well resolved. Recent studies indicated that this tribe could have a complex evolutionary history. Although phylogenetic studies of the tribe have been carried out, most of these phylogenetic reconstructions were based on plastid data, which provide lower phylogenetic resolution compared with nuclear data. In this study, we intended to identify a set of desirable nuclear genes for resolving the phylogeny of the temperate woody bamboos. Using two different methodologies, we identified 209 and 916 genes, respectively, as putative single copy orthologous genes. A total of 112 genes was successfully amplified and sequenced by next-generation sequencing technologies in five species sampled from the tribe. As most of the genes exhibited intra-individual allele heterozygotes, we investigated phylogenetic utility by reconstructing the phylogeny based on individual genes. Discordance among gene trees was observed and, to resolve the conflict, we performed a range of analyses using BUCKy and HybTree. While caution should be taken when inferring a phylogeny from multiple conflicting genes, our analysis indicated that 74 of the 112 investigated genes are potential markers for resolving the phylogeny of the temperate woody bamboos.

So what is a species anyway? A primatological perspective

Since Darwin's time, the question "what a species" has provoked fierce disputes and a tremendous number of publications, from short opinion papers to thick volumes. The debates covered fundamental philosophical questions, such as: Do species exist at all independently of a human observer or are they just a construct of the human mind to categorize nature's organismic diversity and serve as a semantic tool in human communication about biodiversity? or: Are species natural kinds (classes) or individuals that are "born" by speciation, change in course of time, and finally "die" when they go extinct or diverge into new species? Also included was the problem of species as taxa (taxonomic) versus species as products of the speciation process (evolutionary). More pragmatic issues arose, such as: How can we reliably delineate and delimitate species? The great interest in what a species is reflects the importance of "species" as fundamental units in most fields of biology, especially evolutionary biology, ecology, and conservation.

True lemurs…true species - species delimitation using multiple data sources in the brown lemur complex

Background

Species are the fundamental units in evolutionary biology. However, defining them as evolutionary independent lineages requires integration of several independent sources of information in order to develop robust hypotheses for taxonomic classification. Here, we exemplarily propose an integrative framework for species delimitation in the “brown lemur complex” (BLC) of Madagascar, which consists of seven allopatric populations of the genus Eulemur (Primates: Lemuridae), which were sampled extensively across northern, eastern and western Madagascar to collect fecal samples for DNA extraction as well as recordings of vocalizations. Our data base was extended by including museum specimens with reliable identification and locality information for skull shape and pelage color analysis.

Results

Between-group analyses of principal components revealed significant heterogeneity in skull shape, pelage color variation and loud calls across all seven populations. Furthermore, post-hoc statistical tests between pairs of populations revealed considerable discordance among different data sets for different dyads. Despite a high degree of incomplete lineage sorting among nuclear loci, significant exclusive ancestry was found for all populations, except for E. cinereiceps, based on one mitochondrial and three nuclear genetic loci.

Conclusions

Using several independent lines of evidence, our results confirm the species status of the members of the BLC under the general lineage concept of species. More generally, the present analyses demonstrate the importance and value of integrating different kinds of data in delimiting recently evolved radiations.

Phylogenetic resolution power of microsatellites and various single-nucleotide polymorphism types assessed in 10 divergent chicken populations

There has been some debate over the question of which types of DNA variation are most appropriate to accurately reconstruct evolutionary events. We compared the capacity of microsatellites (STRs) and various types of single-nucleotide polymorphism (SNP) loci in the chicken genome. The SNP types differ in their location: in exons, introns and promoters. Genetic distances between all possible pairs of 10 populations were calculated for each marker type. STR loci, which are much more polymorphic than are SNPs, are considered to have occurred at recent time compared with old evolutionary events of SNPs. Using structure software, STR loci assigned individuals to their population much more correctly than did any other marker types, whereas SNPs at promoter regions gave the poorest ascription. Furthermore, 29 STR markers were even better than all 152 SNPs together. Ancient evolutionary events that produced genetic differences between the most distant populations such as Red Jungle Fowl and domestic chicken were detected better by exons and introns than by STR loci and promoters. The significant interactions found between marker types and populations suggest that marker types had different phylogenetic histories, possibly related to a different timescale.

Contrasted evolutionary histories of two Toll-like receptors (Tlr4 and Tlr7) in wild rodents (MURINAE)

BACKGROUND

In vertebrates, it has been repeatedly demonstrated that genes encoding proteins involved in pathogen-recognition by adaptive immunity (e.g. MHC) are subject to intensive diversifying selection. On the other hand, the role and the type of selection processes shaping the evolution of innate-immunity genes are currently far less clear. In this study we analysed the natural variation and the evolutionary processes acting on two genes involved in the innate-immunity recognition of Microbe-Associated Molecular Patterns (MAMPs).

RESULTS

We sequenced genes encoding Toll-like receptor 4 (Tlr4) and 7 (Tlr7), two of the key bacterial- and viral-sensing receptors of innate immunity, across 23 species within the subfamily Murinae. Although we have shown that the phylogeny of both Tlr genes is largely congruent with the phylogeny of rodents based on a comparably sized non-immune sequence dataset, we also identified several potentially important discrepancies. The sequence analyses revealed that major parts of both Tlrs are evolving under strong purifying selection, likely due to functional constraints. Yet, also several signatures of positive selection have been found in both genes, with more intense signal in the bacterial-sensing Tlr4 than in the viral-sensing Tlr7. 92% and 100% of sites evolving under positive selection in Tlr4 and Tlr7, respectively, were located in the extracellular domain. Directly in the Ligand-Binding Region (LBR) of TLR4 we identified two rapidly evolving amino acid residues and one site under positive selection, all three likely involved in species-specific recognition of lipopolysaccharide of gram-negative bacteria. In contrast, all putative sites of LBRTLR7 involved in the detection of viral nucleic acids were highly conserved across rodents. Interspecific differences in the predicted 3D-structure of the LBR of both Tlrs were not related to phylogenetic history, while analyses of protein charges clearly discriminated Rattini and Murini clades.

CONCLUSIONS

In consequence of the constraints given by the receptor protein function purifying selection has been a dominant force in evolution of Tlrs. Nevertheless, our results show that episodic diversifying parasite-mediated selection has shaped the present species-specific variability in rodent Tlrs. The intensity of diversifying selection was higher in Tlr4 than in Tlr7, presumably due to structural properties of their ligands.

Evidence for inter-specific recombination among the mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex

Background

The availability of mitochondrial genomes has allowed for the resolution of numerous questions regarding the evolutionary history of fungi and other eukaryotes. In the Gibberella fujikuroi species complex, the exact relationships among the so-called “African”, “Asian” and “American” Clades remain largely unresolved, irrespective of the markers employed. In this study, we considered the feasibility of using mitochondrial genes to infer the phylogenetic relationships among Fusarium species in this complex. The mitochondrial genomes of representatives of the three Clades (Fusarium circinatum, F. verticillioides and F. fujikuroi) were characterized and we determined whether or not the mitochondrial genomes of these fungi have value in resolving the higher level evolutionary relationships in the complex.

Results

Overall, the mitochondrial genomes of the three species displayed a high degree of synteny, with all the genes (protein coding genes, unique ORFs, ribosomal RNA and tRNA genes) in identical order and orientation, as well as introns that share similar positions within genes. The intergenic regions and introns generally contributed significantly to the size differences and diversity observed among these genomes. Phylogenetic analysis of the concatenated protein-coding dataset separated members of the Gibberella fujikuroi complex from other Fusarium species and suggested that F. fujikuroi (“Asian” Clade) is basal in the complex. However, individual mitochondrial gene trees were largely incongruent with one another and with the concatenated gene tree, because six distinct phylogenetic trees were recovered from the various single gene datasets.

Conclusion

The mitochondrial genomes of Fusarium species in the Gibberella fujikuroi complex are remarkably similar to those of the previously characterized Fusarium species and Sordariomycetes. Despite apparently representing a single replicative unit, all of the genes encoded on the mitochondrial genomes of these fungi do not share the same evolutionary history. This incongruence could be due to biased selection on some genes or recombination among mitochondrial genomes. The results thus suggest that the use of individual mitochondrial genes for phylogenetic inference could mask the true relationships between species in this complex.

Critiquing blind dating: the dangers of over-confident date estimates in comparative genomics

Phylogenomic advances provide more rigorous estimates for the timing of evolutionary divergences than previously available (e.g., Bayesian relaxed-clock estimates with soft fossil constraints). However, because many family-level clades and higher, as well as model species within those clades, have not been included in phylogenomic studies, the literature presents temporal estimates likely harboring substantial errors. Blindly using such dates can substantially retard scientific advancement. We suggest a way forward by conducting analyses that minimize prior assumptions and use large datasets, and demonstrate how using such a phylogenomic approach can lead to significantly more parsimonious conclusions without a good fossil record. We suggest that such an approach calls for research into the biological causes of conflict between molecular and fossil signatures.

Hybrid ancestry of an island subspecies of Galápagos mockingbird explains discordant gene trees

Introgression of genes through hybridization has been proposed to be an important driver of speciation, but in animals this has been shown only in relatively few cases until recently. Additionally, introgressive hybridization among non-sister species leads to a change in the gene tree topology of the concerned loci and thus complicates phylogenetic reconstruction. However, such cases of ancient introgression have been very difficult to demonstrate in birds. Here, we present such an example in an island bird subspecies, the Genovesa mockingbird (Mimus parvulus bauri). We assessed phylogenetic relationships and population structure among mockingbirds of the Galápagos archipelago using mitochondrial and nuclear DNA sequences, autosomal microsatellites, and morphological measurements. Mitochondrial haplotypes of Genovesa mockingbirds clustered closely with the haplotypes from two different species, San Cristóbal (M. melanotis) and Española (M. macdonaldi) mockingbirds. The same pattern was found for some haplotypes of two nuclear gene introns, while the majority of nuclear haplotypes of Genovesa mockingbirds were shared with other populations of the same species (M. parvulus). At 26 autosomal microsatellites, Genovesa mockingbirds grouped with other M. parvulus populations. This pattern shows that Genovesa mockingbirds contain mitochondria and some autosomal alleles that have most likely introgressed from M. melanotis into a largely M. parvulus background, making Genovesa mockingbirds a lineage of mixed ancestry, possibly undergoing speciation. Consistent with this hypothesis, mockingbirds on Genovesa are more clearly differentiated morphologically from other M. parvulus populations than M. melanotis is from M. parvulus.

Phylogenetic signal variation in the genomes of Medicago (Fabaceae)

Genome-scale data offer the opportunity to clarify phylogenetic relationships that are difficult to resolve with few loci, but they can also identify genomic regions with evolutionary history distinct from that of the species history. We collected whole-genome sequence data from 29 taxa in the legume genus Medicago, then aligned these sequences to the Medicago truncatula reference genome to confidently identify 87 596 variable homologous sites. We used this data set to estimate phylogenetic relationships among Medicago species, to investigate the number of sites needed to provide robust phylogenetic estimates and to identify specific genomic regions supporting topologies in conflict with the genome-wide phylogeny. Our full genomic data set resolves relationships within the genus that were previously intractable. Subsampling the data reveals considerable variation in phylogenetic signal and power in smaller subsets of the data. Even when sampling 5000 sites, no random sample of the data supports a topology identical to that of the genome-wide phylogeny. Phylogenetic relationships estimated from 500-site sliding windows revealed genome regions supporting several alternative species relationships among recently diverged taxa, consistent with the expected effects of deep coalescence or introgression in the recent history of Medicago.

Mitochondrial genome sequences reveal deep divergences among Anopheles punctulatus sibling species in Papua New Guinea

Background

Members of the Anopheles punctulatus group (AP group) are the primary vectors of human malaria in Papua New Guinea. The AP group includes 13 sibling species, most of them morphologically indistinguishable. Understanding why only certain species are able to transmit malaria requires a better comprehension of their evolutionary history. In particular, understanding relationships and divergence times among Anopheles species may enable assessing how malaria-related traits (e.g. blood feeding behaviours, vector competence) have evolved.

Methods

DNA sequences of 14 mitochondrial (mt) genomes from five AP sibling species and two species of the Anopheles dirus complex of Southeast Asia were sequenced. DNA sequences from all concatenated protein coding genes (10,770 bp) were then analysed using a Bayesian approach to reconstruct phylogenetic relationships and date the divergence of the AP sibling species.

Results

Phylogenetic reconstruction using the concatenated DNA sequence of all mitochondrial protein coding genes indicates that the ancestors of the AP group arrived in Papua New Guinea 25 to 54 million years ago and rapidly diverged to form the current sibling species.

Conclusion

Through evaluation of newly described mt genome sequences, this study has revealed a divergence among members of the AP group in Papua New Guinea that would significantly predate the arrival of humans in this region, 50 thousand years ago. The divergence observed among the mtDNA sequences studied here may have resulted from reproductive isolation during historical changes in sea-level through glacial minima and maxima. This leads to a hypothesis that the AP sibling species have evolved independently for potentially thousands of generations. This suggests that the evolution of many phenotypes, such as insecticide resistance will arise independently in each of the AP sibling species studied here.

Mathematical properties of the deep coalescence cost

In the minimizing-deep-coalescences (MDC) approach for species tree inference, a tree that has the minimal deep coalescence cost for reconciling a collection of gene trees is taken as an estimate of the species tree topology. The MDC method possesses the desirable Pareto property, and in practice it is quite accurate and computationally efficient. Here, in order to better understand the MDC method, we investigate some properties of the deep coalescence cost. We prove that the unit neighborhood of either a rooted species tree or a rooted gene tree under the deep coalescence cost is exactly the same as the tree's unit neighborhood under the rooted nearest-neighbor interchange (NNI) distance. Next, for a fixed species tree, we obtain the maximum deep coalescence cost across all gene trees as well as the number of gene trees that achieve the maximum cost. We also study corresponding problems for a fixed gene tree.

Gene flow networks among American Aedes aegypti populations

The mosquito Aedes aegypti, the dengue virus vector, has spread throughout the tropics in historical times. While this suggests man-mediated dispersal, estimating contemporary connectivity among populations has remained elusive. Here, we use a large mtDNA dataset and a Bayesian coalescent framework to test a set of hypotheses about gene flow among American Ae. aegypti populations. We assessed gene flow patterns at the continental and subregional (Amazon basin) scales. For the Americas, our data favor a stepping-stone model in which gene flow is higher among adjacent populations but in which, at the same time, North American and southeastern Brazilian populations are directly connected, likely via sea trade. Within Amazonia, the model with highest support suggests extensive gene flow among major cities; Manaus, located at the center of the subregional transport network, emerges as a potentially important connecting hub. Our results suggest substantial connectivity across Ae. aegypti populations in the Americas. As long-distance active dispersal has not been observed in this species, our data support man-mediated dispersal as a major determinant of the genetic structure of American Ae. aegypti populations. The inferred topology of interpopulation connectivity can inform network models of Ae. aegypti and dengue spread.

Phylogenetic species recognition reveals host-specific lineages among poplar rust fungi

Fungal species belonging to the genus Melampsora (Basidiomycota, Pucciniales) comprise rust pathogens that alternate between Salicaceae and other plant hosts. Species delineation and identification are difficult within this group due to the paucity of observable morphological features. Several Melampsora rusts are highly host-specific and this feature has been used for identification at the species level. However, this criterion is not always reliable since different Melampsora rust species can overlap on one host but specialize on a different one. To date, two different species recognition methods are used to recognize and define species within the Melampsora genus: (i) morphological species recognition, which is based solely on morphological criteria; and (ii) ecological species recognition, which combines morphological criteria with host range to recognize and define species. In order to clarify species recognition within the Melampsora genus, we applied phylogenetic species recognition to Melampsora poplar rusts by conducting molecular phylogenetic analyses on 15 Melampsora taxa using six nuclear and mitochondrial loci. By assessing the genealogical concordance between phylogenies, we identified 12 lineages that evolved independently, corresponding to distinct phylogenetic species. All 12 lineages were concordant with host specialization, but only three belonged to strictly defined morphological species. The estimation of the species tree obtained with Bayesian concordance analysis highlighted a potential co-evolutionary history between Melampsora species and their reciprocal aecial host plants. Within the Melampsora speciation process, aecial host may have had a strong effect on ancestral evolution, whereas telial host specificity seems to have evolved more recently. The morphological characters initially used to define species boundaries in the Melampsora genus are not reflective of the evolutionary and genetic relationships among poplar rusts. In order to construct a more meaningful taxonomy, host specificity must be considered an important criterion for delineating and describing species within the genus Melampsora as previously suggested by ecological species recognition.

Phylogenetically diverse groups of Bradyrhizobium isolated from nodules of Crotalaria spp., Indigofera spp., Erythrina brucei and Glycine max growing in Ethiopia

Ethiopian Bradyrhizobium strains isolated from root nodules of Crotalaria spp., Indigofera spp., Erythina brucei and soybean (Glycine max) represented genetically diverse phylogenetic groups of the genus Bradyrhizobium. Strains were characterized using the amplified fragment length polymorphism fingerprinting technique (AFLP) and multilocus sequence analysis (MLSA) of core and symbiotic genes. Based on phylogenetic analyses of concatenated recA-glnII-rpoB-16S rRNA genes sequences, Bradyrhizobium strains were distributed into fifteen phylogenetic groups under B. japonicum and B. elkanii super clades. Some of the isolates belonged to the species B. yuanmingense, B. elkanii and B. japonicum type I. However, the majority of the isolates represented unnamed Bradyrhizobium genospecies and of these, two unique lineages that most likely represent novel Bradyrhizobium species were identified among Ethiopian strains. The nodulation nodA gene sequence analysis revealed that all Ethiopian Bradyrhizobium isolates belonged to nodA sub-clade III.3. Strains were further classified into 14 groups together with strains from Africa, as well as some originating from the other tropical and subtropics regions. Strains were also clustered into 14 groups in nodY/K phylogeny similarly to the nodA tree. The nifH phylogenies of the Ethiopian Bradyrhizobium were generally also congruent with the nodA gene phylogeny, supporting the monophyletic origin of the symbiotic genes in Bradyrhizobium. The phylogenies of nodA and nifH genes were also partially congruent with that inferred from the concatenated core genes sequences, reflecting that the strains obtained their symbiotic genes vertically from their ancestor as well as horizontally from more distantly related Bradyrhizobium species.

Phylogenetic analysis of Babesia species in China based on Cytochrome b (COB) gene

In this study, a mitochondrial marker consisting of an approximately 550-bp region of the Cytochrome b genes (COB) was amplified by polymerase chain reaction (PCR) and sequenced from individual Babesia species. Sequence variation between Babesia species from China was 1.6-30.8%. The constructed phylogenetic tree based on the three unlinked gene sequences (partial COB gene, 18S rDNA and ITS) that evolve at different rates by the method of Neighbor-joining revealed the phylogenetic relationship of Babesia species in China compared with other published corresponding sequences from Babesia species. These data indicate that the 18S rDNA more reliably distinguish the deeper branches among some Babesia species than the partial COB gene and ITS, however, the partial COB gene sequence is better for recognizing close lineages among some Babesia species than the 18S rDNA and ITS sequences. So the combined phylogenetic analysis based on the multiple unlinked loci with different evolving rates can facilitate to establish the more reliable phylogenetic relationship of the Babesia genus. The data could be applicable for the survey of parasite dynamics, epidemiological studies as well as prevention and control of the disease.

A comparative analysis of Y chromosome and mtDNA phylogenies of the Hylobates gibbons

Background

The evolutionary relationships of closely related species have long been of interest to biologists since these species experienced different evolutionary processes in a relatively short period of time. Comparison of phylogenies inferred from DNA sequences with differing inheritance patterns, such as mitochondrial, autosomal, and X and Y chromosomal loci, can provide more comprehensive inferences of the evolutionary histories of species. Gibbons, especially the genus Hylobates, are particularly intriguing as they consist of multiple closely related species which emerged rapidly and live in close geographic proximity. Our current understanding of relationships among Hylobates species is largely based on data from the maternally-inherited mitochondrial DNAs (mtDNAs).

Results

To infer the paternal histories of gibbon taxa, we sequenced multiple Y chromosomal loci from 26 gibbons representing 10 species. As expected, we find levels of sequence variation some five times lower than observed for the mitochondrial genome (mtgenome). Although our Y chromosome phylogenetic tree shows relatively low resolution compared to the mtgenome tree, our results are consistent with the monophyly of gibbon genera suggested by the mtgenome tree. In a comparison of the molecular dating of divergences and on the branching patterns of phylogeny trees between mtgenome and Y chromosome data, we found: 1) the inferred divergence estimates were more recent for the Y chromosome than for the mtgenome, 2) the species H. lar and H. pileatus are monophyletic in the mtgenome phylogeny, respectively, but a H. pileatus individual falls into the H. lar Y chromosome clade.

Conclusions

Based on the ~6.4 kb of Y chromosomal DNA sequence data generated for each of the 26 individuals in this study, we provide molecular inferences on gibbon and particularly on Hylobates evolution complementary to those from mtDNA data. Overall, our results illustrate the utility of comparative studies of loci with different inheritance patterns for investigating potential sex specific processes on the evolutionary histories of closely related taxa, and emphasize the need for further sampling of gibbons of known provenance.

Consensus properties for the deep coalescence problem and their application for scalable tree search

Background

To infer a species phylogeny from unlinked genes, phylogenetic inference methods must confront the biological processes that create incongruence between gene trees and the species phylogeny. Intra-specific gene variation in ancestral species can result in deep coalescence, also known as incomplete lineage sorting, which creates incongruence between gene trees and the species tree. One approach to account for deep coalescence in phylogenetic analyses is the deep coalescence problem, which takes a collection of gene trees and seeks the species tree that implies the fewest deep coalescence events. Although this approach is promising for phylogenetics, the consensus properties of this problem are mostly unknown and analyses of large data sets may be computationally prohibitive.

Results

We prove that the deep coalescence consensus tree problem satisfies the highly desirable Pareto property for clusters (clades). That is, in all instances, each cluster that is present in all of the input gene trees, called a consensus cluster, will also be found in every optimal solution. Moreover, we introduce a new divide and conquer method for the deep coalescence problem based on the Pareto property. This method refines the strict consensus of the input gene trees, thereby, in practice, often greatly reducing the complexity of the tree search and guaranteeing that the estimated species tree will satisfy the Pareto property.

Conclusions

Analyses of both simulated and empirical data sets demonstrate that the divide and conquer method can greatly improve upon the speed of heuristics that do not consider the Pareto consensus property, while also guaranteeing that the proposed solution fulfills the Pareto property. The divide and conquer method extends the utility of the deep coalescence problem to data sets with enormous numbers of taxa.

A multi-locus phylogeny reveals a complex pattern of diversification related to climate and habitat heterogeneity in southern African white-eyes

The recent, rapid radiation of Zosteropidae, coupled with their high levels of colonizing ability and phenotypic diversity, makes species delimitation within this family problematic. Given these problems, challenges to establish the mechanisms driving diversity and speciation within this group have arisen. Four morphologically distinct southern African Zosterops taxa, with a contentious taxonomic past, provide such a challenge. Here, supplemented with morphological and environmental analytical techniques, a combination of mitochondrial and nuclear markers were analyzed using Bayesian and Likelihood methods to determine their speciation patterns and to establish the phylogenetic relationships of these four morphologically diverse southern African Zosterops taxa. Nearly all individuals were phenotypically diagnosable, even those individuals collected in areas of contact between taxa. Localities where two or more taxa co-occur appear to possess intermediate environmental characteristics. Initial Bayesian and Likelihood mitochondrial DNA analyses and Bayesian structure analyses of the combined nuclear markers indicated levels of hybridization in areas of sympatry. A combined mtDNA and nuclear DNA analysis and a species tree analysis (with hybrids excluded) placed Z. pallidus as sister to the other southern African taxa, with Z. senegalensis the putative sister taxon to a clade comprising Z. capensis and Z. virens. The grouping of taxon-specific sampling localities and the apparent intermediate nature of birds from areas of sympatry points toward an influence of habitat type and the associated climatic conditions in driving Zosterops diversification in southern Africa.

Multiple merger gene genealogies in two species: Monophyly, paraphyly, and polyphyly for two examples of Lambda coalescents

Probabilities of monophyly, paraphyly, and polyphyly of two-species gene genealogies are computed for modest sample sizes and compared for two different Λ coalescent processes. Coalescent processes belonging to the Λ coalescent family admit asynchronous multiple mergers of active ancestral lineages. Assigning a timescale to the time of divergence becomes a central issue when different populations have different coalescent processes running on different timescales. Clade probabilities in single populations are also computed, which can be useful for testing for taxonomic distinctiveness of an observed set of monophyletic lineages. The coalescence rates of multiple merger coalescent processes are functions of coalescent parameters. The effect of coalescent parameters on the probabilities studied depends on the coalescent process, and if the population is ancestral or derived. The probability of reciprocal monophyly tends to be somewhat lower, when associated with a Λ coalescent, under the null hypothesis that two groups come from the same population. However, even for fairly recent divergence times, the probability of monophyly tends to be higher as a function of the number of generations for coalescent processes that admit multiple mergers, and is sensitive to the parameter of one of the example processes.

Multiple mitochondrial introgression events and heteroplasmy in trypanosoma cruzi revealed by maxicircle MLST and next generation sequencing

Background

Mitochondrial DNA is a valuable taxonomic marker due to its relatively fast rate of evolution. In Trypanosoma cruzi, the causative agent of Chagas disease, the mitochondrial genome has a unique structural organization consisting of 20–50 maxicircles (∼20 kb) and thousands of minicircles (0.5–10 kb). T. cruzi is an early diverging protist displaying remarkable genetic heterogeneity and is recognized as a complex of six discrete typing units (DTUs). The majority of infected humans are asymptomatic for life while 30–35% develop potentially fatal cardiac and/or digestive syndromes. However, the relationship between specific clinical outcomes and T. cruzi genotype remains elusive. The availability of whole genome sequences has driven advances in high resolution genotyping techniques and re-invigorated interest in exploring the diversity present within the various DTUs.

Methodology/Principal Findings

To describe intra-DTU diversity, we developed a highly resolutive maxicircle multilocus sequence typing (mtMLST) scheme based on ten gene fragments. A panel of 32 TcI isolates was genotyped using the mtMLST scheme, GPI, mini-exon and 25 microsatellite loci. Comparison of nuclear and mitochondrial data revealed clearly incongruent phylogenetic histories among different geographical populations as well as major DTUs. In parallel, we exploited read depth data, generated by Illumina sequencing of the maxicircle genome from the TcI reference strain Sylvio X10/1, to provide the first evidence of mitochondrial heteroplasmy (heterogeneous mitochondrial genomes in an individual cell) in T. cruzi.

Conclusions/Significance

mtMLST provides a powerful approach to genotyping at the sub-DTU level. This strategy will facilitate attempts to resolve phenotypic variation in T. cruzi and to address epidemiologically important hypotheses in conjunction with intensive spatio-temporal sampling. The observations of both general and specific incidences of nuclear-mitochondrial phylogenetic incongruence indicate that genetic recombination is geographically widespread and continues to influence the natural population structure of TcI, a conclusion which challenges the traditional paradigm of clonality in T. cruzi.

Chagas disease, caused by the protozoan parasite Trypanosoma cruzi, is an important public health problem in Latin America. While molecular techniques can differentiate the major T. cruzi genetic lineages, few have sufficient resolution to describe diversity among closely related strains. The online availability of three mitochondrial genomes allowed us to design a multilocus sequence typing (mtMLST) scheme to exploit these rapidly evolving markers. We compared mtMLST with current nuclear typing tools using isolates belonging to the oldest and most widely occurring lineage TcI. T. cruzi is generally believed to reproduce clonally. However, in this study, distinct branching patterns between mitochondrial and nuclear phylogenetic trees revealed multiple incidences of genetic exchange within different geographical populations and major lineages. We also examined Illumina sequencing data from the TcI genome strain which revealed multiple different mitochondrial genomes within an individual parasite (heteroplasmy) that were, however, not sufficiently divergent to represent a major source of typing error. We strongly recommend this combined nuclear and mitochondrial genotyping methodology to reveal cryptic diversity and genetic exchange in T. cruzi. The level of resolution that this mtMLST provides should greatly assist attempts to elucidate the complex interactions between parasite genotype, clinical outcome and disease distribution.

Phylogeny of bradyrhizobia from Chinese cowpea miscellany inferred from 16S rRNA, atpD, glnII, and 16S-23S intergenic spacer sequences

The cowpea (Vigna unguiculata L.), peanut (Arachis hypogaea L.), and mung bean (Vigna radiata L.) belong to a group of plants known as the "cowpea miscellany" plants, which are widely cultivated throughout the tropic and subtropical zones of Africa and Asia. However, the phylogeny of the rhizobial strains that nodulate these plants is poorly understood. Previous studies have isolated a diversity of rhizobial strains from cowpea miscellany hosts and have suggested that, phylogenetically, they are from different species. In this work, the phylogeny of 42 slow-growing rhizobial strains, isolated from root nodules of cowpea, peanut, and mung bean from different geographical regions of China, was investigated using sequences from the 16S rRNA, atpD and glnII genes, and the 16S-23S rRNA intergenic spacer. The indigenous rhizobial strains from the cowpea miscellany could all be placed in the genus Bradyrhizobium , and Bradyrhizobium liaoningense and Bradyrhizobium yuanmingense were the main species. Phylogenies derived from housekeeping genes were consistent with phylogenies generated from the ribosomal gene. Mung bean rhizobia clustered only into B. liaoningense and B. yuanmingense and were phylogenetically less diverse than cowpea and peanut rhizobia. Geographical origin was significantly reflected in the phylogeny of mung bean rhizobia. Most cowpea rhizobia were more closely related to the 3 major groups B. liaoningense, B. yuanmingense, and Bradyrhizobium elkanii than to the minor groups Bradyrhizobium japonicum or Bradyrhizobium canariense . However, most peanut rhizobia were more closely related to the 2 major groups B. liaoningense and B. yuanmingense than to the minor group B. elkanii.

Impact of deep coalescence on the reliability of species tree inference from different types of DNA markers in mammals

An important challenge for phylogenetic studies of closely related species is the existence of deep coalescence and gene tree heterogeneity. However, their effects can vary between species and they are often neglected in phylogenetic analyses. In addition, a practical problem in the reconstruction of shallow phylogenies is to determine the most efficient set of DNA markers for a reliable estimation. To address these questions, we conducted a multilocus simulation study using empirical values of nucleotide diversity and substitution rates obtained from a wide range of mammals and evaluated the performance of both gene tree and species tree approaches to recover the known speciation times and topological relationships. We first show that deep coalescence can be a serious problem, more than usually assumed, for the estimation of speciation times in mammals using traditional gene trees. Furthermore, we tested the performance of different sets of DNA markers in the determination of species trees using a coalescent approach. Although the best estimates of speciation times were obtained, as expected, with the use of an increasing number of nuclear loci, our results show that similar estimations can be obtained with a much lower number of genes and the incorporation of a mitochondrial marker, with its high information content. Thus, the use of the combined information of both nuclear and mitochondrial markers in a species tree framework is the most efficient option to estimate recent speciation times and, consequently, the underlying species tree.

DNA barcoding of recently diverged species: relative performance of matching methods

Recently diverged species are challenging for identification, yet they are frequently of special interest scientifically as well as from a regulatory perspective. DNA barcoding has proven instrumental in species identification, especially in insects and vertebrates, but for the identification of recently diverged species it has been reported to be problematic in some cases. Problems are mostly due to incomplete lineage sorting or simply lack of a 'barcode gap' and probably related to large effective population size and/or low mutation rate. Our objective was to compare six methods in their ability to correctly identify recently diverged species with DNA barcodes: neighbor joining and parsimony (both tree-based), nearest neighbor and BLAST (similarity-based), and the diagnostic methods DNA-BAR, and BLOG. We analyzed simulated data assuming three different effective population sizes as well as three selected empirical data sets from published studies. Results show, as expected, that success rates are significantly lower for recently diverged species (∼75%) than for older species (∼97%) (P<0.00001). Similarity-based and diagnostic methods significantly outperform tree-based methods, when applied to simulated DNA barcode data (P<0.00001). The diagnostic method BLOG had highest correct query identification rate based on simulated (86.2%) as well as empirical data (93.1%), indicating that it is a consistently better method overall. Another advantage of BLOG is that it offers species-level information that can be used outside the realm of DNA barcoding, for instance in species description or molecular detection assays. Even though we can confirm that identification success based on DNA barcoding is generally high in our data, recently diverged species remain difficult to identify. Nevertheless, our results contribute to improved solutions for their accurate identification.

iGLASS: an improvement to the GLASS method for estimating species trees from gene trees

Several methods have been designed to infer species trees from gene trees while taking into account gene tree/species tree discordance. Although some of these methods provide consistent species tree topology estimates under a standard model, most either do not estimate branch lengths or are computationally slow. An exception, the GLASS method of Mossel and Roch, is consistent for the species tree topology, estimates branch lengths, and is computationally fast. However, GLASS systematically overestimates divergence times, leading to biased estimates of species tree branch lengths. By assuming a multispecies coalescent model in which multiple lineages are sampled from each of two taxa at L independent loci, we derive the distribution of the waiting time until the first interspecific coalescence occurs between the two taxa, considering all loci and measuring from the divergence time. We then use the mean of this distribution to derive a correction to the GLASS estimator of pairwise divergence times. We show that our improved estimator, which we call iGLASS, consistently estimates the divergence time between a pair of taxa as the number of loci approaches infinity, and that it is an unbiased estimator of divergence times when one lineage is sampled per taxon. We also show that many commonly used clustering methods can be combined with the iGLASS estimator of pairwise divergence times to produce a consistent estimator of the species tree topology. Through simulations, we show that iGLASS can greatly reduce the bias and mean squared error in obtaining estimates of divergence times in a species tree.

Molecular delineation of species in the coral holobiont

The coral holobiont is a complex assemblage of organisms spanning a diverse taxonomic range including a cnidarian host, as well as various dinoflagellate, prokaryotic and acellular symbionts. With the accumulating information on the molecular diversity of these groups, binomial species classification and a reassessment of species boundaries for the partners in the coral holobiont is a logical extension of this work and will help enhance the capacity for comparative research among studies. To aid in this endeavour, we review the current literature on species diversity for the three best studied partners of the coral holobiont (coral, Symbiodinium, prokaryotes) and provide suggestions for future work on systematics within these taxa. We advocate for an integrative approach to the delineation of species using both molecular genetics in combination with phenetic characters. We also suggest that an a priori set of criteria be developed for each taxonomic group as no one species concept or accompanying set of guidelines is appropriate for delineating all members of the coral holobiont.

Chasing the hare - evaluating the phylogenetic utility of a nuclear single copy gene region at and below species level within the species rich group Peperomia (Piperaceae)

BACKGROUND

The rapidly increasing number of available plant genomes opens up almost unlimited prospects for biology in general and molecular phylogenetics in particular. A recent study took advantage of this data and identified a set of nuclear genes that occur in single copy in multiple sequenced angiosperms. The present study is the first to apply genomic sequence of one of these low copy genes, agt1, as a phylogenetic marker for species-level phylogenetics. Its utility is compared to the performance of several coding and non-coding chloroplast loci that have been suggested as most applicable for this taxonomic level. As a model group, we chose Tildenia, a subgenus of Peperomia (Piperaceae), one of the largest plant genera. Relationships are particularly difficult to resolve within these species rich groups due to low levels of polymorphisms and fast or recent radiation. Therefore, Tildenia is a perfect test case for applying new phylogenetic tools.

RESULTS

We show that the nuclear marker agt1, and in particular the agt1 introns, provide a significantly increased phylogenetic signal compared to chloroplast markers commonly used for low level phylogenetics. 25% of aligned characters from agt1 intron sequence are parsimony informative. In comparison, the introns and spacer of several common chloroplast markers (trnK intron, trnK-psbA spacer, ndhF-rpl32 spacer, rpl32-trnL spacer, psbA-trnH spacer) provide less than 10% parsimony informative characters. The agt1 dataset provides a deeper resolution than the chloroplast markers in Tildenia.

CONCLUSIONS

Single (or very low) copy nuclear genes are of immense value in plant phylogenetics. Compared to other nuclear genes that are members of gene families of all sizes, lab effort, such as cloning, can be kept to a minimum. They also provide regions with different phylogenetic content deriving from coding and non-coding parts of different length. Thus, they can be applied to a wide range of taxonomic levels from family down to population level. As more plant genomes are sequenced, we will obtain increasingly precise information about which genes return to single copy most rapidly following gene duplication and may be most useful across a wide range of plant groups.

Evolution of vacuolar proton pyrophosphatase domains and volutin granules: clues into the early evolutionary origin of the acidocalcisome

Background

Volutin granules appear to be universally distributed and are morphologically and chemically identical to acidocalcisomes, which are electron-dense granular organelles rich in calcium and phosphate, whose functions include storage of phosphorus and various metal ions, metabolism of polyphosphate, maintenance of intracellular pH, osmoregulation and calcium homeostasis. Prokaryotes are thought to differ from eukaryotes in that they lack membrane-bounded organelles. However, it has been demonstrated that as in acidocalcisomes, the calcium and polyphosphate-rich intracellular "volutin granules (polyphosphate bodies)" in two bacterial species, Agrobacterium tumefaciens, and Rhodospirillum rubrum, are membrane bound and that the vacuolar proton-translocating pyrophosphatases (V-H⁺PPases) are present in their surrounding membranes. Volutin granules and acidocalcisomes have been found in organisms as diverse as bacteria and humans.

Results

Here, we show volutin granules also occur in Archaea and are, therefore, present in the three superkingdoms of life (Archaea, Bacteria and Eukarya). Molecular analyses of V-H⁺PPase pumps, which acidify the acidocalcisome lumen and are diagnostic proteins of the organelle, also reveal the presence of this enzyme in all three superkingdoms suggesting it is ancient and universal. Since V-H⁺PPase sequences contained limited phylogenetic signal to fully resolve the ancestral nodes of the tree, we investigated the divergence of protein domains in the V-H⁺PPase molecules. Using Protein family (Pfam) database, we found a domain in the protein, PF03030. The domain is shared by 31 species in Eukarya, 231 in Bacteria, and 17 in Archaea. The universal distribution of the V-H⁺PPase PF03030 domain, which is associated with the V-H⁺PPase function, suggests the domain and the enzyme were already present in the Last Universal Common Ancestor (LUCA).

Conclusion

The importance of the V-H⁺PPase function and the evolutionary dynamics of these domains support the early origin of the acidocalcisome organelle. In particular, the universality of volutin granules and presence of a functional V-H⁺PPase domain in the three superkingdoms of life reveals that the acidocalcisomes may have appeared earlier than the divergence of the superkingdoms. This result is remarkable and highlights the possibility that a high degree of cellular compartmentalization could already have been present in the LUCA.

Reviewers

This article was reviewed by Anthony Poole, Lakshminarayan Iyer and Daniel Kahn

A hitchhikers guide to the Galápagos: co-phylogeography of Galápagos mockingbirds and their parasites

Background

Parasites are evolutionary hitchhikers whose phylogenies often track the evolutionary history of their hosts. Incongruence in the evolutionary history of closely associated lineages can be explained through a variety of possible events including host switching and host independent speciation. However, in recently diverged lineages stochastic population processes, such as retention of ancestral polymorphism or secondary contact, can also explain discordant genealogies, even in fully co-speciating taxa. The relatively simple biogeographic arrangement of the Galápagos archipelago, compared with mainland biomes, provides a framework to identify stochastic and evolutionary informative components of genealogic data in these recently diverged organisms.

Results

Mitochondrial DNA sequences were obtained for four species of Galápagos mockingbirds and three sympatric species of ectoparasites - two louse and one mite species. These data were complemented with nuclear EF1α sequences in selected samples of parasites and with information from microsatellite loci in the mockingbirds. Mitochondrial sequence data revealed differences in population genetic diversity between all taxa and varying degrees of topological congruence between host and parasite lineages. A very low level of genetic variability and lack of congruence was found in one of the louse parasites, which was excluded from subsequent joint analysis of mitochondrial data. The reconciled multi-species tree obtained from the analysis is congruent with both the nuclear data and the geological history of the islands.

Conclusions

The gene genealogies of Galápagos mockingbirds and two of their ectoparasites show strong phylogeographic correlations, with instances of incongruence mostly explained by ancestral genetic polymorphism. A third parasite genealogy shows low levels of genetic diversity and little evidence of co-phylogeny with their hosts. These differences can mostly be explained by variation in life-history characteristics, primarily host specificity and dispersal capabilities. We show that pooling genetic data from organisms living in close ecological association reveals a more accurate phylogeographic history for these taxa. Our results have implications for the conservation and taxonomy of Galápagos mockingbirds and their parasites.

Y-chromosome phylogeny in the evolutionary net of chamois (genus Rupicapra)

Background

The chamois, distributed over most of the medium to high altitude mountain ranges of southern Eurasia, provides an excellent model for exploring the effects of historical and evolutionary events on diversification. Populations have been grouped into two species, Rupicapra pyrenaica from southwestern Europe and R. rupicapra from eastern Europe. The study of matrilineal mitochondrial DNA (mtDNA) and biparentally inherited microsatellites showed that the two species are paraphyletic and indicated alternate events of population contraction and dispersal-hybridization in the diversification of chamois. Here we investigate the pattern of variation of the Y-chromosome to obtain information on the patrilineal phylogenetic position of the genus Rupicapra and on the male-specific dispersal of chamois across Europe.

Results

We analyzed the Y-chromosome of 87 males covering the distribution range of the Rupicapra genus. We sequenced a fragment of the SRY gene promoter and characterized the male specific microsatellites UMN2303 and SRYM18. The SRY promoter sequences of two samples of Barbary sheep (Ammotragus lervia) were also determined and compared with the sequences of Bovidae available in the GenBank. Phylogenetic analysis of the alignment showed the clustering of Rupicapra with Capra and the Ammotragus sequence obtained in this study, different from the previously reported sequence of Ammotragus which groups with Ovis. Within Rupicapra, the combined data define 10 Y-chromosome haplotypes forming two haplogroups, which concur with taxonomic classification, instead of the three clades formed for mtDNA and nuclear microsatellites. The variation shows a west-to-east geographical cline of ancestral to derived alleles.

Conclusions

The phylogeny of the SRY-promoter shows an association between Rupicapra and Capra. The position of Ammotragus needs a reinvestigation. The study of ancestral and derived characters in the Y-chromosome suggests that, contrary to the presumed Asian origin, the paternal lineage of chamois originated in the Mediterranean, most probably in the Iberian Peninsula, and dispersed eastwards through serial funding events during the glacial-interglacial cycles of the Quaternary. The diversity of Y-chromosomes in chamois is very low. The differences in patterns of variation among Y-chromosome, mtDNA and biparental microsatellites reflect the evolutionary characteristics of the different markers as well as the effects of sex-biased dispersal and species phylogeography.

Adaptation, isolation by distance and human-mediated transport determine patterns of gene flow among populations of the disease vector Aedes taeniorhynchus in the Galapagos Islands

The black salt-marsh mosquito (Aedes taeniorhynchus) is the only native mosquito in the Galapagos Islands and potentially a major disease vector for Galapagos wildlife. Little is known about its population structure, or how its dynamics may be influenced by human presence in the archipelago. We used microsatellite data to assess the structure and patterns of A. taeniorhynchus gene flow among and within islands, to identify potential barriers to mosquito dispersal, and to investigate human-aided transport of mosquitoes across the archipelago. Our results show that inter-island migration of A. taeniorhynchus occurs frequently on an isolation by distance basis. High levels of inter-island migration were detected amongst the major ports of the archipelago, strongly suggesting the occurrence of human-aided transport of mosquitoes among islands, underlining the need for strict control measures to avoid the transport of disease vectors between islands. The prevalence of filarial nematode infection in Galapagos flightless cormorants is correlated with the population structure and migration patterns of A. taeniorhynchus, suggesting that A. taeniorhynchus is an important vector of this arthropod-borne parasite in the Galapagos Islands. Therefore mosquito population structure in Galapagos may have the potential to influence mosquito-borne parasite population dynamics, and the subsequent impacts of such pathogens on their host species in the islands.

Phylogeny and divergence times inferred from rps16 sequence data analyses for Tricyrtis (Liliaceae), an endemic genus of north-east Asia

BACKGROUND AND AIMS

Tricyrtis is a genus of monocots with attractive and sophisticated flower shapes and colours, endemic to north-east Asia. There are 18 known species. The highly restricted geographical distribution of the genus is of great interest in terms of both abiotic (continental drift) and biotic (long-distance dispersal) impacts on monocot plant speciation events and their timing, and of evolutionary patterns of diversification leading to the extant taxa. The aims of this study were to (i) predict the time of speciation (divergence) events at infraspecific levels of Tricyrtis, (ii) estimate the rate of evolution of the genus and (iii) provide information on an excellent plant model system in terms of studying loss of biodiversity or extinction of organisms in the dynamic earth environment.

METHODOLOGY

To investigate the divergence time and evolution rate of Tricyrtis, Bayesian Markov chain Monte Carlo (MCMC) analyses were performed by calculating the mean branch lengths of evolutionary paths based on base substitution variations between rps16 intron nucleotide sequences from the 18 known species.

PRINCIPAL RESULTS

Based upon the relaxed molecular clock model test data, a Bayesian phylogenetic inference tree is presented, and the divergence times and rate of evolution of Tricyrtis were estimated. The analyses also suggest that evolution is occurring at the infraspecific level of the genus in a manner that is not strictly clock bound.

CONCLUSIONS

Continental drift may have been the main speciation process giving rise to the current distribution of the taxa of Tricyrtis. The single-locus gene sequence data presented here are a significant step towards an improved future understanding of the molecular evolution of Tricyrtis via multi-locus evaluation.