Multigene Phylogenetics Reveals Temporal Diversification of Major African Malaria Vectors

The major vectors of malaria in sub-Saharan Africa belong to subgenus Cellia. Yet, phylogenetic relationships and temporal diversification among African mosquito species have not been unambiguously determined. Knowledge about vector evolutionary history is crucial for correct interpretation of genetic changes identified through comparative genomics analyses. In this study, we estimated a molecular phylogeny using 49 gene sequences for the African malaria vectors An. gambiae, An. funestus, An. nili, the Asian malaria mosquito An. stephensi, and the outgroup species Culex quinquefasciatus and Aedes aegypti. To infer the phylogeny, we identified orthologous sequences uniformly distributed approximately every 5 Mb in the five chromosomal arms. The sequences were aligned and the phylogenetic trees were inferred using maximum likelihood and neighbor-joining methods. Bayesian molecular dating using a relaxed log normal model was used to infer divergence times. Trees from individual genes agreed with each other, placing An. nili as a basal clade that diversified from the studied malaria mosquito species 47.6 million years ago (mya). Other African malaria vectors originated more recently, and independently acquired traits related to vectorial capacity. The lineage leading to An. gambiae diverged 30.4 mya, while the African vector An. funestus and the Asian vector An. stephensi were the most closely related sister taxa that split 20.8 mya. These results were supported by consistently high bootstrap values in concatenated phylogenetic trees generated individually for each chromosomal arm. Genome-wide multigene phylogenetic analysis is a useful approach for discerning historic relationships among malaria vectors, providing a framework for the correct interpretation of genomic changes across species, and comprehending the evolutionary origins of this ubiquitous and deadly insect-borne disease.     

Molecular identification of yeasts from the order Trichosporonales causing superficial infections

BackgroundThe molecular reclassification of the order Trichosporonales placed the medically relevant Trichosporon species into three genera of the family Trichosporonaceae: Cutaneotrichosporon, Trichosporon, and Apiotrichum. From the clinical and epidemiological standpoint, it is important to identify any species of the family Trichosporonaceae because they present different antifungal susceptibility profiles. In Mexico, little is known about trichosporonosis etiology because the fungi are identified through phenotypic methods.AimsTo identify at a molecular level 12 yeast isolates morfologically compatible with Trichosporon, obtained from patients with superficial infections.MethodsThe yeast isolates were obtained from patients with white piedra, onychomycosis, and hand and foot dermatomycosis, and were identified morphologically and genotypically (sequencing of the IGS1 region and phylogenetic analysis using the Maximum Likelihood Method). The phylogenetic analysis included 40 yeast sequences from the order Trichosporonales and one from Cryptococcus neoformans as outgroup.ResultsBased on the molecular analysis, we identified three (25%) Trichosporon inkin isolates, two (16.7%) Trichosporon asteroides, two (16.7%) Cutaneotrichosporon mucoides, and one each (8.3%) of Trichosporon aquatile, Trichosporon asahii, Apiotrichum montevideense, Cutaneotrichosporon cutaneum, and Cutaneotrichosporon jirovecii.ConclusionsThe molecular characterization of the isolates showed a broad diversity of species within the order Trichosporonales, particularly among onychomycosis. It is essential to identify these yeasts at the species level to delve into their epidemiology.     

Cryptic Genetic Diversity Is Paramount in Small-Bodied Amphibians of the Genus Euparkerella (Anura: Craugastoridae) Endemic to the Brazilian Atlantic Forest

Morphological similarity associated to restricted distributions and low dispersal abilities make the direct developing “Terrarana” frogs of the genus Euparkerella a good model for examining diversification processes. We here infer phylogenetic relationships within the genus Euparkerella, using DNA sequence data from one mitochondrial and four nuclear genes coupled with traditional Bayesian phylogenetic reconstruction approaches and more recent coalescent methods of species tree inference. We also used Bayesian clustering analysis and a recent Bayesian coalescent-based approach specifically to infer species delimitation. The analysis of 39 individuals from the four known Euparkerella species uncovered high levels of genetic diversity, especially within the two previously morphologically-defined E. cochranae and E. brasiliensis. Within these species, the gene trees at five independent loci and trees from combined data (concatenated dataset and the species tree) uncovered six deeply diverged and geographically coherent evolutionary units, which may have diverged between the Miocene and the Pleistocene. These six units were also uncovered in the Bayesian clustering analysis, and supported by the Bayesian coalescent-based species delimitation (BPP), and Genealogical Sorting Index (GSI), providing thus strong evidence for underestimation of the current levels of diversity within Euparkerella. The cryptic diversity now uncovered opens new opportunities to examine the origins and maintenance of microendemism in the context of spatial heterogeneity and/or human induced fragmentation of the highly threatened Brazilian Atlantic forest hotspot.     

The Tree versus the Forest: The Fungal Tree of Life and the Topological Diversity within the Yeast Phylome

A recurrent topic in phylogenomics is the combination of various sequence alignments to reconstruct a tree that describes the evolutionary relationships within a group of species. However, such approach has been criticized for not being able to properly represent the topological diversity found among gene trees. To evaluate the representativeness of species trees based on concatenated alignments, we reconstruct several fungal species trees and compare them with the complete collection of phylogenies of genes encoded in the Saccharomyces cerevisiae genome. We found that, despite high levels of among-gene topological variation, the species trees do represent widely supported phylogenetic relationships. Most topological discrepancies between gene and species trees are concentrated in certain conflicting nodes. We propose to map such information on the species tree so that it accounts for the levels of congruence across the genome. We identified the lack of sufficient accuracy of current alignment and phylogenetic methods as an important source for the topological diversity encountered among gene trees. Finally, we discuss the implications of the high levels of topological variation for phylogeny-based orthology prediction strategies.     

Phylogenetic Relationships and Species Delimitation in Pinus Section Trifoliae Inferrred from Plastid DNA

Recent diversification followed by secondary contact and hybridization may explain complex patterns of intra- and interspecific morphological and genetic variation in the North American hard pines (Pinus section Trifoliae), a group of approximately 49 tree species distributed in North and Central America and the Caribbean islands. We concatenated five plastid DNA markers for an average of 3.9 individuals per putative species and assessed the suitability of the five regions as DNA bar codes for species identification, species delimitation, and phylogenetic reconstruction. The ycf1 gene accounted for the greatest proportion of the alignment (46.9%), the greatest proportion of variable sites (74.9%), and the most unique sequences (75 haplotypes). Phylogenetic analysis recovered clades corresponding to subsections Australes, Contortae, and Ponderosae. Sequences for 23 of the 49 species were monophyletic and sequences for another 9 species were paraphyletic. Morphologically similar species within subsections usually grouped together, but there were exceptions consistent with incomplete lineage sorting or introgression. Bayesian relaxed molecular clock analyses indicated that all three subsections diversified relatively recently during the Miocene. The general mixed Yule-coalescent method gave a mixed model estimate of only 22 or 23 evolutionary entities for the plastid sequences, which corresponds to less than half the 49 species recognized based on morphological species assignments. Including more unique haplotypes per species may result in higher estimates, but low mutation rates, recent diversification, and large effective population sizes may limit the effectiveness of this method to detect evolutionary entities.     

A genome sequence-based approach to taxonomy of the genus Nocardia

The genus Nocardia includes both pathogens and producers of useful secondary metabolites. Although 16S rRNA analysis is required to accurately discriminate among phylogenetic relationships of the Nocardia species, most branches of 16S rRNA-based phylogenetic trees are not reliable. In this study, we performed in silico analyses of the genome sequences of Nocardia species in order to understand their diversity and classification for their identification and applications. Draft genome sequences of 26 Nocardia strains were determined. Phylogenetic trees were prepared on the basis of multilocus sequence analysis of the concatenated sequences of 12 genes (atpD-dnaJ-groL1-groL2-gyrB-recA-rpoA-secA-secY-sodA-trpB-ychF) and a bidirectional best hit. To elucidate the evolutionary relationships of these genes, the genome-to-genome distance was investigated on the basis of the average nucleotide identity, DNA maximal unique matches index, and genome-to-genome distance calculator. The topologies of all phylogenetic trees were found to be essentially similar to each other. Furthermore, whole genome-derived and multiple gene-derived relationships were found to be suitable for extensive intra-genus assessment of the genus Nocardia.     

Multilocus sequence analysis of the actinobacterial genus Kribbella

Multilocus sequence analysis (MLSA) was used to refine the phylogenetic analysis of the genus Kribbella, which currently contains 17 species with validly-published names. Sequences were obtained for the 16S rRNA, gyrB, rpoB, recA, relA and atpD genes for 16 of the 17 type strains of the genus plus seven non-type strains. A five-gene concatenated sequence of 4099 nt was used to examine the phylogenetic relationships between the species of the genus Kribbella. Using the concatenated sequence of the gyrB-rpoB-recA-relA and atpD genes, most Kribbella type strains can be distinguished by a genetic distance of >0.04. Each single-gene tree had an overall topology similar to that of the concatenated sequence tree. The single-gene relA tree, used here for the first time in MLSA of actinobacteria, had good bootstrap support, comparable to the rpoB and atpD gene trees, which had topologies closest to that of the concatenated sequence tree. This illustrates that relA is a useful addition in MLSA studies of the genus Kribbella. We propose that concatenated gyrB-rpoB-recA-relA-atpD gene sequences be used for examining the phylogenetic relationships within the genus Kribbella and for determining the closest phylogenetic relatives to be used for taxonomic comparisons.     

Exploring the tempo of species diversification in legumes

Whatever criteria are used to measure evolutionary success – species numbers, geographic range, ecological abundance, ecological and life history diversity, background diversification rates, or the presence of rapidly evolving clades – the legume family is one of the most successful lineages of flowering plants. Despite this, we still know rather little about the dynamics of lineage and species diversification across the family through the Cenozoic, or about the underlying drivers of diversification. There have been few attempts to estimate net species diversification rates or underlying speciation and extinction rates for legume clades, to test whether among-lineage variation in diversification rates deviates from null expectations, or to locate species diversification rate shifts on specific branches of the legume phylogenetic tree. In this study, time-calibrated phylogenetic trees for a set of species-rich legume clades – Calliandra, Indigofereae, Lupinus, Mimosa and Robinieae – and for the legume family as a whole, are used to explore how we might approach these questions. These clades are analysed using recently developed maximum likelihood and Bayesian methods to detect species diversification rate shifts and test for among-lineage variation in speciation, extinction and net diversification rates. Possible explanations for rate shifts in terms of extrinsic factors and/or intrinsic trait evolution are discussed. In addition, several methodological issues and limitations associated with these analyses are highlighted emphasizing the potential to improve our understanding of the evolutionary dynamics of legume diversification by using much more densely sampled phylogenetic trees that integrate information across broad taxonomic, geographical and temporal levels.     

Cophenetic correlation analysis as a strategy to select phylogenetically informative proteins: an example from the fungal kingdom

Background  

The construction of robust and well resolved phylogenetic trees is important for our understanding of many, if not all biological processes, including speciation and origin of higher taxa, genome evolution, metabolic diversification, multicellularity, origin of life styles, pathogenicity and so on. Many older phylogenies were not well supported due to insufficient phylogenetic signal present in the single or few genes used in phylogenetic reconstructions. Importantly, single gene phylogenies were not always found to be congruent. The phylogenetic signal may, therefore, be increased by enlarging the number of genes included in phylogenetic studies. Unfortunately, concatenation of many genes does not take into consideration the evolutionary history of each individual gene. Here, we describe an approach to select informative phylogenetic proteins to be used in the Tree of Life (TOL) and barcoding projects by comparing the cophenetic correlation coefficients (CCC) among individual protein distance matrices of proteins, using the fungi as an example. The method demonstrated that the quality and number of concatenated proteins is important for a reliable estimation of TOL. Approximately 40–45 concatenated proteins seem needed to resolve fungal TOL.     

Multilocus phylogenetic analysis of the genus Aeromonas

A broad multilocus phylogenetic analysis (MLPA) of the representative diversity of a genus offers the opportunity to incorporate concatenated inter-species phylogenies into bacterial systematics. Recent analyses based on single housekeeping genes have provided coherent phylogenies of Aeromonas. However, to date, a multi-gene phylogenetic analysis has never been tackled. In the present study, the intra- and inter-species phylogenetic relationships of 115 strains representing all Aeromonas species described to date were investigated by MLPA. The study included the independent analysis of seven single gene fragments (gyrB, rpoD, recA, dnaJ, gyrA, dnaX, and atpD), and the tree resulting from the concatenated 4705 bp sequence. The phylogenies obtained were consistent with each other, and clustering agreed with the Aeromonas taxonomy recognized to date. The highest clustering robustness was found for the concatenated tree (i.e. all Aeromonas species split into 100% bootstrap clusters). Both possible chronometric distortions and poor resolution encountered when using single-gene analysis were buffered in the concatenated MLPA tree. However, reliable phylogenetic species delineation required an MLPA including several “bona fide” strains representing all described species.     

Hierarchies of evolutionary radiation in the world’s most species rich vertebrate group,the Neotropical Pristimantis leaf litter frogs

The Neotropical leaf litter frog genus Pristimantis is very species-rich, with 526 species described to date, but the full extent of its diversity is much higher and remains unknown. This study explores the phylogenetic processes and resulting evolutionary patterns of diversification in Pristimantis. Given the well-recognised failure of morphology- and community-based species groups to describe diversity within the genus, we apply a new test for the presence and phylogenetic distribution of higher evolutionary units. We developed a phylogeny based on 260 individuals encompassing 149 Pristimantis presumed species, sampled at mitochondrial and nuclear genes (3718 base pair alignment), combining new and available sequence data. Our phylogeny broadly agrees with previous studies, both in topology and age estimates, with the origin of Pristimantis at 28.97 (95% HDP =21.59 – 37.33) million years ago (MYA). New taxa that we add to the genus, which had not previously been included in Pristimantis phylogenies, suggest considerable diversity remains to be described. We assessed patterns of lineage origin and recovered 14 most likely (95% CI: 13–19) phylogenetic clusters or higher evolutionary significant units (hESUs) within Pristimantis. Diversification rates decrease towards the present following a density-dependent pattern for Pristimantis overall and for most hESU clusters, reflecting historical evolutionary radiation. The timing of diversification suggests that geological events in the Miocene, such as Andes orogenesis and Pebas system formation and drainage, may have had a direct or indirect impact on the evolution of Pristimantis and thus contributed to the origins of evolutionary independent phylogenetic clusters.     

Aliarcobacter vitoriensis sp. nov., isolated from carrot and urban wastewater

Two isolates, one recovered from a carrot and another one from urban wastewater, were characterized using a polyphasic approach. Phylogenetic analysis based on 16S rRNA gene sequences revealed that both isolates clustered together, and were most closely related to Aliarcobacter lanthieri. Multilocus phylogenetic analysis (MLPA) using the concatenated sequences of five housekeeping genes (atpA, gyrA, gyrB, hsp60 and rpoB) suggested that these isolates formed a distinct phylogenetic lineage among the genera derived from the former genus Arcobacter. Whole-genome sequence, in silico DNA-DNA hybridization (isDDH) and the average nucleotide identity (ANI) value between the genome of strain F199^T and those of related species confirmed that these isolates represent a novel species. These strains can be differentiated from its phylogenetically closest species A. lanthieri by its inability to growth on 1% glycine and by their enzyme activity of esterase lipase (C8) and acid phosphatase. Our results, by the application of a polyphasic analysis, confirmed that these two isolates represent a novel species of the genus Aliarcobacter, for which the name Aliarcobacter vitoriensis sp. nov. is proposed. The type strain is F199^T (=CECT 9230^T=LMG 30050^T).     

Correlating Molecular Phylogeny with Venom Apparatus Occurrence in Panamic Auger Snails (Terebridae)

Central to the discovery of neuroactive compounds produced by predatory marine snails of the superfamily Conoidea (cone snails, terebrids, and turrids) is identifying those species with a venom apparatus. Previous analyses of western Pacific terebrid specimens has shown that some Terebridae groups have secondarily lost their venom apparatus. In order to efficiently characterize terebrid toxins, it is essential to devise a key for identifying which species have a venom apparatus. The findings presented here integrate molecular phylogeny and the evolution of character traits to infer the presence or absence of the venom apparatus in the Terebridae. Using a combined dataset of 156 western and 33 eastern Pacific terebrid samples, a phylogenetic tree was constructed based on analyses of 16S, COI and 12S mitochondrial genes. The 33 eastern Pacific specimens analyzed represent four different species: Acus strigatus, Terebra argyosia, T. ornata, and T. cf. formosa. Anatomical analysis was congruent with molecular characters, confirming that species included in the clade Acus do not have a venom apparatus, while those in the clade Terebra do. Discovery of the association between terebrid molecular phylogeny and the occurrence of a venom apparatus provides a useful tool for effectively identifying the terebrid lineages that may be investigated for novel pharmacological active neurotoxins, enhancing conservation of this important resource, while providing supplementary information towards understanding terebrid evolutionary diversification.     

Trypanosoma cruzi: New insights on ecophylogeny and hybridization by multigene sequencing of three nuclear and one maxicircle genes

Natural populations of Trypanosoma cruzi are structured into five genetic lineages, T. cruzi I and T. cruzi II a to e, as the result of clonal evolution with rare genetic recombination events. To explore more in depth these phenomenons, a multigene sequencing approach was used, for the first time in the case of T. cruzi. Three nuclear loci and a maxicircle locus were sequenced on 18 T. cruzi stocks. Sequences were used to build phylogenetic trees from each locus and from concatenated sequences of all loci.The data confirmed the hybrid origin of DTUs IId and IIe, as the result of an ancient genetic recombination between strains pertaining to IIb and IIc. The data confirmed also a hybrid origin of DTUs IIa and IIc. Contrary to previous reports, we failed to detect mosaic genes. The phylogenetic relationship between DTUs and the respective roles of recombination and selection were tested.     

Key Processes for Cheirolophus (Asteraceae) Diversification on Oceanic Islands Inferred from AFLP Data

The radiation of the genus Cheirolophus (Asteraceae) in Macaronesia constitutes a spectacular case of rapid diversification on oceanic islands. Twenty species – nine of them included in the IUCN Red List of Threatened Species – have been described to date inhabiting the Madeiran and Canarian archipelagos. A previous phylogenetic study revealed that the diversification of Cheirolophus in Macaronesia started less than 2 Ma. As a result of such an explosive speciation process, limited phylogenetic resolution was reported, mainly due to the low variability of the employed molecular markers. In the present study, we used highly polymorphic AFLP markers to i) evaluate species'' boundaries, ii) infer their evolutionary relationships and iii) investigate the patterns of genetic diversity in relation to the potential processes likely involved in the radiation of Cheirolophus. One hundred and seventy-two individuals representing all Macaronesian Cheirolophus species were analysed using 249 AFLP loci. Our results suggest that geographic isolation played an important role in this radiation process. This was likely driven by the combination of poor gene flow capacity and a good ability for sporadic long-distance colonisations. In addition, we also found some traces of introgression and incipient ecological adaptation, which could have further enhanced the extraordinary diversification of Cheirolophus in Macaronesia. Last, we hypothesize that current threat categories assigned to Macaronesian Cheirolophus species do not reflect their respective evolutionary relevance, so future evaluations of their conservation status should take into account the results presented here.     

The Complete Mitochondrial Genome of Aix galericulata and Tadorna ferruginea: Bearings on Their Phylogenetic Position in the Anseriformes

Aix galericulata and Tadorna ferruginea are two Anatidae species representing different taxonomic groups of Anseriformes. We used a PCR-based method to determine the complete mtDNAs of both species, and estimated phylogenetic trees based on the complete mtDNA alignment of these and 14 other Anseriforme species, to clarify Anseriform phylogenetics. Phylogenetic trees were also estimated using a multiple sequence alignment of three mitochondrial genes (Cyt b, ND2, and COI) from 68 typical species in GenBank, to further clarify the phylogenetic relationships of several groups among the Anseriformes. The new mtDNAs are circular molecules, 16,651 bp (Aix galericulata) and 16,639 bp (Tadorna ferruginea) in length, containing the 37 typical genes, with an identical gene order and arrangement as those of other Anseriformes. Comparing the protein-coding genes among the mtDNAs of 16 Anseriforme species, ATG is generally the start codon, TAA is the most frequent stop codon, one of three, TAA, TAG, and T-, commonly observed. All tRNAs could be folded into canonical cloverleaf secondary structures except for tRNASer (AGY) and tRNALeu (CUN), which are missing the "DHU" arm.Phylogenetic relationships demonstrate that Aix galericula and Tadorna ferruginea are in the same group, the Tadorninae lineage, based on our analyses of complete mtDNAs and combined gene data. Molecular phylogenetic analysis suggests the 68 species of Anseriform birds be divided into three families: Anhimidae, Anatidae, and Anseranatidae. The results suggest Anatidae birds be divided into five subfamilies: Anatinae, Tadorninae, Anserinae, Oxyurinae, and Dendrocygninae. Oxyurinae and Dendrocygninae should not belong to Anserinae, but rather represent independent subfamilies. The Anatinae includes species from the tribes Mergini, Somaterini, Anatini, and Aythyini. The Anserinae includes species from the tribes Anserini and Cygnini.     

Molecular markers based on LTR retrotransposons BARE-1 and Jeli uncover different strata of evolutionary relationships in diploid wheats

Molecular markers based on retrotransposon insertions are widely used for various applications including phylogenetic analysis. Multiple cases were described where retrotransposon-based markers, namely sequence-specific amplification polymorphism (SSAP), were superior to other marker types in resolving the phylogenetic relationships due to their higher variability and informativeness. However, the patterns of evolutionary relationships revealed by SSAP may be dependent on the underlying retrotransposon activity in different periods of time. Hence, the proper choice of retrotransposon family is essential for obtaining significant results. We compared the phylogenetic trees for a diverse set of diploid A-genome wheat species (Triticum boeoticum, T. urartu and T. monococcum) based on two unrelated retrotransposon families, BARE-1 and Jeli. BARE-1 belongs to Copia class and has a uniform distribution between common wheat (T. aestivum) genomes of different origin (A, B and D), indicating similar activity in the respective diploid genome donors. Gypsy-class family Jeli was found by us to be an A-genome retrotransposon with >70% copies residing in A genome of hexaploid common wheat, suggesting a burst of transposition in the history of A-genome progenitors. The results indicate that a higher Jeli transpositional activity was associated with T. urartu versus T. boeoticum speciation, while BARE-1 produced more polymorphic insertions during subsequent intraspecific diversification; as an outcome, each retrotransposon provides more informative markers at the corresponding level of phylogenetic relationships. We conclude that multiple retroelement families should be analyzed for an image of evolutionary relationships to be solid and comprehensive.     

New Pseudomonas spp. Are Pathogenic to Citrus

Five putative novel Pseudomonas species shown to be pathogenic to citrus have been characterized in a screening of 126 Pseudomonas strains isolated from diseased citrus leaves and stems in northern Iran. The 126 strains were studied using a polyphasic approach that included phenotypic characterizations and phylogenetic multilocus sequence analysis. The pathogenicity of these strains against 3 cultivars of citrus is demonstrated in greenhouse and field studies. The strains were initially grouped phenotypically and by their partial rpoD gene sequences into 11 coherent groups in the Pseudomonas fluorescens phylogenetic lineage. Fifty-three strains that are representatives of the 11 groups were selected and analyzed by partial sequencing of their 16S rRNA and gyrB genes. The individual and concatenated partial sequences of the three genes were used to construct the corresponding phylogenetic trees. The majority of the strains were identified at the species level: P. lurida (5 strains), P. monteilii (2 strains), P. moraviensis (1 strain), P. orientalis (16 strains), P. simiae (7 strains), P. syringae (46 strains, distributed phylogenetically in at least 5 pathovars), and P. viridiflava (2 strains). This is the first report of pathogenicity on citrus of P. orientalis, P. simiae, P. lurida, P. moraviensis and P. monteilii strains. The remaining 47 strains that could not be identified at the species level are considered representatives of at least 5 putative novel Pseudomonas species that are not yet described.     

Inferring species phylogenies from multiple genes: concatenated sequence tree versus consensus gene tree

Phylogenetic trees from multiple genes can be obtained in two fundamentally different ways. In one, gene sequences are concatenated into a super-gene alignment, which is then analyzed to generate the species tree. In the other, phylogenies are inferred separately from each gene, and a consensus of these gene phylogenies is used to represent the species tree. Here, we have compared these two approaches by means of computer simulation, using 448 parameter sets, including evolutionary rate, sequence length, base composition, and transition/transversion rate bias. In these simulations, we emphasized a worst-case scenario analysis in which 100 replicate datasets for each evolutionary parameter set (gene) were generated, and the replicate dataset that produced a tree topology showing the largest number of phylogenetic errors was selected to represent that parameter set. Both randomly selected and worst-case replicates were utilized to compare the consensus and concatenation approaches primarily using the neighbor-joining (NJ) method. We find that the concatenation approach yields more accurate trees, even when the sequences concatenated have evolved with very different substitution patterns and no attempts are made to accommodate these differences while inferring phylogenies. These results appear to hold true for parsimony and likelihood methods as well. The concatenation approach shows >95% accuracy with only 10 genes. However, this gain in accuracy is sometimes accompanied by reinforcement of certain systematic biases, resulting in spuriously high bootstrap support for incorrect partitions, whether we employ site, gene, or a combined bootstrap resampling approach. Therefore, it will be prudent to report the number of individual genes supporting an inferred clade in the concatenated sequence tree, in addition to the bootstrap support.