Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstruction

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

Trichoderma koningiopsis (Hypocreaceae) has the smallest mitogenome of the genus Trichoderma

Introduction

Fungal mitogenomes exhibit remarkable variation in conformation, size, gene content, arrangement and expression, including their intergenic spacers and introns.

Methods

The complete mitochondrial genome sequence of the mycoparasitic fungus Trichoderma koningiopsis was determined using the Illumina next-generation sequencing technology. We used data from our recent Illumina NGS-based project of T. koningiopsis genome sequencing to study its mitochondrial genome. The mitogenome was assembled, annotated, and compared with other fungal mitogenomes.

Results

T. koningiopsis strain POS7 mitogenome is a circular molecule of 27,560 bp long with a GC content of 27.80%. It harbors the whole complement of the 14 conserved mitochondrial protein-coding genes (PCG) such as atp6, atp8, atp9, cox1, cox2, cox3, cob, nad1, nad2, nad3, nad4, nad4L, nad5, and nad6, also found in the same gene order to other Hypocreales. The mitogenome also contains 26 transfer RNA genes (tRNAs), 5 of them with more than one copy. Other genes also present in the assembled mitochondrial genome are a small rRNA subunit and a large rRNA subunit containing ribosomal protein S3 gene. Despite the small genome size, two introns were detected in the T. koningiopsis POS7 mitogenome, one of them in cox3 gene and the other in rnl gene, accounting 7.34% of this mitogenome with a total size of 2,024 bp. A phylogenetic analysis was done using the 14 PCGs genes of T. koningiopsis strain POS7 mitogenome to compare them with those from other fungi of the Subphyla Pezizomycotina and Saccharomycotina. T. koningiopsis strain POS7 was clustered together with other representatives of Trichoderma lineage, within the Hypocreales group, which is also supported by previous phylogenetic studies based on nuclear markers.

Discussion

The mitochondrial genome of T. koningiopsis POS7 will allow further investigations into the taxonomy, phylogenetics, conservation genetics, and evolutionary biology of this important genus as well as other closely related species.

A real data-driven simulation strategy to select an imputation method for mixed-type trait data

Missing observations in trait datasets pose an obstacle for analyses in myriad biological disciplines. Considering the mixed results of imputation, the wide variety of available methods, and the varied structure of real trait datasets, a framework for selecting a suitable imputation method is advantageous. We invoked a real data-driven simulation strategy to select an imputation method for a given mixed-type (categorical, count, continuous) target dataset. Candidate methods included mean/mode imputation, k-nearest neighbour, random forests, and multivariate imputation by chained equations (MICE). Using a trait dataset of squamates (lizards and amphisbaenians; order: Squamata) as a target dataset, a complete-case dataset consisting of species with nearly complete information was formed for the imputation method selection. Missing data were induced by removing values from this dataset under different missingness mechanisms: missing completely at random (MCAR), missing at random (MAR), and missing not at random (MNAR). For each method, combinations with and without phylogenetic information from single gene (nuclear and mitochondrial) or multigene trees were used to impute the missing values for five numerical and two categorical traits. The performances of the methods were evaluated under each missing mechanism by determining the mean squared error and proportion falsely classified rates for numerical and categorical traits, respectively. A random forest method supplemented with a nuclear-derived phylogeny resulted in the lowest error rates for the majority of traits, and this method was used to impute missing values in the original dataset. Data with imputed values better reflected the characteristics and distributions of the original data compared to complete-case data. However, caution should be taken when imputing trait data as phylogeny did not always improve performance for every trait and in every scenario. Ultimately, these results support the use of a real data-driven simulation strategy for selecting a suitable imputation method for a given mixed-type trait dataset.

Mitogenomic Comparison of the Mole Crickets Gryllotalpidae with the Phylogenetic Implications (Orthoptera: Ensifera)

Owing to limited molecular data, the phylogenetic position of the family Gryllotalpidae is still controversial in the infraorder Gryllidea. Mitochondrial genome (mitogenome) plays a crucial role in reconstructing phylogenetic relationships and revealing the molecular evolution of insects. However, only four mitogenomes have been reported in Gryllotalpidae to date. Herein, we obtained the first mitogenomes of Gryllotalpa henana Cai & Niu, 1998 and the Chinese G. orientalis Burmeister, 1838, made a detailed comparison of all mitogenomes available in Gryllotalpidae and reconstructed the phylogeny of Gryllidea based on mitogenomes using Bayesian inference (BI) and maximum likelihood (ML) methods. The results show that the complete mitogenome sequences of G. henana (15,504 bp) and G. orientalis (15,497 bp) are conserved, both exhibiting the double-stranded circular structure, typical gene content and the ancestral insect gene arrangement. The complete mitogenome of G.henana exhibits the lowest average AT content ever detected in Gryllotalpidae, and even Gryllidea. The gene nad2 of both species has atypical initiation codon GTG. All tRNAs exhibit typical clover-leaf structure, except for trnS1 lacking the dihydrouridine (DHU) arm. A potential stem-loop structure, containing a (T)_n(TC)₂(T)_n sequence, is detected in the control region of all gryllotalpids investigated and is likely related to the replication initiation of the minority strand. The phylogenetic analyses recover the six families of Gryllidea as Gryllotalpidae + (Myrmecophilidae + (Mogoplistidae + (Trigonidiidae + (Phalangopsidae + Gryllidae)))), similar to the trees based on transcriptomic and mitogenomic data. However, the trees are slightly different from the multilocus phylogenies, which show the sister-group relationship of Gryllotalpidae and Myrmecophilidae. The contradictions between mitogenomic and multilocus trees are briefly discussed.

Revised phylogeny from complete mitochondrial genomes of phyllostomid bats resolves subfamilial classification

Classically, molecular phylogenetic trees of Phyllostomidae have been inferred using a combination of a few mitochondrial and nuclear markers. However, there is still uncertainty in the relationships, especially among deep clades within the family. In this study, we provide newly sequenced complete mitochondrial genomes from 26 bat species, including genomes of 23 species reported here for the first time. By carefully analysing these genomes using maximum likelihood and Bayesian methods and different ingroup and outgroup samples, partition schemes and data types, we investigated the robustness and sensitivity of our phylogenetic results. The optimal topologies were those inferred from the complete data matrix of nucleotides, with complex and highly parameterized substitution models and partition schemes. Our results show a statistically robust picture of the evolutionary relationships between phyllostomid subfamilies and clarify hitherto uncertain relationships of Lonchorhininae and Macrotinae.

Rapid divergence of a gamete recognition gene promoted macroevolution of Eutheria

BACKGROUND

Speciation genes contribute disproportionately to species divergence, but few examples exist, especially in vertebrates. Here we test whether Zan, which encodes the sperm acrosomal protein zonadhesin that mediates species-specific adhesion to the egg's zona pellucida, is a speciation gene in placental mammals.

RESULTS

Genomic ontogeny reveals that Zan arose by repurposing of a stem vertebrate gene that was lost in multiple lineages but retained in Eutheria on acquiring a function in egg recognition. A 112-species Zan sequence phylogeny, representing 17 of 19 placental Orders, resolves all species into monophyletic groups corresponding to recognized Orders and Suborders, with <5% unsupported nodes. Three other rapidly evolving germ cell genes (Adam2, Zp2, and Prm1), a paralogous somatic cell gene (TectA), and a mitochondrial gene commonly used for phylogenetic analyses (Cytb) all yield trees with poorer resolution than the Zan tree and inferior topologies relative to a widely accepted mammalian supertree. Zan divergence by intense positive selection produces dramatic species differences in the protein's properties, with ordinal divergence rates generally reflecting species richness of placental Orders consistent with expectations for a speciation gene that acts across a wide range of taxa. Furthermore, Zan's combined phylogenetic utility and divergence exceeds those of all other genes known to have evolved in Eutheria by positive selection, including the only other mammalian speciation gene, Prdm9.

CONCLUSIONS

Species-specific egg recognition conferred by Zan's functional divergence served as a mode of prezygotic reproductive isolation that promoted the extraordinary adaptive radiation and success of Eutheria.

Draft Genome Assembly of an Iconic Arctic Species: Muskox (Ovibos moschatus)

Muskoxen (Ovibos moschatus) are Arctic species within the Caprinae subfamily that are economically and culturally significant to northern Indigenous communities. Low genetic diversity from repeated genetic bottlenecks, coupled with the effects of Arctic warming (e.g., heat stress, changing forage, pathogen range expansions), present conservation concerns for this species. Reference genome assemblies enhance our ecological and evolutionary understanding of species (which in turn aid conservation efforts). Herein, we provide a full draft reference genome of muskox using Illumina Hiseq data and cross-species scaffolding. The final reference assembly yielded a genome of 2,621,890,883 bp in length, a scaffold N50 of ~13.2 million, and an annotation identifying ~19.3 k genes. The muskox genome assembly and annotation were then used to reconstruct a phylogenetic tree which estimated muskoxen diverged from other ungulate species~12 Mya. To gain insight into the demographic history of muskoxen we also performed pairwise sequentially Markovian coalescent (PSMC) that identified two population bottlenecks coinciding with major glaciation events contributing to the notoriously low genetic variation observed in muskoxen. Overall, this genome assembly provides a foundation for future population genomic studies, such as latitudinal analyses, to explore the capacity of muskoxen to adapt to rapidly changing environments.

Phylogenetic Signal Dissection of Heterogeneous 28S and 16S rRNA Genes in Spinicaudata (Branchiopoda, Diplostraca)

It is still a challenge to reconstruct the deep phylogenetic relationships within spinicaudatans, and there are several different competing hypotheses regarding the interrelationships among Eocyzicidae, Cyzicidae s. s., Leptestheriidae, and Limnadiidae of the Suborder Spinicaudata. In order to explore the source of the inconsistencies, we focus on the sequence variation and the structure model of two rRNA genes based on extensive taxa sampling. The comparative sequence analysis revealed heterogeneity across species and the existence of conserved motifs in all spinicaudatan species. The level of intraspecific heterogeneity differed among species, which suggested that some species might have undergone a relaxed concerted evolution with respect to the 28S rRNA gene. The Bayesian analyses were performed on nuclear (28S rRNA, EF1α) and mitochondrial (16S rRNA, COI) genes. Further, we investigated compositional heterogeneity between lineages and assessed the potential for phylogenetic noise compared to signal in the combined data set. Reducing the non-phylogenetic signals and application of optimal rRNA model recovered a topology congruent with inference from the transcriptome data, whereby Limnadiidae was placed as a sister group to Leptestheriidae + Eocyzicidae with high support (topology I). Tests of alternative hypotheses provided implicit support for four competing topologies, and topology I was the best.

Lower statistical support with larger datasets: insights from the Ochrophyta radiation

It is commonly assumed that increasing the number of characters has the potential to resolve evolutionary radiations. Here, we studied photosynthetic stramenopiles (Ochrophyta) using alignments of heterogeneous origin mitochondrion, plastid, and nucleus. Surprisingly while statistical support for the relationships between the six major Ochrophyta lineages increases when comparing the mitochondrion (6,762 sites) and plastid (21,692 sites) trees, it decreases in the nuclear (209,105 sites) tree. Statistical support is not simply related to the data set size but also to the quantity of phylogenetic signal available at each position and our ability to extract it. Here, we show that this ability for current phylogenetic methods is limited, because conflicting results were obtained when varying taxon sampling. Even though the use of a better fitting model improved signal extraction and reduced the observed conflicts, the plastid data set provided higher statistical support for the ochrophyte radiation than the larger nucleus data set. We propose that the higher support observed in the plastid tree is due to an acceleration of the evolutionary rate in one short deep internal branch, implying that more phylogenetic signal per position is available to resolve the Ochrophyta radiation in the plastid than in the nuclear data set. Our work therefore suggests that, in order to resolve radiations, beyond the obvious use of data sets with more positions, we need to continue developing models of sequence evolution that better extract the phylogenetic signal and design methods to search for genes/characters that contain more signal specifically for short internal branches.

Autosomal, sex-linked and mitochondrial loci resolve evolutionary relationships among wrens in the genus Campylorhynchus

Although there is general consensus that sampling of multiple genetic loci is critical in accurate reconstruction of species trees, the exact numbers and the best types of molecular markers remain an open question. In particular, the phylogenetic utility of sex-linked loci is underexplored. Here, we sample all species and 70% of the named diversity of the New World wren genus Campylorhynchus using sequences from 23 loci, to evaluate the effects of linkage on efficiency in recovering a well-supported tree for the group. At a tree-wide level, we found that most loci supported fewer than half the possible clades and that sex-linked loci produced similar resolution to slower-coalescing autosomal markers, controlling for locus length. By contrast, we did find evidence that linkage affected the efficiency of recovery of individual relationships; as few as two sex-linked loci were necessary to resolve a selection of clades with long to medium subtending branches, whereas 4-6 autosomal loci were necessary to achieve comparable results. These results support an expanded role for sampling of the avian Z chromosome in phylogenetic studies, including target enrichment approaches. Our concatenated and species tree analyses represent significant improvements in our understanding of diversification in Campylorhynchus, and suggest a relatively complex scenario for its radiation across the Miocene/Pliocene boundary, with multiple invasions of South America.

A chromosome-level assembly of the black tiger shrimp (Penaeus monodon) genome facilitates the identification of growth-associated genes

To salvage marine ecosystems from fishery overexploitation, sustainable and efficient aquaculture must be emphasized. The knowledge obtained from available genome sequence of marine organisms has accelerated marine aquaculture in many cases. The black tiger shrimp (Penaeus monodon) is one of the most prominent cultured penaeid shrimps (Crustacean) with an average annual global production of half a million tons in the last decade. However, its currently available genome assemblies lack the contiguity and completeness required for accurate genome annotation due to the highly repetitive nature of the genome and technical difficulty in extracting high-quality, high-molecular weight DNA. Here, we report the first chromosome-level whole-genome assembly of P. monodon. The combination of long-read Pacific Biosciences (PacBio) and long-range Chicago and Hi-C technologies enabled a successful assembly of this first high-quality genome sequence. The final assembly covered 2.39 Gb (92.3% of the estimated genome size) and contained 44 pseudomolecules, corresponding to the haploid chromosome number. Repetitive elements occupied a substantial portion of the assembly (62.5%), the highest of the figures reported among crustacean species. The availability of this high-quality genome assembly enabled the identification of genes associated with rapid growth in the black tiger shrimp through the comparison of hepatopancreas transcriptome of slow-growing and fast-growing shrimps. The results highlighted several growth-associated genes. Our high-quality genome assembly provides an invaluable resource for genetic improvement and breeding penaeid shrimp in aquaculture. The availability of P. monodon genome enables analyses of ecological impact, environment adaptation and evolution, as well as the role of the genome to protect the ecological resources by promoting sustainable shrimp farming.

Historical biogeography, systematics, and integrative taxonomy of the non-Ethiopian speckled pelage brush-furred rats (Lophuromys flavopunctatus group)

BACKGROUND

The speckled-pelage brush-furred rats (Lophuromys flavopunctatus group) have been difficult to define given conflicting genetic, morphological, and distributional records that combine to obscure meaningful accounts of its taxonomic diversity and evolution. In this study, we inferred the systematics, phylogeography, and evolutionary history of the L. flavopunctatus group using maximum likelihood and Bayesian phylogenetic inference, divergence times, historical biogeographic reconstruction, and morphometric discriminant tests. We compiled comprehensive datasets of three loci (two mitochondrial [mtDNA] and one nuclear) and two morphometric datasets (linear and geometric) from across the known range of the genus Lophuromys.

RESULTS

The mtDNA phylogeny supported the division of the genus Lophuromys into three primary groups with nearly equidistant pairwise differentiation: one group corresponding to the subgenus Kivumys (Kivumys group) and two groups corresponding to the subgenus Lophuromys (L. sikapusi group and L. flavopunctatus group). The L. flavopunctatus group comprised the speckled-pelage brush-furred Lophuromys endemic to Ethiopia (Ethiopian L. flavopunctatus members [ETHFLAVO]) and the non-Ethiopian ones (non-Ethiopian L. flavopunctatus members [NONETHFLAVO]) in deeply nested relationships. There were distinctly geographically structured mtDNA clades among the NONETHFLAVO, which were incongruous with the nuclear tree where several clades were unresolved. The morphometric datasets did not systematically assign samples to meaningful taxonomic units or agree with the mtDNA clades. The divergence dating and ancestral range reconstructions showed the NONETHFLAVO colonized the current ranges over two independent dispersal events out of Ethiopia in the early Pleistocene.

CONCLUSION

The phylogenetic associations and divergence times of the L. flavopunctatus group support the hypothesis that paleoclimatic impacts and ecosystem refugia during the Pleistocene impacted the evolutionary radiation of these rodents. The overlap in craniodental variation between distinct mtDNA clades among the NONETHFLAVO suggests unraveling underlying ecomorphological drivers is key to reconciling taxonomically informative morphological characters. The genus Lophuromys requires a taxonomic reassessment based on extensive genomic evidence to elucidate the patterns and impacts of genetic isolation at clade contact zones.

Nuclear phylotranscriptomics and phylogenomics support numerous polyploidization events and hypotheses for the evolution of rhizobial nitrogen-fixing symbiosis in Fabaceae

Fabaceae are the third largest angiosperm family, with 765 genera and ∼19 500 species. They are important both economically and ecologically, and global Fabaceae crops are intensively studied in part for their nitrogen-fixing ability. However, resolution of the intrasubfamilial Fabaceae phylogeny and divergence times has remained elusive, precluding a reconstruction of the evolutionary history of symbiotic nitrogen fixation in Fabaceae. Here, we report a highly resolved phylogeny using >1500 nuclear genes from newly sequenced transcriptomes and genomes of 391 species, along with other datasets, for a total of 463 legumes spanning all 6 subfamilies and 333 of 765 genera. The subfamilies are maximally supported as monophyletic. The clade comprising subfamilies Cercidoideae and Detarioideae is sister to the remaining legumes, and Duparquetioideae and Dialioideae are successive sisters to the clade of Papilionoideae and Caesalpinioideae. Molecular clock estimation revealed an early radiation of subfamilies near the K/Pg boundary, marked by mass extinction, and subsequent divergence of most tribe-level clades within ∼15 million years. Phylogenomic analyses of thousands of gene families support 28 proposed putative whole-genome duplication/whole-genome triplication events across Fabaceae, including those at the ancestors of Fabaceae and five of the subfamilies, and further analyses supported the Fabaceae ancestral polyploidy. The evolution of rhizobial nitrogen-fixing nodulation in Fabaceae was probed by ancestral character reconstruction and phylogenetic analyses of related gene families and the results support the hypotheses of one or two switch(es) to rhizobial nodulation followed by multiple losses. Collectively, these results provide a foundation for further morphological and functional evolutionary analyses across Fabaceae.

Asterid Phylogenomics/Phylotranscriptomics Uncover Morphological Evolutionary Histories and Support Phylogenetic Placement for Numerous Whole-Genome Duplications

Asterids are one of the most successful angiosperm lineages, exhibiting extensive morphological diversity and including a number of important crops. Despite their biological prominence and value to humans, the deep asterid phylogeny has not been fully resolved, and the evolutionary landscape underlying their radiation remains unknown. To resolve the asterid phylogeny, we sequenced 213 transcriptomes/genomes and combined them with other data sets, representing all accepted orders and nearly all families of asterids. We show fully supported monophyly of asterids, Berberidopsidales as sister to asterids, monophyly of all orders except Icacinales, Aquifoliales, and Bruniales, and monophyly of all families except Icacinaceae and Ehretiaceae. Novel taxon placements benefited from the expanded sampling with living collections from botanical gardens, resolving hitherto uncertain relationships. The remaining ambiguous placements here are likely due to limited sampling and could be addressed in the future with relevant additional taxa. Using our well-resolved phylogeny as reference, divergence time estimates support an Aptian (Early Cretaceous) origin of asterids and the origin of all orders before the Cretaceous-Paleogene boundary. Ancestral state reconstruction at the family level suggests that the asterid ancestor was a woody terrestrial plant with simple leaves, bisexual, and actinomorphic flowers with free petals and free anthers, a superior ovary with a style, and drupaceous fruits. Whole-genome duplication (WGD) analyses provide strong evidence for 33 WGDs in asterids and one in Berberidopsidales, including four suprafamilial and seven familial/subfamilial WGDs. Our results advance the understanding of asterid phylogeny and provide numerous novel evolutionary insights into their diversification and morphological evolution.

The number of species of degus (genus Octodon) is currently underestimated: An appraisal of species limits and their phylogenetic relationships (Rodentia: Hystricomorpha: Octodontidae)

As currently understood, the genus Octodon contains five species degus, lunatus, bridgesii, pacificus, and ricardojeda. Previous phylogenetic studies suggest that genus specific diversity is underestimated. In order to evaluate the taxonomic diversity of Octodon, we implemented unilocus (cytochrome-b) and multilocus (cytochrome-b + 4 nuclear genes) species delimitation methods. Octodon degus was recovered as a sister of the other species of the genus. The unilocus bGMYC and mPTP methods, based on cytochrome-b sequences, delimits 11 and 7 candidate species respectively, and both methods fail to recognize O. pacificus from O. ricardojeda. Results of the multilocus analysis (BPP) vary as a function of the dataset used. When the five genes are used 11 species are delimited, while eight species are delimited when only the nuclear genes are used. Octodon bridgesii is shown as comprising at least two species (one on the Pacific coast and the typical form found on the Andean slopes), while O. ricardojeda may comprise two species (one on the Chilean side of the Andes and the other in Argentina). Likewise, both multilocus matrices recover O. pacificus as a distinct species. This shows that species diversity of Octodon is underestimated. Remarkably, many of the delimited species based on genetic data are morphologically differentiated in cranio-dental characteristics. However, a pair of species has not achieved morphological differentiation, being cryptic species. Finally, the incongruence between mitochondrial and nuclear phylogenies suggests that processes such as incomplete lineage sorting and/or introgression have been present during the radiation of the genus.

Enhancing mitogenomic phylogeny and resolving the relationships of extinct megafaunal placental mammals

Mitochondrial genomes provided the first widely used sequences that were sufficiently informative to resolve relationships among animals across a wide taxonomic domain, from within species to between phyla. However, mitogenome studies supported several anomalous relationships and fell partly out of favour as sequencing multiple, independent nuclear loci proved to be highly effective. A tendency to blame mitochondrial DNA (mtDNA) has overshadowed efforts to understand and ameliorate underlying model misspecification. Here we find that influential assessments of the infidelity of mitogenome phylogenies have often been overstated, but nevertheless, substitution saturation and compositional non-stationarity substantially mislead reconstruction. We show that RY coding the mtDNA, excluding protein-coding 3rd codon sites, partitioning models based on amino acid hydrophobicity and enhanced taxon sampling improve the accuracy of mitogenomic phylogeny reconstruction for placental mammals, almost to the level of multi-gene nuclear datasets. Indeed, combined analysis of mtDNA with 3-fold longer nuclear sequence data either maintained or improved upon the nuclear support for all generally accepted clades, even those that mtDNA alone did not favour, thus indicating "hidden support". Confident mtDNA phylogeny reconstruction is especially important for understanding the evolutionary dynamics of mitochondria themselves, and for merging extinct taxa into the tree of life, with ancient DNA often only accessible as mtDNA. Our ancient mtDNA analyses lend confidence to the relationships of three extinct megafaunal taxa: glyptodonts are nested within armadillos, the South American ungulate, Macrauchenia is sister to horses and rhinoceroses, and sabre-toothed and scimitar cats are the monophyletic sister-group of modern cats.

Structure, gene order, and nucleotide composition of mitochondrial genomes in parasitic lice from Amblycera

Parasitic lice have unique mitochondrial (mt) genomes characterized by rearranged gene orders, variable genome structures, and less AT content compared to most other insects. However, relatively little is known about the mt genomes of Amblycera, the suborder sister to all other parasitic lice. Comparing among nine different genera (including representative of all seven families), we show that Amblycera have variable and highly rearranged mt genomes. Some genera have fragmented genomes that vary considerably in length, whereas others have a single mt chromosome. Notably, these genomes are more AT-biased than most other lice. We also recover genus-level phylogenetic relationships among Amblycera that are consistent with those reported from large nuclear datasets, indicating that mt sequences are reliable for reconstructing evolutionary relationships in Amblycera. However, gene order data cannot reliably recover these same relationships. Overall, our results suggest that the mt genomes of lice, already know to be distinctive, are even more variable than previously thought.

Molecular evidence of hybridization between pig and human Ascaris indicates an interbred species complex infecting humans

Human ascariasis is a major neglected tropical disease caused by the nematode Ascaris lumbricoides. We report a 296 megabase (Mb) reference-quality genome comprised of 17,902 protein-coding genes derived from a single, representative Ascaris worm. An additional 68 worms were collected from 60 human hosts in Kenyan villages where pig husbandry is rare. Notably, the majority of these worms (63/68) possessed mitochondrial genomes that clustered closer to the pig parasite Ascaris suum than to A. lumbricoides. Comparative phylogenomic analyses identified over 11 million nuclear-encoded SNPs but just two distinct genetic types that had recombined across the genomes analyzed. The nuclear genomes had extensive heterozygosity, and all samples existed as genetic mosaics with either A. suum-like or A. lumbricoides-like inheritance patterns supporting a highly interbred Ascaris species genetic complex. As no barriers appear to exist for anthroponotic transmission of these 'hybrid' worms, a one-health approach to control the spread of human ascariasis will be necessary.

Digging for the spiny rat and hutia phylogeny using a gene capture approach, with the description of a new mammal subfamily

Next generation sequencing (NGS) and genomic database mining allow biologists to gather and select large molecular datasets well suited to address phylogenomics and molecular evolution questions. Here we applied this approach to a mammal family, the Echimyidae, for which generic relationships have been difficult to recover and often referred to as a star phylogeny. These South-American spiny rats represent a family of caviomorph rodents exhibiting a striking diversity of species and life history traits. Using a NGS exon capture protocol, we isolated and sequenced ca. 500 nuclear DNA exons for 35 species belonging to all major echimyid and capromyid clades. Exons were carefully selected to encompass as much diversity as possible in terms of rate of evolution, heterogeneity in the distribution of site-variation and nucleotide composition. Supermatrix inferences and coalescence-based approaches were subsequently applied to infer this family's phylogeny. The inferred topologies were the same for both approaches, and support was maximal for each node, entirely resolving the ambiguous relationships of previous analyses. Fast-evolving nuclear exons tended to yield more reliable phylogenies, as slower-evolving sequences were not informative enough to disentangle the short branches of the Echimyidae radiation. Based on this resolved phylogeny and on molecular and morphological evidence, we confirm the rank of the Caribbean hutias - formerly placed in the Capromyidae family - as Capromyinae, a clade nested within Echimyidae. We also name and define Carterodontinae, a new subfamily of Echimyidae, comprising the extant monotypic genus Carterodon from Brazil, which is the closest living relative of West Indies Capromyinae.

Complete mitochondrial genome of bamboo grasshopper, Ceracris fasciata, and the phylogenetic analyses and divergence time estimation of Caelifera (Orthoptera)

The bamboo grasshopper Ceracris fasciata is regarded as a major pest species because of the damage it causes to bamboo, and its classification within the families and subfamilies of the suborder Caelifera remains unclear. Thus, we attempted to resolve these questions using molecular biology methods and analyses. Our results are as follows: (1) the complete mitochondrial genome of C. fasciata is 15,569 bp in length. The mitochondrial genome contains a standard set of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and an A + T-rich region in the same order as those of the other analysed Caeliferan species. The putative start codon for the COX1 gene in C. fasciata is ACC, although it is not defined in other genes. The presence of tandem repeats of different sizes in the A + T-rich region may lead to size differences in other mitochondrial genomes. The mitochondrial genome of C. fasciata harbours the typical 37 genes and an A + T-rich region, and it shows similar characteristics to those of other grasshopper species. Characterization of the mitochondrial genome has enriched our knowledge of the mitochondrial genomes of Orthoptera around the world. Therefore, the phylogenetic relationships in Orthoptera can be re-examined. (2) In phylogenetic analyses, the monophyly of Orthoptera and its two suborders (Caelifera and Ensifera) has been consistently recovered based on most of the datasets selected, regardless of the optimal criteria. Our results do not support the monophyly of the subfamily Oedipodinae of Caelifera. We found that Phlaeoba albonema of the Acridinae is sorted into a clade with Ceracris in all our phylogenetic trees, and field experiments show that Phlaeoba always lives with Ceracris in the same ecotopes. Therefore, we suggest that Phlaeoba should be classified as a member of the Oedipodinae. We found that C. fasciata always clustered with Ceracris kiangsu, and both were sisters to Ceracris versicolor. Therefore, the genetic relationship between C. fasciata and C. kiangsu is closer than that between C. fasciata and C. versicolor. (3) The oldest estimated time of divergence of Ensifera in this context was determined to be 146.16 million years ago (Mya), or around the late Jurassic or early Cretaceous. We estimated that katydids (Grylloidea) likely diverged from other groups in the early Cretaceous. According to our divergence time analyses, we concluded that the ancestral Acrididae probably originated in the early Paleogene, and it is likely that the major diversification events happened at the middle Paleogene, well into the next geologic time. We estimated that crickets (Tettigoniidae) likely diverged from other groups in the early Cretaceous. Acrididae and Romaleinae group, Pyrgacrididae and Ommexechidae group, the youngest two clades we observed, were estimated to have diverged 58.79 Mya, between the middle and early Paleogene. C. versicolor is a sister to the group containing C. kiangsu and C. fasciata. First, C. versicolor diverged from the sister group (C. kiangsu + C. fasciata) around 44.81 Mya, and then the C. kiangsu and C. fasciata group separated at 43.04 Mya.

Using a multi-gene approach to infer the complicated phylogeny and evolutionary history of lorises (Order Primates: Family Lorisidae)

Extensive phylogenetic studies have found robust phylogenies are modeled by using a multi-gene approach and sampling from the majority of the taxa of interest. Yet, molecular studies focused on the lorises, a cryptic primate family, have often relied on one gene, or just mitochondrial DNA, and many were unable to include all four genera in the analyses, resulting in inconclusive phylogenies. Past phylogenetic loris studies resulted in lorises being monophyletic, paraphyletic, or an unresolvable trichotomy with the closely related galagos. The purpose of our study is to improve our understanding of loris phylogeny and evolutionary history by using a multi-gene approach. We used the mitochondrial genes cytochrome b, and cytochrome c oxidase subunit 1, along with a nuclear intron (recombination activating gene 2) and nuclear exon (the melanocortin 1 receptor). Maximum Likelihood and Bayesian phylogenetic analyses were conducted based on data from each locus, as well as on the concatenated sequences. The robust, concatenated results found lorises to be a monophyletic family (Lorisidae) (PP ≥ 0.99) with two distinct subfamilies: the African Perodictinae (PP ≥ 0.99) and the Asian Lorisinae (PP ≥ 0.99). Additionally, from these analyses all four genera were all recovered as monophyletic (PP ≥ 0.99). Some of our single-gene analyses recovered monophyly, but many had discordances, with some showing paraphyly or a deep-trichotomy. Bayesian partitioned analyses inferred the most recent common ancestors of lorises emerged ∼42 ± 6 million years ago (mya), the Asian Lorisinae separated ∼30 ± 9 mya, and Perodictinae arose ∼26 ± 10 mya. These times fit well with known historical tectonic shifts of the area, as well as with the sparse loris fossil record. Additionally, our results agree with previous multi-gene studies on Lorisidae which found lorises to be monophyletic and arising ∼40 mya (Perelman et al., 2011; Pozzi et al., 2014). By taking a multi-gene approach, we were able to recover a well-supported, monophyletic loris phylogeny and inferred the evolutionary history of this cryptic family.

Genus-level phylogeny of cephalopods using molecular markers: current status and problematic areas

Comprising more than 800 extant species, the class Cephalopoda (octopuses, squid, cuttlefish, and nautiluses) is a fascinating group of marine conchiferan mollusks. Recently, the first cephalopod genome (of Octopus bimaculoides) was published, providing a genomic framework, which will enable more detailed investigations of cephalopod characteristics, including developmental, morphological, and behavioural traits. Meanwhile, a robust phylogeny of the members of the subclass Coleoidea (octopuses, squid, cuttlefishes) is crucial for comparative and evolutionary studies aiming to investigate the group's traits and innovations, but such a phylogeny has proven very challenging to obtain. Here, we present the results of phylogenetic inference at the genus level using mitochondrial and nuclear marker sequences available from public databases. Topologies are presented which show support for (1) the monophyly of the two main superorders, Octobrachia and Decabrachia, and (2) some of the interrelationships at the family level. We have mapped morphological characters onto the tree and conducted molecular dating analyses, obtaining congruent results with previous estimates of divergence in major lineages. Our study also identifies unresolved phylogenetic relationships within the cephalopod phylogeny and insufficient taxonomic sampling among squids excluding the Loliginidae in the Decabrachia and within the Order Cirromorphida in the Octobrachia. Genomic and transcriptomic resources should enable resolution of these issues in the relatively near future. We provide our alignment as an open access resource, to allow other researchers to reconstruct phylogenetic trees upon this work in the future.

Molecular phylogeny and patterns of diversification in syngnathid fishes

The family Syngnathidae is a large and diverse clade of morphologically unique bony fishes, with 57 genera and 300 described species of seahorses, pipefishes, pipehorses, and seadragons. They primarily inhabit shallow coastal waters in temperate and tropical oceans, and are characterized by a fused jaw, male brooding, and extraordinary crypsis. Phylogenetic relationships within the Syngnathidae remain poorly resolved due to lack of generic taxon sampling, few diagnostic morphological characters, and limited molecular data. The phylogenetic placement of the threatened, commercially exploited seahorses remains a topic of intense interest, with conflicting topologies based on morphology and predominantly mitochondrial genetic data. In this study, we integrate eight nuclear and mitochondrial markers and 17 morphological characters to investigate the phylogenetic structure of the family Syngnathidae at the generic level. We include 91 syngnathid species representing 48 of the 57 recognized genera, all major ocean basins, and a broad array of temperate and tropical habitats including rocky and coral reefs, sand and silt, mangroves, seagrass beds, estuaries, and rivers. Maximum likelihood and Bayesian analyses of 5160bp from eight loci produced high congruence among alternate topologies, defining well-supported and sometimes novel clades. We present a hypothesis that confirms a deep phylogenetic split between lineages with trunk- or tail-brood pouch placement, and provides significant new insights into the morphological evolution and biogeography of this highly derived fish clade. Based on the fundamental division between lineages - the tail brooding "Urophori" and the trunk brooding "Gastrophori" - we propose a revision of Syngnathidae classification into only two subfamilies: the Nerophinae and the Syngnathinae. We find support for distinct principal clades within the trunk-brooders and tail-brooders, the latter of which include seahorses, seadragons, independent lineages of pipehorses, and clades that originated in southern Australia and the Western Atlantic. We suggest the seahorse genus Hippocampus is of Indo-Pacific origin and its sister clade is an unexpected grouping of several morphologically disparate Indo-Pacific genera, including the Pacific pygmy pipehorses. Taxonomic revision is required for multiple genera, particularly to reflect deep evolutionary splits in nominal lineages from the Atlantic versus the Indo-Pacific.

Mammal madness: is the mammal tree of life not yet resolved?

Most molecular phylogenetic studies place all placental mammals into four superordinal groups, Laurasiatheria (e.g. dogs, bats, whales), Euarchontoglires (e.g. humans, rodents, colugos), Xenarthra (e.g. armadillos, anteaters) and Afrotheria (e.g. elephants, sea cows, tenrecs), and estimate that these clades last shared a common ancestor 90-110 million years ago. This phylogeny has provided a framework for numerous functional and comparative studies. Despite the high level of congruence among most molecular studies, questions still remain regarding the position and divergence time of the root of placental mammals, and certain 'hard nodes' such as the Laurasiatheria polytomy and Paenungulata that seem impossible to resolve. Here, we explore recent consensus and conflict among mammalian phylogenetic studies and explore the reasons for the remaining conflicts. The question of whether the mammal tree of life is or can be ever resolved is also addressed.This article is part of the themed issue 'Dating species divergences using rocks and clocks'.

Nuclear introns outperform mitochondrial DNA in inter-specific phylogenetic reconstruction: Lessons from horseshoe bats (Rhinolophidae: Chiroptera)

Despite many studies illustrating the perils of utilising mitochondrial DNA in phylogenetic studies, it remains one of the most widely used genetic markers for this purpose. Over the last decade, nuclear introns have been proposed as alternative markers for phylogenetic reconstruction. However, the resolution capabilities of mtDNA and nuclear introns have rarely been quantified and compared. In the current study we generated a novel ∼5kb dataset comprising six nuclear introns and a mtDNA fragment. We assessed the relative resolution capabilities of the six intronic fragments with respect to each other, when used in various combinations together, and when compared to the traditionally used mtDNA. We focused on a major clade in the horseshoe bat family (Afro-Palaearctic clade; Rhinolophidae) as our case study. This old, widely distributed and speciose group contains a high level of conserved morphology. This morphological stasis renders the reconstruction of the phylogeny of this group with traditional morphological characters complex. We sampled multiple individuals per species to represent their geographic distributions as best as possible (122 individuals, 24 species, 68 localities). We reconstructed the species phylogeny using several complementary methods (partitioned Maximum Likelihood and Bayesian and Bayesian multispecies-coalescent) and made inferences based on consensus across these methods. We computed pairwise comparisons based on Robinson-Foulds tree distance metric between all Bayesian topologies generated (27,000) for every gene(s) and visualised the tree space using multidimensional scaling (MDS) plots. Using our supported species phylogeny we estimated the ancestral state of key traits of interest within this group, e.g. echolocation peak frequency which has been implicated in speciation. Our results revealed many potential cryptic species within this group, even in taxa where this was not suspected a priori and also found evidence for mtDNA introgression. We demonstrated that by using just two introns one can recover a better supported species tree than when using the mtDNA alone, despite the shorter overall length of the combined introns. Additionally, when combining any single intron with mtDNA, we showed that the result is highly similar to the mtDNA gene tree and far from the true species tree and therefore this approach should be avoided. We caution against the indiscriminate use of mtDNA in phylogenetic studies and advocate for pilot studies to select nuclear introns. The selection of marker type and number is a crucial step that is best based on critical examination of preliminary or previously published data. Based on our findings and previous publications, we recommend the following markers to recover phylogenetic relationships between recently diverged taxa (<20 My) in bats and other mammals: ACOX2, COPS7A, BGN, ROGDI and STAT5A.

Supramolecular Organization of Respiratory Complexes

Since the discovery of the existence of superassemblies between mitochondrial respiratory complexes, such superassemblies have been the object of a passionate debate. It is accepted that respiratory supercomplexes are structures that occur in vivo, although which superstructures are naturally occurring and what could be their functional role remain open questions. The main difficulty is to make compatible the existence of superassemblies with the corpus of data that drove the field to abandon the early understanding of the physical arrangement of the mitochondrial respiratory chain as a compact physical entity (the solid model). This review provides a nonexhaustive overview of the evolution of our understanding of the structural organization of the electron transport chain from the original idea of a compact organization to a view of freely moving complexes connected by electron carriers. Today supercomplexes are viewed not as a revival of the old solid model but rather as a refined revision of the fluid model, which incorporates a new layer of structural and functional complexity.

Mitochondrial data are not suitable for resolving placental mammal phylogeny

Mitochondrial data have traditionally been used in reconstructing a variety of species phylogenies. The low rates of recombination and thorough characterization of mitochondrial data across vertebrate species make it a particularly attractive phylogenetic marker. The relatively low number of fully sequenced mammal genomes and the lack of extensive sampling within Superorders have posed a serious problem for reaching agreement on the placement mammal species. The use of mitochondrial data sequences from large numbers of mammals could serve to circumvent the taxon-sampling deficit. Here we assess the suitability of mitochondrial data as a phylogenetic marker in mammal phylogenetics. MtDNA datasets of mammal origin have been filtered as follows: (i) we have sampled sparsely across the phylogenetic tree, (ii) we have constrained our sampling to genes with high taxon coverage, (iii) we have categorised rates across sites in a phylogeny independent manner and have removed fast evolving sites, and (iv), we have sampled from very shallow divergence times to reduce phylogenetic conflict. However, topologies obtained using these filters are not consistent with previous studies and are discordant across different genes. Individual mitochondrial genes, and indeed all mitochondrial genes analysed as a supermatrix, resulted in poor resolution of the species phylogeny. Overall, our study highlights the limitations of mitochondrial data, not only for resolving deep divergences and but also for shallow divergences in the mammal phylogeny.

The phylogeny of the Orthoptera (Insecta) as deduced from mitogenomic gene sequences

Background

The phylogeny of the Orthoptera was analyzed based on 6 datasets from 47 orthopteran mitochondrial genomes (mitogenomes). The phylogenetic signals in the mitogenomes were rigorously examined under analytical regimens of maximum likelihood (ML) and Bayesian inference (BI), along with how gene types and different partitioning schemes influenced the phylogenetic reconstruction within the Orthoptera. The monophyly of the Orthoptera and its two suborders (Caelifera and Ensifera) was consistently recovered in the analyses based on most of the datasets we selected, regardless of the optimality criteria.

Results

When the seven NADH dehydrogenase subunits were concatenated into a single alignment (NADH) and were analyzed; a near-identical topology to the traditional morphological analysis was recovered, especially for BI_NADH. In both the concatenated cytochrome oxidase (COX) subunits and COX + cytochrome b (Cyt b) datasets, the small extent of sequence divergence seemed to be helpful for resolving relationships among major Orthoptera lineages (between suborders or among superfamilies). The conserved and variable domains of ribosomal (r)RNAs performed poorly when respectively analyzed but provided signals at some taxonomic levels.

Conclusions

Our findings suggest that the best phylogenetic inferences can be made when moderately divergent nucleotide data from mitogenomes are analyzed, and that the NADH dataset was suited for studying orthopteran phylogenetic relationships at different taxonomic levels, which may have been due to the larger amount of DNA sequence data and the larger number of phylogenetically informative sites.

Summary of Laurasiatheria (mammalia) phylogeny

Laurasiatheria is one of the richest and most diverse superorders of placental mammals. Because this group had a rapid evolutionary radiation, the phylogenetic relationships among the six orders of Laurasiatheria remain a subject of heated debate and several issues related to its phylogeny remain open. Reconstructing the true phylogenetic relationships of Laurasiatheria is a significant case study in evolutionary biology due to the diversity of this suborder and such research will have significant implications for biodiversity conservation. We review the higher-level (inter-ordinal) phylogenies of Laurasiatheria based on previous cytogenetic, morphological and molecular data, and discuss the controversies of its phylogenetic relationship. This review aims to outline future researches on Laurasiatheria phylogeny and adaptive evolution.

The bee tree of life: a supermatrix approach to apoid phylogeny and biogeography

BACKGROUND

Bees are the primary pollinators of angiosperms throughout the world. There are more than 16,000 described species, with broad variation in life history traits such as nesting habitat, diet, and social behavior. Despite their importance as pollinators, the evolution of bee biodiversity is understudied: relationships among the seven families of bees remain controversial, and no empirical global-level reconstruction of historical biogeography has been attempted. Morphological studies have generally suggested that the phylogeny of bees is rooted near the family Colletidae, whereas many molecular studies have suggested a root node near (or within) Melittidae. Previous molecular studies have focused on a relatively small sample of taxa (~150 species) and genes (seven at most). Public databases contain an enormous amount of DNA sequence data that has not been comprehensively analysed in the context of bee evolution.

RESULTS

We downloaded, aligned, concatenated, and analysed all available protein-coding nuclear gene DNA sequence data in GenBank as of October, 2011. Our matrix consists of 20 genes, with over 17,000 aligned nucleotide sites, for over 1,300 bee and apoid wasp species, representing over two-thirds of bee genera. Whereas the matrix is large in terms of number of genes and taxa, there is a significant amount of missing data: only ~15% of the matrix is populated with data. The placement of the root as well as relationships between Andrenidae and other bee families remain ambiguous, as several alternative maximum-likelihood estimates fall within the statistically credible set. However, we recover strong bootstrap support for relationships among many families and for their monophyly. Ancestral geographic range reconstruction suggests a likely origin of bees in the southern hemisphere, with Melittidae ancestrally located within Africa, and Halictidae, Colletidae, and Apidae within the New World.

CONCLUSIONS

Our study affirms the monophyly of each bee family, sister-taxa relationships between Apidae and Megachilidae (the 'long-tongued bees'), between Colletidae and Stenotritidae, and between Colletidae + Stenotritidae and Halictidae. Our analyses reject a Colletidae-basal hypothesis for family-level relationships and instead support Melittidae as sister to the remaining bees. Southern hemisphere vicariance likely played an important role in early diversification within many bee families.

Phylogenetic signal in the acoustic parameters of the advertisement calls of four clades of anurans

BACKGROUND

Anuran vocalizations, especially their advertisement calls, are largely species-specific and can be used to identify taxonomic affiliations. Because anurans are not vocal learners, their vocalizations are generally assumed to have a strong genetic component. This suggests that the degree of similarity between advertisement calls may be related to large-scale phylogenetic relationships. To test this hypothesis, advertisement calls from 90 species belonging to four large clades (Bufo, Hylinae, Leptodactylus, and Rana) were analyzed. Phylogenetic distances were estimated based on the DNA sequences of the 12S mitochondrial ribosomal RNA gene, and, for a subset of 49 species, on the rhodopsin gene. Mean values for five acoustic parameters (coefficient of variation of root-mean-square amplitude, dominant frequency, spectral flux, spectral irregularity, and spectral flatness) were computed for each species. We then tested for phylogenetic signal on the body-size-corrected residuals of these five parameters, using three statistical tests (Moran's I, Mantel, and Blomberg's K) and three models of genetic distance (pairwise distances, Abouheif's proximities, and the variance-covariance matrix derived from the phylogenetic tree).

RESULTS

A significant phylogenetic signal was detected for most acoustic parameters on the 12S dataset, across statistical tests and genetic distance models, both for the entire sample of 90 species and within clades in several cases. A further analysis on a subset of 49 species using genetic distances derived from rhodopsin and from 12S broadly confirmed the results obtained on the larger sample, indicating that the phylogenetic signals observed in these acoustic parameters can be detected using a variety of genetic distance models derived either from a variable mitochondrial sequence or from a conserved nuclear gene.

CONCLUSIONS

We found a robust relationship, in a large number of species, between anuran phylogenetic relatedness and acoustic similarity in the advertisement calls in a taxon with no evidence for vocal learning, even after correcting for the effect of body size. This finding, covering a broad sample of species whose vocalizations are fairly diverse, indicates that the intense selection on certain call characteristics observed in many anurans does not eliminate all acoustic indicators of relatedness. Our approach could potentially be applied to other vocal taxa.

Reconstructing the phylogenetic history of long-term effective population size and life-history traits using patterns of amino acid replacement in mitochondrial genomes of mammals and birds

The nearly neutral theory, which proposes that most mutations are deleterious or close to neutral, predicts that the ratio of nonsynonymous over synonymous substitution rates (dN/dS), and potentially also the ratio of radical over conservative amino acid replacement rates (Kr/Kc), are negatively correlated with effective population size. Previous empirical tests, using life-history traits (LHT) such as body-size or generation-time as proxies for population size, have been consistent with these predictions. This suggests that large-scale phylogenetic reconstructions of dN/dS or Kr/Kc might reveal interesting macroevolutionary patterns in the variation in effective population size among lineages. In this work, we further develop an integrative probabilistic framework for phylogenetic covariance analysis introduced previously, so as to estimate the correlation patterns between dN/dS, Kr/Kc, and three LHT, in mitochondrial genomes of birds and mammals. Kr/Kc displays stronger and more stable correlations with LHT than does dN/dS, which we interpret as a greater robustness of Kr/Kc, compared with dN/dS, the latter being confounded by the high saturation of the synonymous substitution rate in mitochondrial genomes. The correlation of Kr/Kc with LHT was robust when controlling for the potentially confounding effects of nucleotide compositional variation between taxa. The positive correlation of the mitochondrial Kr/Kc with LHT is compatible with previous reports, and with a nearly neutral interpretation, although alternative explanations are also possible. The Kr/Kc model was finally used for reconstructing life-history evolution in birds and mammals. This analysis suggests a fairly large-bodied ancestor in both groups. In birds, life-history evolution seems to have occurred mainly through size reduction in Neoavian birds, whereas in placental mammals, body mass evolution shows disparate trends across subclades. Altogether, our work represents a further step toward a more comprehensive phylogenetic reconstruction of the evolution of life-history and of the population-genetics environment.

TreeParser-aided Klee diagrams display taxonomic clusters in DNA barcode and nuclear gene datasets

Indicator vector analysis of a nucleotide sequence alignment generates a compact heat map, called a Klee diagram, with potential insight into clustering patterns in evolution. However, so far this approach has examined only mitochondrial cytochrome c oxidase I (COI) DNA barcode sequences. To further explore, we developed TreeParser, a freely-available web-based program that sorts a sequence alignment according to a phylogenetic tree generated from the dataset. We applied TreeParser to nuclear gene and COI barcode alignments from birds and butterflies. Distinct blocks in the resulting Klee diagrams corresponded to species and higher-level taxonomic divisions in both groups, and this enabled graphic comparison of phylogenetic information in nuclear and mitochondrial genes. Our results demonstrate TreeParser-aided Klee diagrams objectively display taxonomic clusters in nucleotide sequence alignments. This approach may help establish taxonomy in poorly studied groups and investigate higher-level clustering which appears widespread but not well understood.

The utility of mitochondrial and y chromosome phylogenetic data to improve correction for population stratification

Genome-wide association (GWA) studies have become a standard approach for discovering and validating genomic polymorphisms putatively associated with phenotypes of interest. Accounting for population structure in GWA studies is critical to attain unbiased parameter measurements and control Type I error. One common approach to accounting for population structure is to include several principal components derived from the entire autosomal dataset, which reflects population structure signal. However, knowing which components to include is subjective and generally not conclusive. We examined how phylogenetic signal from mitochondrial DNA (mtDNA) and chromosome Y (chr:Y) markers is concordant with principal component data based on autosomal markers to determine whether mtDNA and chr:Y phylogenetic data can help guide principal component selection. Using HAPMAP and other original data from individuals of multiple ancestries, we examined the relationships of mtDNA and chr:Y phylogenetic signal with the autosomal PCA using best subset logistic regression. We show that while the two approaches agree at times, this is independent of the component order and not completely represented in the Eigen values. Additionally, we use simulations to demonstrate that our approach leads to a slightly reduced Type I error rate compared to the standard approach. This approach provides preliminary evidence to support the theoretical concept that mtDNA and chr:Y data can be informative in locating the PCs that are most associated with evolutionary history of populations that are being studied, although the utility of such information will depend on the specific situation.

[Summary of phylogeny in family Felidae of Carnivora]

Felidae (cats) is one of the strict carnivorous groups in the order Carnivora, many of which are most familiar and spectacular to us. They are the top predators in the world. Thirty-six of 37 living cat species are considered as either "endangered" or "threatened". The relationships among species of the family Felidae, which evolved recently and rapidly, are difficult to resolve, and have been the subject of debate. Construction of a reliable Felidae phylogeny will be of evolutionarily significance and conservation value. In this paper, we summarized phylogeny of Felidae, including cytological, morphological and molecular evidence, and pointed out the existing phylogenetic problems. This review is expected to guide future researches of Felidae phylogeny, and to lay a theoretic foundation for the protection of this animal group.

The complete mitochondrial genome of the black mud crab, Scylla serrata (Crustacea: Brachyura: Portunidae) and its phylogenetic position among (pan)crustaceans

The black mud crab, Scylla serrata (Forskål 1775), is the most economically important edible crab in South-East Asia. In the present study, the complete mitochondrial genome of black mud crab, S. serrata, was determined with the sequential polymerase chain reaction and primer walking sequencing. The complete mitochondrial genome was 15,721 bp in length with an A+T content of 69.2 % and contained 37 mitochondrial genes (13 protein coding genes (PCGs), 2 ribosomal RNA genes and 22 transfer RNA genes) and a control region (CR). The analysis of the CR sequence shows that it contains a multitude of repetitive fragments which can fold into hairpin-like or secondary structures and conserved elements as in other arthropods. The gene order of S. serrata mainly retains as the pancrustacean ground pattern, except for a single translocation of trnH. The gene arrangement of S. serrata appears to be a typical feature of portunid crabs. Phylogenetic analyses with concatenated amino acid sequences of 12 PCGs establishes that S. serrata in a well-supported monophyletic Portunidae and is consistent with previous morphological classification. Moreover, the phylogenomic results strongly support monophyletic Pancrustacea (Hexapoda plus "Crustaceans"). Within Pancrustacea, this study identifies Malacostraca + Entomostraca and Branchiopoda as the sister group to Hexapoda, which confirms that "Crustacea" is not monophyletic. Cirripedia + Remipedia appear to be a basal lineage of Pancrustacea. The present study also provides considerable data for the application of both population and phylogenetic studies of other crab species.

Efficiency of nuclear and mitochondrial markers recovering and supporting known amniote groups

We have analysed the efficiency of all mitochondrial protein coding genes and six nuclear markers (Adora3, Adrb2, Bdnf, Irbp, Rag2 and Vwf) in reconstructing and statistically supporting known amniote groups (murines, rodents, primates, eutherians, metatherians, therians). The efficiencies of maximum likelihood, Bayesian inference, maximum parsimony, neighbor-joining and UPGMA were also evaluated, by assessing the number of correct and incorrect recovered groupings. In addition, we have compared support values using the conservative bootstrap test and the Bayesian posterior probabilities. First, no correlation was observed between gene size and marker efficiency in recovering or supporting correct nodes. As expected, tree-building methods performed similarly, even UPGMA that, in some cases, outperformed other most extensively used methods. Bayesian posterior probabilities tend to show much higher support values than the conservative bootstrap test, for correct and incorrect nodes. Our results also suggest that nuclear markers do not necessarily show a better performance than mitochondrial genes. The so-called dependency among mitochondrial markers was not observed comparing genome performances. Finally, the amniote groups with lowest recovery rates were therians and rodents, despite the morphological support for their monophyletic status. We suggest that, regardless of the tree-building method, a few carefully selected genes are able to unfold a detailed and robust scenario of phylogenetic hypotheses, particularly if taxon sampling is increased.

Sustained high levels of neuregulin-1 in the longest-lived rodents; a key determinant of rodent longevity

Naked mole-rats (Heterocephalus glaber), the longest-lived rodents, live 7-10 times longer than similarly sized mice and exhibit normal activities for approximately 75% of their lives. Little is known about the mechanisms that allow them to delay the aging process and live so long. Neuregulin-1 (NRG-1) signaling is critical for normal brain function during both development and adulthood. We hypothesized that long-lived species will maintain higher levels of NRG-1 and that this contributes to their sustained brain function and concomitant maintenance of normal activity. We monitored the levels of NRG-1 and its receptor ErbB4 in H. glaber at different ages ranging from 1 day to 26 years and found that levels of NRG-1 and ErbB4 were sustained throughout development and adulthood. In addition, we compared seven rodent species with widely divergent (4-32 year) maximum lifespan potential (MLSP) and found that at a physiologically equivalent age, the longer-lived rodents had higher levels of NRG-1 and ErbB4. Moreover, phylogenetic independent contrast analyses revealed that this significant strong correlation between MLSP and NRG-1 levels was independent of phylogeny. These results suggest that NRG-1 is an important factor contributing to divergent species MLSP through its role in maintaining neuronal integrity.

Statistics and truth in phylogenomics

Phylogenomics refers to the inference of historical relationships among species using genome-scale sequence data and to the use of phylogenetic analysis to infer protein function in multigene families. With rapidly decreasing sequencing costs, phylogenomics is becoming synonymous with evolutionary analysis of genome-scale and taxonomically densely sampled data sets. In phylogenetic inference applications, this translates into very large data sets that yield evolutionary and functional inferences with extremely small variances and high statistical confidence (P value). However, reports of highly significant P values are increasing even for contrasting phylogenetic hypotheses depending on the evolutionary model and inference method used, making it difficult to establish true relationships. We argue that the assessment of the robustness of results to biological factors, that may systematically mislead (bias) the outcomes of statistical estimation, will be a key to avoiding incorrect phylogenomic inferences. In fact, there is a need for increased emphasis on the magnitude of differences (effect sizes) in addition to the P values of the statistical test of the null hypothesis. On the other hand, the amount of sequence data available will likely always remain inadequate for some phylogenomic applications, for example, those involving episodic positive selection at individual codon positions and in specific lineages. Again, a focus on effect size and biological relevance, rather than the P value, may be warranted. Here, we present a theoretical overview and discuss practical aspects of the interplay between effect sizes, bias, and P values as it relates to the statistical inference of evolutionary truth in phylogenomics.

Dynamic evolution of base composition: causes and consequences in avian phylogenomics

Resolving the phylogenetic relationships among birds is a classical problem in systematics, and this is particularly so when it comes to understanding the relationships among Neoaves. Previous phylogenetic inference of birds has been limited to mitochondrial genomes or a few nuclear genes. Here, we apply deep brain transcriptome sequencing of nine bird species (several passerines, hummingbirds, dove, parrot, and emu), using next-generation sequencing technology to understand features of transcriptome evolution in birds and how this affects phylogenetic inference, and combine with data from two bird species using first generation technology. The phylogenomic data matrix comprises 1,995 genes and a total of 0.77 Mb of exonic sequence. First, we find an unexpected heterogeneity in the evolution of base composition among avian lineages. There is a pronounced increase in guanine + cytosine (GC) content in the third codon position in several independent lineages, with the strongest effect seen in passerines. Second, we evaluate the effect of GC content variation on phylogenetic reconstruction. We find important inconsistencies between the topologies obtained with or without taking GC variation into account, each supporting different conclusions of past studies and also influencing hypotheses on the evolution of the trait of vocal learning. Third, we demonstrate a link between GC content evolution and recombination rate and, focusing on the zebra finch lineage, find that recombination seems to drive GC content. Although we cannot reveal the causal relationships, this observation is consistent with the model of GC-biased gene conversion. Finally, we use this unparalleled amount of avian sequence data to study the rate of molecular evolution, calibrated by fossil evidence and augmented with data from alligator transcriptome sequencing. There is a 2- to 3-fold variation in substitution rate among lineages with passerines being the most rapidly evolving and ratites the slowest. This study illustrates the potential of next-generation sequencing for phylogenomic studies but also the pitfalls when using genome-wide data with heterogeneous base composition.

Phylogeography of the dark kangaroo mouse, Microdipodops megacephalus: cryptic lineages and dispersal routes in North America's Great Basin

AIM: The rodent genus Microdipodops (kangaroo mice) includes two sand-obligate endemics of the Great Basin Desert: M. megacephalus and M. pallidus. The dark kangaroo mouse, M. megacephalus, is distributed throughout the Great Basin and our principal aims were to formulate phylogenetic hypotheses for this taxon and make phylogeographical comparisons with its congener. LOCATION: The Great Basin Desert of western North America. METHODS: DNA sequence data from three mitochondrial genes were examined from 186 individuals of M. megacephalus, representing 47 general localities. Phylogenetic inference was used to analyse the sequence data. Directional analysis of phylogeographical patterns was used to examine haplotype sharing patterns and recover routes of gene exchange. Haplotype-area curves were constructed to evaluate the relationship between genetic variation and distributional island size for M. megacephalus and M. pallidus. RESULTS: Microdipodops megacephalus is a rare desert rodent (trapping success was 2.67%). Temporal comparison of trapping data shows that kangaroo mice are becoming less abundant in the study area. The distribution has changed slightly since the 1930s but many northern populations now appear to be small, fragmented, or locally extinct. Four principal phylogroups (the Idaho isolate and the western, central and eastern clades) are evident; mean sequence divergence between phylogroups for cytochrome b is c. 8%. Data from haplotype sharing show two trends: a north-south trend and a web-shaped trend. Analyses of haplotype-area curves reveal significant positive relationships. MAIN CONCLUSIONS: The four phylogroups of M. megacephalus appear to represent morphologically cryptic species; in comparison, a companion study revealed two cryptic lineages in M. pallidus. Estimated divergence times of the principal clades of M. megacephalus (c. 2-4 Ma) indicate that these kangaroo mice were Pleistocene invaders into the Great Basin coincident with the formation of sandy habitats. The north-south and web patterns from directional analyses reveal past routes of gene flow and provide evidence for source-sink population regulation. The web pattern was not seen in the companion study of M. pallidus. Significant haplotype-area curves indicate that the distributional islands are now in approximate genetic equilibrium. The patterns described here are potentially useful to conservation biologists and wildlife managers and may serve as a model for other sand-obligate organisms of the Great Basin.

Potential pitfalls of modelling ribosomal RNA data in phylogenetic tree reconstruction: evidence from case studies in the Metazoa

Background

Failure to account for covariation patterns in helical regions of ribosomal RNA (rRNA) genes has the potential to misdirect the estimation of the phylogenetic signal of the data. Furthermore, the extremes of length variation among taxa, combined with regional substitution rate variation can mislead the alignment of rRNA sequences and thus distort subsequent tree reconstructions. However, recent developments in phylogenetic methodology now allow a comprehensive integration of secondary structures in alignment and tree reconstruction analyses based on rRNA sequences, which has been shown to correct some of these problems. Here, we explore the potentials of RNA substitution models and the interactions of specific model setups with the inherent pattern of covariation in rRNA stems and substitution rate variation among loop regions.

Results

We found an explicit impact of RNA substitution models on tree reconstruction analyses. The application of specific RNA models in tree reconstructions is hampered by interaction between the appropriate modelling of covarying sites in stem regions, and excessive homoplasy in some loop regions. RNA models often failed to recover reasonable trees when single-stranded regions are excessively homoplastic, because these regions contribute a greater proportion of the data when covarying sites are essentially downweighted. In this context, the RNA6A model outperformed all other models, including the more parametrized RNA7 and RNA16 models.

Conclusions

Our results depict a trade-off between increased accuracy in estimation of interdependencies in helical regions with the risk of magnifying positions lacking phylogenetic signal. We can therefore conclude that caution is warranted when applying rRNA covariation models, and suggest that loop regions be independently screened for phylogenetic signal, and eliminated when they are indistinguishable from random noise. In addition to covariation and homoplasy, other factors, like non-stationarity of substitution rates and base compositional heterogeneity, can disrupt the signal of ribosomal RNA data. All these factors dictate sophisticated estimation of evolutionary pattern in rRNA data, just as other molecular data require similarly complicated (but different) corrections.

A multigenic perspective on phylogenetic relationships in the largest family of salamanders, the Plethodontidae

Despite several recent studies, the phylogeny of plethodontid salamanders is not yet fully resolved and the phylogenetic positions of several key genera, especially Aneides, Hemidactylium, Hydromantes and Karsenia, are contentious. Here we present a combined dataset of complete mitochondrial genomes and three nuclear loci for 20 species (16 genera) of plethodontids, representing all major clades in the family. The combined dataset without mitochondrial third codon positions provides a fully resolved, statistically well-supported tree. In this topology two major clades are recovered. A northern clade includes Aneides, Desmognathus, Ensatina, Hydromantes, Karsenia, Phaeognathus and Plethodon, with Plethodon being the sister taxon to the rest of the clade. Hydromantes and Karsenia are sister taxa, and Aneides is recovered as the sister taxon to Ensatina. Desmognathus+Phaeognathus form the sister taxon to Aneides+Ensatina. An eastern/southern clade comprises two subclades. One subclade, the spelerpines (Eurycea, Gyrinophilus, Pseudotriton, Stereochilus, Urspelerpes) is the sister taxon to a subclade comprising Hemidactylium, Batrachoseps and the tropical plethodontids (represented by Bolitoglossa, Nototriton and Thorius). In this topology Hemidactylium is well-supported as the sister taxon to Batrachoseps. Only when mitochondrial third codon positions are included using maximum likelihood analysis is Hemidactylium recovered as the sister taxon to Batrachoseps+tropical genera. Hypothesis testing of alternative topologies supports these conclusions. On the basis of these results we propose a conservative taxonomy for Plethodontidae.

Mitogenomic phylogenetic analyses of the Delphinidae with an emphasis on the Globicephalinae

BACKGROUND

Previous DNA-based phylogenetic studies of the Delphinidae family suggest it has undergone rapid diversification, as characterised by unresolved and poorly supported taxonomic relationships (polytomies) for some of the species within this group. Using an increased amount of sequence data we test between alternative hypotheses of soft polytomies caused by rapid speciation, slow evolutionary rate and/or insufficient sequence data, and hard polytomies caused by simultaneous speciation within this family. Combining the mitogenome sequences of five new and 12 previously published species within the Delphinidae, we used Bayesian and maximum-likelihood methods to estimate the phylogeny from partitioned and unpartitioned mitogenome sequences. Further ad hoc tests were then conducted to estimate the support for alternative topologies.

RESULTS

We found high support for all the relationships within our reconstructed phylogenies, and topologies were consistent between the Bayesian and maximum-likelihood trees inferred from partitioned and unpartitioned data. Resolved relationships included the placement of the killer whale (Orcinus orca) as sister taxon to the rest of the Globicephalinae subfamily, placement of the Risso's dolphin (Grampus griseus) within the Globicephalinae subfamily, removal of the white-beaked dolphin (Lagenorhynchus albirostris) from the Delphininae subfamily and the placement of the rough-toothed dolphin (Steno bredanensis) as sister taxon to the rest of the Delphininae subfamily rather than within the Globicephalinae subfamily. The additional testing of alternative topologies allowed us to reject all other putative relationships, with the exception that we were unable to reject the hypothesis that the relationship between L. albirostris and the Globicephalinae and Delphininae subfamilies was polytomic.

CONCLUSION

Despite their rapid diversification, the increased sequence data yielded by mitogenomes enables the resolution of a strongly supported, bifurcating phylogeny, and a chronology of the divergences within the Delphinidae family. This highlights the benefits and potential application of large mitogenome datasets to resolve long-standing phylogenetic uncertainties.