What can the mitochondrial genome reveal about higher‐level phylogeny of the molluscan class Cephalopoda?

We present the first analysis of cephalopod mitochondrial gene order and construct phylogenies based on gene order using Bayesian, distance, and parsimony analysis methods. Analyses included all species of cephalopod for which the whole mitochondrial genome has been sequenced. Where resolution was obtained, these analyses supported division of Neocoleoidea, in which all recent coleoid Cephalopoda can be placed, into Octopodiformes and Decapodiformes. For the same taxa, we also constructed a phylogeny in a maximum likelihood framework based on amino‐acid coded sequence data of all mitochondrial protein coding genes. As well as supporting Octopodiformes and Decapodiformes, amino‐acid analyses established support for Teuthoidea (Oegospida and Myopsida) to the exclusion of Sepiidae, and supported a monophyletic Oegopsida. Partial mitochondrial sequences of additional higher‐level taxa for which whole genome data were not available were subsequently included in the amino‐acid analysis to provide additional information on phylogeny. Spirulida was found to be basal amongst Decapodiformes. Mapping of morphological characters onto our phylogeny and consideration of palaeontological evidence suggests that our phylogeny reflects true evolutionary relationships. © 2011 The Linnean Society of London, Zoological Journal of the Linnean Society, 2011, 161 , 573–586.     

Phylogenomic networks

Phylogenomics is aimed at studying functional and evolutionary aspects of genome biology using phylogenetic analysis of whole genomes. Current approaches to genome phylogenies are commonly founded in terms of phylogenetic trees. However, several evolutionary processes are non tree-like in nature, including recombination and lateral gene transfer (LGT). Phylogenomic networks are a special type of phylogenetic network reconstructed from fully sequenced genomes. The network model, comprising genomes connected by pairwise evolutionary relations, enables the reconstruction of both vertical and LGT events. Modeling genome evolution in the form of a network enables the use of an extensive toolbox developed for network research. The structural properties of phylogenomic networks open up fundamentally new insights into genome evolution.     

Development of taxon‐specific markers for high‐throughput screening of microbial‐feeding nematodes

In an effort to assess the taxonomic identity of large‐scale samplings of nematodes from the Konza Tallgrass Prairie, we sequenced a portion of the 18S rRNA gene and its associated internally transcribed spacer (ITS1) from 74 nematodes encompassing four taxonomic families. From these sequences, we have developed a series of molecular probes to distinguish 16 distinct microbivorous nematode taxa. Using a combination of low power microscopy and taxon‐specific real‐time probes, the 74 nematodes were correctly assigned to their respective taxonomic groups. This optimized method provides a high‐throughput assay to determine nematode identities across larger data sets.     

Small‐ and large‐scale heterogeneity in genetic variation across the collard flycatcher genome: implications for estimating genetic diversity in nonmodel organisms

Population genetic studies in nonmodel organisms are often hampered by a lack of reference genomes that are essential for whole‐genome resequencing. In the light of this, genotyping methods have been developed to effectively eliminate the need for a reference genome, such as genotyping by sequencing or restriction site‐associated DNA sequencing (RAD‐seq). However, what remains relatively poorly studied is how accurately these methods capture both average and variation in genetic diversity across an organism's genome. In this issue of Molecular Ecology Resources, Dutoit et al. (2016) use whole‐genome resequencing data from the collard flycatcher to assess what factors drive heterogeneity in nucleotide diversity across the genome. Using these data, they then simulate how well different sequencing designs, including RAD sequencing, could capture most of the variation in genetic diversity. They conclude that for evolutionary and conservation‐related studies focused on the estimating genomic diversity, researchers should emphasize the number of loci analysed over the number of individuals sequenced.     

Appearance of new tetraspanin genes during vertebrate evolution

A detailed phylogenetic analysis of tetraspanins from 10 fully sequenced metazoan genomes and several fungal and protist genomes gives insight into their evolutionary origins and organization. Our analysis suggests that the superfamily can be divided into four large families. These four families-the CD family, CD63 family, uroplakin family, and RDS family-are further classified as consisting of several ortholog groups. The clustering of several ortholog groups together, such as the CD9/Tsp2/CD81 cluster, suggests functional relatedness of those ortholog groups. The fact that our studies are based on whole genome analysis enabled us to estimate not only the phylogenetic relationships among the tetraspanins, but also the first appearance in the tree of life of certain tetraspanin ortholog groups. Taken together, our data suggest that the tetraspanins are derived from a single (or a few) ancestral gene(s) through sequence divergence, rather than convergence, and that the majority of tetraspanins found in the human genome are vertebrate (21 instances), tetrapod (4 instances), or mammalian (6 instances) inventions.     

Arthropod Phylogenetics in Light of Three Novel Millipede (Myriapoda: Diplopoda) Mitochondrial Genomes with Comments on the Appropriateness of Mitochondrial Genome Sequence Data for Inferring Deep Level Relationships

Background

Arthropods are the most diverse group of eukaryotic organisms, but their phylogenetic relationships are poorly understood. Herein, we describe three mitochondrial genomes representing orders of millipedes for which complete genomes had not been characterized. Newly sequenced genomes are combined with existing data to characterize the protein coding regions of myriapods and to attempt to reconstruct the evolutionary relationships within the Myriapoda and Arthropoda.

Results

The newly sequenced genomes are similar to previously characterized millipede sequences in terms of synteny and length. Unique translocations occurred within the newly sequenced taxa, including one half of the Appalachioria falcifera genome, which is inverted with respect to other millipede genomes. Across myriapods, amino acid conservation levels are highly dependent on the gene region. Additionally, individual loci varied in the level of amino acid conservation. Overall, most gene regions showed low levels of conservation at many sites. Attempts to reconstruct the evolutionary relationships suffered from questionable relationships and low support values. Analyses of phylogenetic informativeness show the lack of signal deep in the trees (i.e., genes evolve too quickly). As a result, the myriapod tree resembles previously published results but lacks convincing support, and, within the arthropod tree, well established groups were recovered as polyphyletic.

Conclusions

The novel genome sequences described herein provide useful genomic information concerning millipede groups that had not been investigated. Taken together with existing sequences, the variety of compositions and evolution of myriapod mitochondrial genomes are shown to be more complex than previously thought. Unfortunately, the use of mitochondrial protein-coding regions in deep arthropod phylogenetics appears problematic, a result consistent with previously published studies. Lack of phylogenetic signal renders the resulting tree topologies as suspect. As such, these data are likely inappropriate for investigating such ancient relationships.     

Synechococcus: 3 billion years of global dominance

Cyanobacteria are among the most important primary producers on the Earth. However, the evolutionary forces driving cyanobacterial species diversity remain largely enigmatic due to both their distinction from macro‐organisms and an undersampling of sequenced genomes. Thus, we present a new genome of a Synechococcus‐like cyanobacterium from a novel evolutionary lineage. Further, we analyse all existing 16S rRNA sequences and genomes of Synechococcus‐like cyanobacteria. Chronograms showed extremely polyphyletic relationships in Synechococcus, which has not been observed in any other cyanobacteria. Moreover, most Synechococcus lineages bifurcated after the Great Oxidation Event, including the most abundant marine picoplankton lineage. Quantification of horizontal gene transfer among 70 cyanobacterial genomes revealed significant differences among studied genomes. Horizontal gene transfer levels were not correlated with ecology, genome size or phenotype, but were correlated with the age of divergence. All findings were synthetized into a novel model of cyanobacterial evolution, characterized by serial convergence of the features, that is multicellularity and ecology.     

DETECTING EVOLUTIONARY TRANSFER OF GENES USING PhIGs1

Organisms have acquired plastids by convoluted paths that have provided multiple opportunities for gene transfer into a host nucleus from intracellular organisms, including the cyanobacterial ancestor of plastids, the proteobacterial ancestor of mitochondria, and both green and red algae whose engulfment has led to secondary acquisition of plastids. These gene movements are most accurately demonstrated by building phylogenetic trees that identify the evolutionary origin of each gene, and one effective tool for this is “PhIGs” (Phylogenetically Inferred Groups; http://PhIGs.org ), a set of databases and computer tools with a Web interface for whole‐genome evolutionary analysis. PhIGs takes as input gene sets of completely sequenced genomes, builds clusters of genes using a novel, graph‐based approach, and reconstructs the evolutionary relationships among all gene families. The user can view and download the sequence alignments, compare intron‐exon structures, and follow links to functional genomic databases. Currently, PhIGs contains 652,756 genes from 45 genomes grouped into 61,059 gene families. Graphical displays show the relative positions of these genes among genomes. PhIGs has been used to detect the evolutionary transfer of hundreds of genes from cyanobacteria and red algae into oömycete nuclear genomes, revealing that even though they have no plastids, their ancestors did, having secondarily acquired them from an intracellular red alga. A great number of genomes are soon to become available that are relevant to our broader understanding of the movement of genes among intracellular compartments after engulfing other organisms, and PhIGs will be an effective tool to interpret these gene movements.     

Phylogenomic analysis of UDP‐dependent glycosyltransferases provides insights into the evolutionary landscape of glycosylation in plant metabolism

Glycosylated metabolites generated by UDP‐dependent glycosyltransferases (UGTs) play critical roles in plant interactions with the environment as well as human and animal nutrition. The evolution of plant UGTs has previously been explored, but with a limited taxon sampling. In this study, 65 fully sequenced plant genomes were analyzed, and stringent criteria for selection of candidate UGTs were applied to ensure a more comprehensive taxon sampling and reliable sequence inclusion. In addition to revealing the overall evolutionary landscape of plant UGTs, the phylogenomic analysis also resolved the phylogenetic association of UGTs from free‐sporing plants and gymnosperms, and identified an additional UGT group (group R) in seed plants. Furthermore, lineage‐specific expansions and contractions of UGT groups were detected in angiosperms, with the total number of UGTs per genome remaining constant generally. The loss of group Q UGTs in Poales and Brassicales, rather than functional convergence in the group Q containing species, was supported by a gene tree of group Q UGTs sampled from many species, and further corroborated by the absence of group Q homologs on the syntenic chromosomal regions in Arabidopsis thaliana (Brassicales). Branch‐site analyses of the group Q UGT gene tree allowed for identification of branches and amino acid sites that experienced episodic positive selection. The positively selected sites are located on the surface of a representative group Q UGT (PgUGT95B2), away from the active site, suggesting their role in protein folding/stability or protein–protein interactions.     

Transposable elements and viruses as factors in adaptation and evolution: an expansion and strengthening of the TE‐Thrust hypothesis

In addition to the strong divergent evolution and significant and episodic evolutionary transitions and speciation we previously attributed to TE‐Thrust, we have expanded the hypothesis to more fully account for the contribution of viruses to TE‐Thrust and evolution. The concept of symbiosis and holobiontic genomes is acknowledged, with particular emphasis placed on the creativity potential of the union of retroviral genomes with vertebrate genomes. Further expansions of the TE‐Thrust hypothesis are proposed regarding a fuller account of horizontal transfer of TEs, the life cycle of TEs, and also, in the case of a mammalian innovation, the contributions of retroviruses to the functions of the placenta. The possibility of drift by TE families within isolated demes or disjunct populations, is acknowledged, and in addition, we suggest the possibility of horizontal transposon transfer into such subpopulations. “Adaptive potential” and “evolutionary potential” are proposed as the extremes of a continuum of “intra‐genomic potential” due to TE‐Thrust. Specific data is given, indicating “adaptive potential” being realized with regard to insecticide resistance, and other insect adaptations. In this regard, there is agreement between TE‐Thrust and the concept of adaptation by a change in allele frequencies. Evidence on the realization of “evolutionary potential” is also presented, which is compatible with the known differential survivals, and radiations of lineages. Collectively, these data further suggest the possibility, or likelihood, of punctuated episodes of speciation events and evolutionary transitions, coinciding with, and heavily underpinned by, intermittent bursts of TE activity.     

Genomic signatures and co‐occurrence patterns of the ultra‐small Saccharimonadia (phylum CPR/Patescibacteria) suggest a symbiotic lifestyle

The size of bacterial genomes is often associated with organismal metabolic capabilities determining ecological breadth and lifestyle. The recently proposed Candidate Phyla Radiation (CPR)/Patescibacteria encompasses mostly unculturable bacterial taxa with relatively small genome sizes with potential for co‐metabolism interdependencies. As yet, little is known about the ecology and evolution of CPR, particularly with respect to how they might interact with other taxa. Here, we reconstructed two novel genomes (namely, Candidatus Saccharibacter sossegus and Candidatus Chaer renensis) of taxa belonging to the class Saccharimonadia within the CPR/Patescibacteria using metagenomes obtained from acid mine drainage (AMD). By testing the hypothesis of genome streamlining or symbiotic lifestyle, our results revealed clear signatures of gene losses in these genomes, such as those associated with de novo biosynthesis of essential amino acids, nucleotides, fatty acids and cofactors. In addition, co‐occurrence analysis provided evidence supporting potential symbioses of these organisms with Hydrotalea sp. in the AMD system. Together, our findings provide a better understanding of the ecology and evolution of CPR/Patescibacteria and highlight the importance of genome reconstruction for studying metabolic interdependencies between unculturable Saccharimonadia representatives.     

Organization of plastid genomes in the freshwater red algal order Batrachospermales (Rhodophyta)

Little is known about genome organization in members of the order Batrachospermales, and the infra‐ordinal relationship remains unresolved. Plastid (cp) genomes of seven members of the freshwater red algal order Batrachospermales were sequenced, with the following aims: (i) to describe the characteristics of cp genomes and compare these with other red algal groups; (ii) to infer the phylogenetic relationships among these members to better understand the infra‐ordinal classification. Cp genomes of Batrachospermales are large, with several cases of gene loss, they are gene‐dense (high gene content for the genome size and short intergenic regions) and have highly conserved gene order. Phylogenetic analyses based on concatenated nucleotide genome data roughly supports the current taxonomic system for the order. Comparative analyses confirm data for members of the class Florideophyceae that cp genomes in Batrachospermales is highly conserved, with little variation in gene composition. However, relevant new features were revealed in our study: genome sizes in members of Batrachospermales are close to the lowest values reported for Florideophyceae; differences in cp genome size within the order are large in comparison with other orders (Ceramiales, Gelidiales, Gracilariales, Hildenbrandiales, and Nemaliales); and members of Batrachospermales have the lowest number of protein‐coding genes among the Florideophyceae. In terms of gene loss, apcF, which encodes the allophycocyanin beta subunit, is absent in all sequenced taxa of Batrachospermales. We reinforce that the interordinal relationships between the freshwater orders Batrachospermales and Thoreales within the Nemaliophycidae is not well resolved due to limited taxon sampling.     

Biological consequences of ancient gene acquisition and duplication in the large genome of Candidatus Solibacter usitatus Ellin6076

Members of the bacterial phylum Acidobacteria are widespread in soils and sediments worldwide, and are abundant in many soils. Acidobacteria are challenging to culture in vitro, and many basic features of their biology and functional roles in the soil have not been determined. Candidatus Solibacter usitatus strain Ellin6076 has a 9.9 Mb genome that is approximately 2-5 times as large as the other sequenced Acidobacteria genomes. Bacterial genome sizes typically range from 0.5 to 10 Mb and are influenced by gene duplication, horizontal gene transfer, gene loss and other evolutionary processes. Our comparative genome analyses indicate that the Ellin6076 large genome has arisen by horizontal gene transfer via ancient bacteriophage and/or plasmid-mediated transduction, and widespread small-scale gene duplications, resulting in an increased number of paralogs. Low amino acid sequence identities among functional group members, and lack of conserved gene order and orientation in regions containing similar groups of paralogs, suggest that most of the paralogs are not the result of recent duplication events. The genome sizes of additional cultured Acidobacteria strains were estimated using pulsed-field gel electrophoresis to determine the prevalence of the large genome trait within the phylum. Members of subdivision 3 had larger genomes than those of subdivision 1, but none were as large as the Ellin6076 genome. The large genome of Ellin6076 may not be typical of the phylum, and encodes traits that could provide a selective metabolic, defensive and regulatory advantage in the soil environment.     

A targeted next‐generation sequencing toolkit for exon‐based cichlid phylogenomics

Cichlid fishes (family Cichlidae) are models for evolutionary and ecological research. Massively parallel sequencing approaches have been successfully applied to study relatively recent diversification in groups of African and Neotropical cichlids, but such technologies have yet to be used for addressing larger‐scale phylogenetic questions of cichlid evolution. Here, we describe a process for identifying putative single‐copy exons from five African cichlid genomes and sequence the targeted exons for a range of divergent (>tens of millions of years) taxa with probes designed from a single reference species (Oreochromis niloticus, Nile tilapia). Targeted sequencing of 923 exons across 10 cichlid species that represent the family's major lineages and geographic distribution resulted in a complete taxon matrix of 564 exons (649 549 bp), representing 559 genes. Maximum likelihood and Bayesian analyses in both species tree and concatenation frameworks yielded the same fully resolved and highly supported topology, which matched the expected backbone phylogeny of the major cichlid lineages. This work adds to the body of evidence that it is possible to use a relatively divergent reference genome for exon target design and successful capture across a broad phylogenetic range of species. Furthermore, our results show that the use of a third‐party laboratory coupled with accessible bioinformatics tools makes such phylogenomics projects feasible for research groups that lack direct access to genomic facilities. We expect that these resources will be used in further cichlid evolution studies and hope the protocols and identified targets will also be useful for phylogenetic studies of a wider range of organisms.     

From Gene Trees to Organismal Phylogeny in Prokaryotes:The Case of the γ-Proteobacteria

The rapid increase in published genomic sequences for bacteria presents the first opportunity to reconstruct evolutionary events on the scale of entire genomes. However, extensive lateral gene transfer (LGT) may thwart this goal by preventing the establishment of organismal relationships based on individual gene phylogenies. The group for which cases of LGT are most frequently documented and for which the greatest density of complete genome sequences is available is the γ-Proteobacteria, an ecologically diverse and ancient group including free-living species as well as pathogens and intracellular symbionts of plants and animals. We propose an approach to multigene phylogeny using complete genomes and apply it to the case of the γ-Proteobacteria. We first applied stringent criteria to identify a set of likely gene orthologs and then tested the compatibilities of the resulting protein alignments with several phylogenetic hypotheses. Our results demonstrate phylogenetic concordance among virtually all (203 of 205) of the selected gene families, with each of the exceptions consistent with a single LGT event. The concatenated sequences of the concordant families yield a fully resolved phylogeny. This topology also received strong support in analyses aimed at excluding effects of heterogeneity in nucleotide base composition across lineages. Our analysis indicates that single-copy orthologous genes are resistant to horizontal transfer, even in ancient bacterial groups subject to high rates of LGT. This gene set can be identified and used to yield robust hypotheses for organismal phylogenies, thus establishing a foundation for reconstructing the evolutionary transitions, such as gene transfer, that underlie diversity in genome content and organization.     

From Gene Trees to Organismal Phylogeny in Prokaryotes:The Case of the γ-Proteobacteria

The rapid increase in published genomic sequences for bacteria presents the first opportunity to reconstruct evolutionary events on the scale of entire genomes. However, extensive lateral gene transfer (LGT) may thwart this goal by preventing the establishment of organismal relationships based on individual gene phylogenies. The group for which cases of LGT are most frequently documented and for which the greatest density of complete genome sequences is available is the γ-Proteobacteria, an ecologically diverse and ancient group including free-living species as well as pathogens and intracellular symbionts of plants and animals. We propose an approach to multigene phylogeny using complete genomes and apply it to the case of the γ-Proteobacteria. We first applied stringent criteria to identify a set of likely gene orthologs and then tested the compatibilities of the resulting protein alignments with several phylogenetic hypotheses. Our results demonstrate phylogenetic concordance among virtually all (203 of 205) of the selected gene families, with each of the exceptions consistent with a single LGT event. The concatenated sequences of the concordant families yield a fully resolved phylogeny. This topology also received strong support in analyses aimed at excluding effects of heterogeneity in nucleotide base composition across lineages. Our analysis indicates that single-copy orthologous genes are resistant to horizontal transfer, even in ancient bacterial groups subject to high rates of LGT. This gene set can be identified and used to yield robust hypotheses for organismal phylogenies, thus establishing a foundation for reconstructing the evolutionary transitions, such as gene transfer, that underlie diversity in genome content and organization.     

Poxvirus genome evolution by gene gain and loss

The poxviruses (Poxviridae) are a family of viruses with double-stranded DNA genomes and substantial numbers (often >200) of genes per genome. We studied the patterns of gene gain and loss over the evolutionary history of 17 poxvirus complete genomes. A phylogeny based on gene family presence/absence showed good agreement with families based on concatenated amino acid sequences of conserved single-copy genes. Gene duplications in poxviruses were often lineage specific, and the most extensively duplicated viral gene families were found in only a few of the genomes analyzed. A total of 34 gene families were found to include a member in at least one of the poxvirus genomes analyzed and at least one animal genome; in 16 (47%) of these families, there was evidence of recent horizontal gene transfer (HGT) from host to virus. Gene families with evidence of HGT included several involved in host immune defense mechanisms (the MHC class I, interleukin-10, interleukin-24, interleukin-18, the interferon gamma receptor, and tumor necrosis factor receptor II) and others (glutaredoxin and glutathione peroxidase) involved in resistance of cells to oxidative stress. Thus "capture" of host genes by HGT has been a recurrent feature of poxvirus evolution and has played an important role in adapting the virus to survive host antiviral defense mechanisms.     

Defense mechanisms against herbivory in Picea: sequence evolution and expression regulation of gene family members in the phenylpropanoid pathway

Background

The insect order Neuroptera encompasses more than 5,700 described species. To date, only three neuropteran mitochondrial genomes have been fully and one partly sequenced. Current knowledge on neuropteran mitochondrial genomes is limited, and new data are strongly required. In the present work, the mitochondrial genome of the ascalaphid owlfly Libelloides macaronius is described and compared with the known neuropterid mitochondrial genomes: Megaloptera, Neuroptera and Raphidioptera. These analyses are further extended to other endopterygotan orders.

Results

The mitochondrial genome of L. macaronius is a circular molecule 15,890 bp long. It includes the entire set of 37 genes usually present in animal mitochondrial genomes. The gene order of this newly sequenced genome is unique among Neuroptera and differs from the ancestral type of insects in the translocation of trnC. The L. macaronius genome shows the lowest A+T content (74.50%) among known neuropterid genomes. Protein-coding genes possess the typical mitochondrial start codons, except for cox1, which has an unusual ACG. Comparisons among endopterygotan mitochondrial genomes showed that A+T content and AT/GC-skews exhibit a broad range of variation among 84 analyzed taxa. Comparative analyses showed that neuropterid mitochondrial protein-coding genes experienced complex evolutionary histories, involving features ranging from codon usage to rate of substitution, that make them potential markers for population genetics/phylogenetics studies at different taxonomic ranks. The 22 tRNAs show variable substitution patterns in Neuropterida, with higher sequence conservation in genes located on the α strand. Inferred secondary structures for neuropterid rrnS and rrnL genes largely agree with those known for other insects. For the first time, a model is provided for domain I of an insect rrnL. The control region in Neuropterida, as in other insects, is fast-evolving genomic region, characterized by AT-rich motifs.

Conclusions

The new genome shares many features with known neuropteran genomes but differs in its low A+T content. Comparative analysis of neuropterid mitochondrial genes showed that they experienced distinct evolutionary patterns. Both tRNA families and ribosomal RNAs show composite substitution pathways. The neuropterid mitochondrial genome is characterized by a complex evolutionary history.