Mitochondrial genomes of five Hyphessobrycon tetras and their phylogenetic implications

To date, the taxonomic status and phylogenetic affinities within Hyphessobrycon, even among other genera in Characidae, remain unclear. Here, we determined five new mitochondrial genomes (mitogenomes) of Hyphessobrycon species (H. elachys, H. flammeus, H. pulchripinnis, H. roseus, and H. sweglesi). The mitogenomes were all classical circular structures, with lengths ranging from 16,008 to 17,224 bp. The type of constitutive genes and direction of the coding strand that appeared in the mitogenomes were identical to those of other species in Characidae. The highest value of the Ka/Ks ratio within 13 protein‐coding genes (PCGs) was found in ND2 with 0.83, suggesting that they were subject to purifying selection in the Hyphessobrycon genus. Comparison of the control region sequences among seven Hyphessobrycon fish revealed that repeat units differ in length and copy number across different species, which led to sharp differences in mitogenome sizes. Phylogenetic trees based on the 13 PCGs did not support taxonomic relationships, as the Hyphessobrycon fish mixed with those from other genera. These data were combined to explore higher level relationships within Characidae and could aid in the understanding of the evolution of this group.     

Complete mitochondrial genomes of five skippers (Lepidoptera: Hesperiidae) and phylogenetic reconstruction of Lepidoptera

We sequenced mitogenomes of five skippers (family Hesperiidae, Lepidoptera) to obtain further insight into the characteristics of butterfly mitogenomes and performed phylogenetic reconstruction using all available gene sequences (PCGs, rRNAs, and tRNAs) from 85 species (20 families in eight superfamilies). The general genomic features found in the butterflies also were found in the five skippers: a high A + T composition (79.3%–80.9%), dominant usage of TAA stop codon, similar skewness pattern in both strands, consistently length intergenic spacer sequence between tRNA^Gln and ND2 (64–87 bp), conserved ATACTAA motif between tRNA^{Ser (UCN)} and ND1, and characteristic features of the A + T-rich region (the ATAGA motif, varying length of poly-T stretch, and poly-A stretch). The start codon for COI was CGA in four skippers as typical, but Lobocla bifasciatus evidently possessed canonical ATG as start codon. All species had the ancestral arrangement tRNA^Asn/tRNA^{Ser (AGN)}, instead of the rearrangement tRNA^{Ser (AGN)}/tRNA^Asn, found in another skipper species (Erynnis). Phylogenetic analyses using all available genes (PCGs, rRNAS, and tRNAs) yielded the consensus superfamilial relationships ((((((Bombycoidea + Noctuoidea + Geometroidea) + Pyraloidea) + Papilionoidea) + Tortricoidea) + Yponomeutoidea) + Hepialoidea), confirming the validity of Macroheterocera (Bombycoidea, Noctuoidea, and Geometroidea in this study) and its sister relationship to Pyraloidea. Within Rhopalocera (butterflies and skippers) the familial relationships (Papilionidae + (Hesperiidae + (Pieridae + ((Lycaenidae + Riodinidae) + Nymphalidae)))) were strongly supported in all analyses (0.98–1 by BI and 96–100 by ML methods), rendering invalid the superfamily status for Hesperioidea. On the other hand, current mitogenome-based phylogeny did not find consistent superfamilial relationships among Noctuoidea, Geometroidea, and Bombycoidea and the familial relationships within Bombycoidea between analyses, requiring further taxon sampling in future studies.     

Isochore chromosome maps of eukaryotic genomes

Analytical DNA ultracentrifugation revealed that eukaryotic genomes are mosaics of isochores: long DNA segments (>300 kb on average) relatively homogeneous in G+C. Important genome features are dependent on this isochore structure, e.g. genes are found predominantly in the GC-richest isochore classes. However, no reliable method is available to rigorously partition the genome sequence into relatively homogeneous regions of different composition, thereby revealing the isochore structure of chromosomes at the sequence level. Homogeneous regions are currently ascertained by plain statistics on moving windows of arbitrary length, or simply by eye on G+C plots. On the contrary, the entropic segmentation method is able to divide a DNA sequence into relatively homogeneous, statistically significant domains. An early version of this algorithm only produced domains having an average length far below the typical isochore size. Here we show that an improved segmentation method, specifically intended to determine the most statistically significant partition of the sequence at each scale, is able to identify the boundaries between long homogeneous genome regions displaying the typical features of isochores. The algorithm precisely locates classes II and III of the human major histocompatibility complex region, two well-characterized isochores at the sequence level, the boundary between them being the first isochore boundary experimentally characterized at the sequence level. The analysis is then extended to a collection of human large contigs. The relatively homogeneous regions we find show many of the features (G+C range, relative proportion of isochore classes, size distribution, and relationship with gene density) of the isochores identified through DNA centrifugation. Isochore chromosome maps, with many potential applications in genomics, are then drawn for all the completely sequenced eukaryotic genomes available.     

Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes

Background  

Molecular evolutionary studies share the common goal of elucidating historical relationships, and the common challenge of adequately sampling taxa and characters. Particularly at low taxonomic levels, recent divergence, rapid radiations, and conservative genome evolution yield limited sequence variation, and dense taxon sampling is often desirable. Recent advances in massively parallel sequencing make it possible to rapidly obtain large amounts of sequence data, and multiplexing makes extensive sampling of megabase sequences feasible. Is it possible to efficiently apply massively parallel sequencing to increase phylogenetic resolution at low taxonomic levels?     

Intraspecific phylogenetic congruence among multiple symbiont genomes

Eukaryotes often form intimate endosymbioses with prokaryotic organisms. Cases in which these symbionts are transmitted cytoplasmically to host progeny create the potential for co-speciation or congruent evolution among the distinct genomes of these partners. If symbionts do not move horizontally between different eukaryotic hosts, strict phylogenetic congruence of their genomes is predicted and should extend to relationships within a single host species. Conversely, even rare 'host shifts' among closely related lineages should yield conflicting tree topologies at the intraspecific level. Here, we investigate the historical associations among four symbiotic genomes residing within an aphid host: the mitochondrial DNA of Uroleucon ambrosiae aphids, the bacterial chromosome of their Buchnera bacterial endosymbionts, and two plasmids associated with Buchnera. DNA sequence polymorphisms provided a significant phylogenetic signal and no homoplasy for each data set, yielding completely and significantly congruent phylogenies for these four genomes and no evidence of horizontal transmission. This study thus provides the first evidence for strictly vertical transmission and 'co-speciation' of symbiotic organisms at the intraspecific level, and represents the lowest phylogenetic level at which such coevolution has been demonstrated. These results may reflect the obligate nature of this intimate mutualism and indicate opportunities for adaptive coevolution among linked symbiont genomes.     

Increasing phylogenetic resolution at low taxonomic levels using massively parallel sequencing of chloroplast genomes

Background

Mycoheterotrophic plants are considered to associate very specifically with fungi. Mycoheterotrophic orchids are mostly associated with ectomycorrhizal fungi in temperate regions, or with saprobes or parasites in tropical regions. Although most mycoheterotrophic orchids occur in the tropics, few studies have been devoted to them, and the main conclusions about their specificity have hitherto been drawn from their association with ectomycorrhizal fungi in temperate regions.

Results

We investigated three Asiatic Neottieae species from ectomycorrhizal forests in Thailand. We found that all were associated with ectomycorrhizal fungi, such as Thelephoraceae, Russulaceae and Sebacinales. Based on ¹³C enrichment of their biomass, they probably received their organic carbon from these fungi, as do mycoheterotrophic Neottieae from temperate regions. Moreover, ¹³C enrichment suggested that some nearby green orchids received part of their carbon from fungi too. Nevertheless, two of the three orchids presented a unique feature for mycoheterotrophic plants: they were not specifically associated with a narrow clade of fungi. Some orchid individuals were even associated with up to nine different fungi.

Conclusion

Our results demonstrate that some green and mycoheterotrophic orchids in tropical regions can receive carbon from ectomycorrhizal fungi, and thus from trees. Our results reveal the absence of specificity in two mycoheterotrophic orchid-fungus associations in tropical regions, in contrast to most previous studies of mycoheterotrophic plants, which have been mainly focused on temperate orchids.     

A simple method for phylogenomic inference using the information of gene content of genomes

Many studies have been contributed to the inferences of phylogenies. Some studies are based on a single-gene (family), and some are based on entire genome data. In this paper, we propose a total loss genome distance approach based on gene content information to inferring phylogenies. Through various simulations, we demonstrate and evaluate the proposed approach. We compare it with some other approaches built upon gene content or extended gene content. Overall, the proposed approach performs equally well as the other methods do and is more efficient than some of the methods. We apply our approach to 34 microbial complete genomes from COG. The reconstructed tree agrees with the results from other approaches and the tree supports the concept of universal trees.     

Neogastropod phylogenetic relationships based on entire mitochondrial genomes

Background  

The Neogastropoda is a highly diversified group of predatory marine snails (Gastropoda: Caenogastropoda). Traditionally, its monophyly has been widely accepted based on several morphological synapomorphies mostly related with the digestive system. However, recent molecular phylogenetic studies challenged the monophyly of Neogastropoda due to the inclusion of representatives of other caenogastropod lineages (e.g. Littorinimorpha) within the group. Neogastropoda has been classified into up to six superfamilies including Buccinoidea, Muricoidea, Olivoidea, Pseudolivoidea, Conoidea, and Cancellarioidea. Phylogenetic relationships among neogastropod superfamilies remain unresolved.     

Reconstruction of certain phylogenetic networks from the genomes at their leaves

A network N is a rooted acyclic digraph. A base-set X for N is a subset of vertices including the root (or outgroup), all leaves, and all vertices of outdegree 1. A simple model of evolution is considered in which all characters are binary and in which back-mutations occur only at hybrid vertices. It is assumed that the genome is known for each member of the base-set X. If the network is known and is assumed to be “normal,” then it is proved that the genome of every vertex is uniquely determined and can be explicitly reconstructed. Under additional hypotheses involving time-consistency and separation of the hybrid vertices, the network itself can also be reconstructed from the genomes of all members of X. An explicit polynomial-time procedure is described for performing the reconstruction.     

Physical maps of the genomes of three Bacillus cereus strains.

NotI restriction maps of the chromosomes from Bacillus cereus ATCC 10876, ATCC 11778, and the B. cereus type strain ATCC 14579 have been established and compared with the previously established map of B. cereus ATCC 10987. Between 10 and 14 NotI fragments were observed, ranging from 15 to 1,300 kb, in digests of DNA from the various strains. The sizes of the genomes varied between 5.4 and 6.3 Mb. The maps were constructed by hybridization of 42 random probes, prepared from B. cereus ATCC 10987 libraries, to fragments from partial and complete NotI digests, separated by pulsed-field gel electrophoresis. Nine probes were specific for ATCC 10987 only. Probes for five B. subtilis and five B. cereus genes were also used. The NotI restriction fragment patterns of the four strains were strikingly different.     

DNA sequence organization in the genomes of five marine invertebrates

The arrangement of repetitive and non-repetitive sequence was studied in the genomic DNA of the oyster (Crassostrea virginica), the surf clam (Spisula solidissima), the horseshoe crab (Limulus polyphemus), a nemertean worm (Cerebratulus lacteus) and a jelly-fish (Aurelia aurita). Except for the jellyfish these animals belong to the protostomial branch of animal evolution, for which little information regarding DNA sequence organization has previously been available. The reassociation kinetics of short (250-300 nucleotide) and long (2,000-3,000 nucleotide) DNA fragments was studied by the hydroxyapatite method. It was shown that in each case a major fraction of the DNA consists of single copy sequences less than about 3,000 nucleotides in length, interspersed with short repetitive sequences. The lengths of the repetitive sequences were estimated by optical hyperchromicity and S1 nuclease measurements made on renaturation products. All the genomes studied include a prominent fraction of interspersed repetitive sequences about 300 nucleotides in length, as well as longer repetitive sequence regions.     

ProGenExpress: Visualization of quantitative data on prokaryotic genomes

Background  

The integration of genomic information with quantitative experimental data is a key component of systems biology. An increasing number of microbial genomes are being sequenced, leading to an increasing amount of data from post-genomics technologies. The genomes of prokaryotes contain many structures of interest, such as operons, pathogenicity islands and prophage sequences, whose behaviour is of interest during infection and disease. There is a need for simple and novel tools to display and analyse data from these integrated datasets, and we have developed ProGenExpress as a tool for visualising arbitrarily complex numerical data in the context of prokaryotic genomes.     

Completing the uric acid degradation pathway through phylogenetic comparison of whole genomes

Mammals that degrade uric acid are not affected by gout or urate kidney stones. It is not fully understood how they convert uric acid into the much more soluble allantoin. Until recently, it had long been thought that urate oxidase was the only enzyme responsible for this conversion. However, detailed studies of the mechanism and regiochemistry of urate oxidation have called this assumption into question, suggesting the existence of other distinct enzymatic activities. Through phylogenetic genome comparison, we identify here two genes that share with urate oxidase a common history of loss or gain events. We show that the two proteins encoded by mouse genes catalyze two consecutive steps following urate oxidation to 5-hydroxyisourate (HIU): hydrolysis of HIU to give 2-oxo-4-hydroxy-4-carboxy-5-ureidoimidazoline (OHCU) and decarboxylation of OHCU to give S-(+)-allantoin. Urate oxidation produces racemic allantoin on a time scale of hours, whereas the full enzymatic complement produces dextrorotatory allantoin on a time scale of seconds. The use of these enzymes in association with urate oxidase could improve the therapy of hyperuricemia.     

Mitochondrial genomes in the Hymenoptera and their utility as phylogenetic markers

Abstract.  In this study, we assessed the ability of mitochondrial genome sequences to recover a test phylogeny of five hymenopteran taxa from which phylogenetic relationships are well accepted. Our analyses indicated that the test phylogeny was well recovered in all nucleotide Bayesian analyses when all the available holometabolan (i.e. outgroup) taxa were included, but only in Bayesian analyses excluding third codon positions when only the hymenopteran representatives and a single outgroup were included. This result suggests that taxon sampling of the outgroup might be as important as taxon sampling of the ingroup when recovering hymenopteran phylogenetic relationships using whole mitochondrial genomes. Parsimony analyses were more sensitive to both taxon sampling and the analytical model than Bayesian analyses, and analyses using the protein dataset did not recover the test phylogeny. In general, mitochondrial genomes did not resolve the position of the Hymenoptera within the Holometabola with confidence, suggesting that an increased taxon sampling, both within the Holometabola and among outgroups, is necessary.     

Considerations for the inclusion of 2x mammalian genomes in phylogenetic analyses

A response to 2x genomes - depth does matter by MC Milinkovitch, R Helaers, E Depiereux, AC Tzika and T Gabaldón. Genome Biol 2010, 11:R16.     

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.     

On the information content of discrete phylogenetic characters

Journal of Mathematical Biology - Phylogenetic inference aims to reconstruct the evolutionary relationships of different species based on genetic (or other) data. Discrete characters are a...     

The isochore patterns of mammalian genomes and their phylogenetic implications

The compositional distributions of high molecular weight DNA fragments from 20 species belonging to 9 out of the 17 eutherian orders were investigated by analytical CsCl density gradient centrifugation and by preparative fractionation in Cs₂SO₄/BAMD density gradients followed by analysis of the fractions in CsCl. These compositional distributions reflect those of the isochores making up the corresponding genomes. A “general distribution” was found in species belonging to eight mammalian orders. A “myomorph distribution” was found in Myomorpha, but not in the other rodent infraorders Sciuromorpha and Histricomorpha, which share the general distribution. Two other distributions were found in a megachiropteran (but not in microchiropteran, which, again, shares the general distribution) and in pangolin (a species from the only genus of the order Pholidota), respectively. The main difference between the general distribution and all other distributions is that the former contains sizable amounts (6–10%) of GC-rich isochores (detected as DNA fragments equal to, or higher than, 1.710 g/cm³ in modal buoyant density), which are scarce, or absent, in the other distributions. This difference is remarkable because gene concentrations in mammalian genomes are paralleled by GC levels, the highest gene concentrations being present in the GC-richest isochores. The compositional distributions of mammalian genomes reported here shed light on mammalian phylogeny. Indeed, all orders investigated, with the exception of Pholidota, seem to share a common ancestor. The compositional patterns of the megachiropteran and of Myomorpha may be derived from the general pattern or have independent origins.