Concatenated Analysis Sheds Light on Early Metazoan Evolution and Fuels a Modern “Urmetazoon” Hypothesis

For more than a century, the origin of metazoan animals has been debated. One aspect of this debate has been centered on what the hypothetical “urmetazoon” bauplan might have been. The morphologically most simply organized metazoan animal, the placozoan Trichoplax adhaerens, resembles an intriguing model for one of several “urmetazoon” hypotheses: the placula hypothesis. Clear support for a basal position of Placozoa would aid in resolving several key issues of metazoan-specific inventions (including, for example, head–foot axis, symmetry, and coelom) and would determine a root for unraveling their evolution. Unfortunately, the phylogenetic relationships at the base of Metazoa have been controversial because of conflicting phylogenetic scenarios generated while addressing the question. Here, we analyze the sum of morphological evidence, the secondary structure of mitochondrial ribosomal genes, and molecular sequence data from mitochondrial and nuclear genes that amass over 9,400 phylogenetically informative characters from 24 to 73 taxa. Together with mitochondrial DNA genome structure and sequence analyses and Hox-like gene expression patterns, these data (1) provide evidence that Placozoa are basal relative to all other diploblast phyla and (2) spark a modernized “urmetazoon” hypothesis.     

Concatenated Analysis Sheds Light on Early Metazoan Evolution and Fuels a Modern “Urmetazoon” Hypothesis

For more than a century, the origin of metazoan animals has been debated. One aspect of this debate has been centered on what the hypothetical “urmetazoon” bauplan might have been. The morphologically most simply organized metazoan animal, the placozoan Trichoplax adhaerens, resembles an intriguing model for one of several “urmetazoon” hypotheses: the placula hypothesis. Clear support for a basal position of Placozoa would aid in resolving several key issues of metazoan-specific inventions (including, for example, head–foot axis, symmetry, and coelom) and would determine a root for unraveling their evolution. Unfortunately, the phylogenetic relationships at the base of Metazoa have been controversial because of conflicting phylogenetic scenarios generated while addressing the question. Here, we analyze the sum of morphological evidence, the secondary structure of mitochondrial ribosomal genes, and molecular sequence data from mitochondrial and nuclear genes that amass over 9,400 phylogenetically informative characters from 24 to 73 taxa. Together with mitochondrial DNA genome structure and sequence analyses and Hox-like gene expression patterns, these data (1) provide evidence that Placozoa are basal relative to all other diploblast phyla and (2) spark a modernized “urmetazoon” hypothesis.     

Systematic Identification of Gene Families for Use as “Markers” for Phylogenetic and Phylogeny-Driven Ecological Studies of Bacteria and Archaea and Their Major Subgroups

With the astonishing rate that genomic and metagenomic sequence data sets are accumulating, there are many reasons to constrain the data analyses. One approach to such constrained analyses is to focus on select subsets of gene families that are particularly well suited for the tasks at hand. Such gene families have generally been referred to as “marker” genes. We are particularly interested in identifying and using such marker genes for phylogenetic and phylogeny-driven ecological studies of microbes and their communities (e.g., construction of species trees, phylogenetic based assignment of metagenomic sequence reads to taxonomic groups, phylogeny-based assessment of alpha- and beta-diversity of microbial communities from metagenomic data). We therefore refer to these as PhyEco (for phylogenetic and phylogenetic ecology) markers. The dual use of these PhyEco markers means that we needed to develop and apply a set of somewhat novel criteria for identification of the best candidates for such markers. The criteria we focused on included universality across the taxa of interest, ability to be used to produce robust phylogenetic trees that reflect as much as possible the evolution of the species from which the genes come, and low variation in copy number across taxa.We describe here an automated protocol for identifying potential PhyEco markers from a set of complete genome sequences. The protocol combines rapid searching, clustering and phylogenetic tree building algorithms to generate protein families that meet the criteria listed above. We report here the identification of PhyEco markers for different taxonomic levels including 40 for “all bacteria and archaea”, 114 for “all bacteria (greatly expanding on the ∼30 commonly used), and 100 s to 1000 s for some of the individual phyla of bacteria. This new list of PhyEco markers should allow much more detailed automated phylogenetic and phylogenetic ecology analyses of these groups than possible previously.     

Discriminating Micropathogen Lineages and Their Reticulate Evolution through Graph Theory-Based Network Analysis: The Case of Trypanosoma cruzi,the Agent of Chagas Disease

Micropathogens (viruses, bacteria, fungi, parasitic protozoa) share a common trait, which is partial clonality, with wide variance in the respective influence of clonality and sexual recombination on the dynamics and evolution of taxa. The discrimination of distinct lineages and the reconstruction of their phylogenetic history are key information to infer their biomedical properties. However, the phylogenetic picture is often clouded by occasional events of recombination across divergent lineages, limiting the relevance of classical phylogenetic analysis and dichotomic trees. We have applied a network analysis based on graph theory to illustrate the relationships among genotypes of Trypanosoma cruzi, the parasitic protozoan responsible for Chagas disease, to identify major lineages and to unravel their past history of divergence and possible recombination events. At the scale of T. cruzi subspecific diversity, graph theory-based networks applied to 22 isoenzyme loci (262 distinct Multi-Locus-Enzyme-Electrophoresis -MLEE) and 19 microsatellite loci (66 Multi-Locus-Genotypes -MLG) fully confirms the high clustering of genotypes into major lineages or “near-clades”. The release of the dichotomic constraint associated with phylogenetic reconstruction usually applied to Multilocus data allows identifying putative hybrids and their parental lineages. Reticulate topology suggests a slightly different history for some of the main “near-clades”, and a possibly more complex origin for the putative hybrids than hitherto proposed. Finally the sub-network of the near-clade T. cruzi I (28 MLG) shows a clustering subdivision into three differentiated lesser near-clades (“Russian doll pattern”), which confirms the hypothesis recently proposed by other investigators. The present study broadens and clarifies the hypotheses previously obtained from classical markers on the same sets of data, which demonstrates the added value of this approach. This underlines the potential of graph theory-based network analysis for describing the nature and relationships of major pathogens, thereby opening stimulating prospects to unravel the organization, dynamics and history of major micropathogen lineages.     

The “Naked Coral” Hypothesis Revisited – Evidence for and Against Scleractinian Monophyly

The relationship between Scleractinia and Corallimorpharia, Orders within Anthozoa distinguished by the presence of an aragonite skeleton in the former, is controversial. Although classically considered distinct groups, some phylogenetic analyses have placed the Corallimorpharia within a larger Scleractinia/Corallimorpharia clade, leading to the suggestion that the Corallimorpharia are “naked corals” that arose via skeleton loss during the Cretaceous from a Scleractinian ancestor. Scleractinian paraphyly is, however, contradicted by a number of recent phylogenetic studies based on mt nucleotide (nt) sequence data. Whereas the “naked coral” hypothesis was based on analysis of the sequences of proteins encoded by a relatively small number of mt genomes, here a much-expanded dataset was used to reinvestigate hexacorallian phylogeny. The initial observation was that, whereas analyses based on nt data support scleractinian monophyly, those based on amino acid (aa) data support the “naked coral” hypothesis, irrespective of the method and with very strong support. To better understand the bases of these contrasting results, the effects of systematic errors were examined. Compared to other hexacorallians, the mt genomes of “Robust” corals have a higher (A+T) content, codon usage is far more constrained, and the proteins that they encode have a markedly higher phenylalanine content, leading us to suggest that mt DNA repair may be impaired in this lineage. Thus the “naked coral” topology could be caused by high levels of saturation in these mitochondrial sequences, long-branch effects or model violations. The equivocal results of these extensive analyses highlight the fundamental problems of basing coral phylogeny on mitochondrial sequence data.     

Conflicting phylogenetic signals at the base of the metazoan tree

A phylogenetic framework is essential for under-standing the origin and evolution of metazoan development. Despite a number of recent molecular studies and a rich fossil record of sponges and cnidarians, the evolutionary relationships of the early branching metazoan groups to each other and to a putative outgroup, the choanoflagellates, remain uncertain. This situation may be the result of the limited amount of phylogenetic information found in single genes and the small number of relevant taxa surveyed. To alleviate the effect of these analytical factors in the phylogenetic recons-truction of early branching metazoan lineages, we cloned multiple protein-coding genes from two choanoflagellates and diverse sponges, cnidarians, and a ctenophore. Comparisons of sequences for alpha-tubulin, beta-tubulin, elongation factor 2, HSP90, and HSP70 robustly support the hypothesis that choanoflagellates are closely affiliated with animals. However, analyses of single and concatenated amino acid sequences fail to resolve the relationships either between early branching metazoan groups or between Metazoa and choano-flagellates. We demonstrate that variable rates of evolution among lineages, sensitivity of the analyses to taxon selection, and conflicts in the phylogenetic signal contained in different amino acid sequences obscure the phylogenetic associations among the early branching Metazoa. These factors raise concerns about the ability to resolve the phylogenetic history of animals with molecular sequences. A consensus view of animal evolution may require investigations of genome-scale characters.     

Non-LTR R2 Element Evolutionary Patterns: Phylogenetic Incongruences,Rapid Radiation and the Maintenance of Multiple Lineages

Retrotransposons of the R2 superclade specifically insert within the 28S ribosomal gene. They have been isolated from a variety of metazoan genomes and were found vertically inherited even if their phylogeny does not always agree with that of the host species. This was explained with the diversification/extinction of paralogous lineages, being proved the absence of horizontal transfer. We here analyze the widest available collection of R2 sequences, either newly isolated from recently sequenced genomes or drawn from public databases, in a phylogenetic framework. Results are congruent with previous analyses, but new important issues emerge. First, the N-terminal end of the R2-B clade protein, so far unknown, presents a new zinc fingers configuration. Second, the phylogenetic pattern is consistent with an ancient, rapid radiation of R2 lineages: being the estimated time of R2 origin (850–600 Million years ago) placed just before the metazoan Cambrian explosion, the wide element diversity and the incongruence with the host phylogeny could be attributable to the sudden expansion of available niches represented by host’s 28S ribosomal genes. Finally, we detect instances of coexisting multiple R2 lineages showing a non-random phylogenetic pattern, strongly similar to that of the “library” model known for tandem repeats: a collection of R2s were present in the ancestral genome and then differentially activated/repressed in the derived species. Models for activation/repression as well as mechanisms for sequence maintenance are also discussed within this framework.     

The mitochondrial genome of the sea anemone Metridium senile (Cnidaria): introns,a paucity of tRNA genes,and a near-standard genetic code.

The circular, 17,443 nucleotide-pair mitochondrial (mt) DNA molecule of the sea anemone, Metridium senile (class Anthozoa, phylum Cnidaria) is presented. This molecule contains genes for 13 energy pathway proteins and two ribosomal (r) RNAs but, relative to other metazoan mtDNAs, has two unique features: only two transfer RNAs (tRNA(f-Met) and tRNA(Trp)) are encoded, and the cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 5 (ND5) genes each include a group I intron. The COI intron encodes a putative homing endonuclease, and the ND5 intron contains the molecule''s ND1 and ND3 genes. Most of the unusual characteristics of other metazoan mtDNAs are not found in M. senile mtDNA: unorthodox translation initiation codons and partial translation termination codons are absent, the use of TGA to specify tryptophan is the only genetic code modification, and both encoded tRNAs have primary and secondary structures closely resembling those of standard tRNAs. Also, with regard to size and secondary structure potential, the mt-s-rRNA and mt-1-rRNA have the least deviation from Escherichia coli 16S and 23S rRNAs of all known metazoan mt-rRNAs. These observations indicate that most of the genetic variations previously reported in metazoan mtDNAs developed after Cnidaria diverged from the common ancestral line of all other Metazoa.     

Estimating Tempo and Mode of Y Chromosome Turnover: Explaining Y Chromosome Loss With the Fragile Y Hypothesis

Chromosomal sex determination is phylogenetically widespread, having arisen independently in many lineages. Decades of theoretical work provide predictions about sex chromosome differentiation that are well supported by observations in both XY and ZW systems. However, the phylogenetic scope of previous work gives us a limited understanding of the pace of sex chromosome gain and loss and why Y or W chromosomes are more often lost in some lineages than others, creating XO or ZO systems. To gain phylogenetic breadth we therefore assembled a database of 4724 beetle species’ karyotypes and found substantial variation in sex chromosome systems. We used the data to estimate rates of Y chromosome gain and loss across a phylogeny of 1126 taxa estimated from seven genes. Contrary to our initial expectations, we find that highly degenerated Y chromosomes of many members of the suborder Polyphaga are rarely lost, and that cases of Y chromosome loss are strongly associated with chiasmatic segregation during male meiosis. We propose the “fragile Y” hypothesis, that recurrent selection to reduce recombination between the X and Y chromosome leads to the evolution of a small pseudoautosomal region (PAR), which, in taxa that require XY chiasmata for proper segregation during meiosis, increases the probability of aneuploid gamete production, with Y chromosome loss. This hypothesis predicts that taxa that evolve achiasmatic segregation during male meiosis will rarely lose the Y chromosome. We discuss data from mammals, which are consistent with our prediction.     

Revised Timeline and Distribution of the Earliest Diverged Human Maternal Lineages in Southern Africa

The oldest extant human maternal lineages include mitochondrial haplogroups L0d and L0k found in the southern African click-speaking forager peoples broadly classified as Khoesan. Profiling these early mitochondrial lineages allows for better understanding of modern human evolution. In this study, we profile 77 new early-diverged complete mitochondrial genomes and sub-classify another 105 L0d/L0k individuals from southern Africa. We use this data to refine basal phylogenetic divergence, coalescence times and Khoesan prehistory. Our results confirm L0d as the earliest diverged lineage (∼172 kya, 95%CI: 149–199 kya), followed by L0k (∼159 kya, 95%CI: 136–183 kya) and a new lineage we name L0g (∼94 kya, 95%CI: 72–116 kya). We identify two new L0d1 subclades we name L0d1d and L0d1c4/L0d1e, and estimate L0d2 and L0d1 divergence at ∼93 kya (95%CI:76–112 kya). We concur the earliest emerging L0d1’2 sublineage L0d1b (∼49 kya, 95%CI:37–58 kya) is widely distributed across southern Africa. Concomitantly, we find the most recent sublineage L0d2a (∼17 kya, 95%CI:10–27 kya) to be equally common. While we agree that lineages L0d1c and L0k1a are restricted to contemporary inland Khoesan populations, our observed predominance of L0d2a and L0d1a in non-Khoesan populations suggests a once independent coastal Khoesan prehistory. The distribution of early-diverged human maternal lineages within contemporary southern Africans suggests a rich history of human existence prior to any archaeological evidence of migration into the region. For the first time, we provide a genetic-based evidence for significant modern human evolution in southern Africa at the time of the Last Glacial Maximum at between ∼21–17 kya, coinciding with the emergence of major lineages L0d1a, L0d2b, L0d2d and L0d2a.     

Mitonuclear Coevolution,but not Nuclear Compensation,Drives Evolution of OXPHOS Complexes in Bivalves

In Metazoa, four out of five complexes involved in oxidative phosphorylation (OXPHOS) are formed by subunits encoded by both the mitochondrial (mtDNA) and nuclear (nuDNA) genomes, leading to the expectation of mitonuclear coevolution. Previous studies have supported coadaptation of mitochondria-encoded (mtOXPHOS) and nuclear-encoded OXPHOS (nuOXPHOS) subunits, often specifically interpreted with regard to the “nuclear compensation hypothesis,” a specific form of mitonuclear coevolution where nuclear genes compensate for deleterious mitochondrial mutations due to less efficient mitochondrial selection. In this study, we analyzed patterns of sequence evolution of 79 OXPHOS subunits in 31 bivalve species, a taxon showing extraordinary mtDNA variability and including species with “doubly uniparental” mtDNA inheritance. Our data showed strong and clear signals of mitonuclear coevolution. NuOXPHOS subunits had concordant topologies with mtOXPHOS subunits, contrary to previous phylogenies based on nuclear genes lacking mt interactions. Evolutionary rates between mt and nuOXPHOS subunits were also highly correlated compared with non-OXPHO-interacting nuclear genes. Nuclear subunits of chimeric OXPHOS complexes (I, III, IV, and V) also had higher dN/dS ratios than Complex II, which is formed exclusively by nuDNA-encoded subunits. However, we did not find evidence of nuclear compensation: mitochondria-encoded subunits showed similar dN/dS ratios compared with nuclear-encoded subunits, contrary to most previously studied bilaterian animals. Moreover, no site-specific signals of compensatory positive selection were detected in nuOXPHOS genes. Our analyses extend the evidence for mitonuclear coevolution to a new taxonomic group, but we propose a reconsideration of the nuclear compensation hypothesis.     

Evolution of Pentameric Ligand-Gated Ion Channels: Pro-Loop Receptors

Pentameric ligand-gated ion channels (pLGICs) are ubiquitous neurotransmitter receptors in Bilateria, with a small number of known prokaryotic homologues. Here we describe a new inventory and phylogenetic analysis of pLGIC genes across all kingdoms of life. Our main finding is a set of pLGIC genes in unicellular eukaryotes, some of which are metazoan-like Cys-loop receptors, and others devoid of Cys-loop cysteines, like their prokaryotic relatives. A number of such “Cys-less” receptors also appears in invertebrate metazoans. Together, those findings draw a new distribution of pLGICs in eukaryotes. A broader distribution of prokaryotic channels also emerges, including a major new archaeal taxon, Thaumarchaeota. More generally, pLGICs now appear nearly ubiquitous in major taxonomic groups except multicellular plants and fungi. However, pLGICs are sparsely present in unicellular taxa, suggesting a high rate of gene loss and a non-essential character, contrasting with their essential role as synaptic receptors of the bilaterian nervous system. Multiple alignments of these highly divergent sequences reveal a small number of conserved residues clustered at the interface between the extracellular and transmembrane domains. Only the “Cys-loop” proline is absolutely conserved, suggesting the more fitting name “Pro loop” for that motif, and “Pro-loop receptors” for the superfamily. The infered molecular phylogeny shows a Cys-loop and a Cys-less clade in eukaryotes, both containing metazoans and unicellular members. This suggests new hypotheses on the evolutionary history of the superfamily, such as a possible origin of the Cys-loop cysteines in an ancient unicellular eukaryote. Deeper phylogenetic relationships remain uncertain, particularly around the split between bacteria, archaea, and eukaryotes.     

Comprehensive Characterization of Human Genome Variation by High Coverage Whole-Genome Sequencing of Forty Four Caucasians

Whole genome sequencing studies are essential to obtain a comprehensive understanding of the vast pattern of human genomic variations. Here we report the results of a high-coverage whole genome sequencing study for 44 unrelated healthy Caucasian adults, each sequenced to over 50-fold coverage (averaging 65.8×). We identified approximately 11 million single nucleotide polymorphisms (SNPs), 2.8 million short insertions and deletions, and over 500,000 block substitutions. We showed that, although previous studies, including the 1000 Genomes Project Phase 1 study, have catalogued the vast majority of common SNPs, many of the low-frequency and rare variants remain undiscovered. For instance, approximately 1.4 million SNPs and 1.3 million short indels that we found were novel to both the dbSNP and the 1000 Genomes Project Phase 1 data sets, and the majority of which (∼96%) have a minor allele frequency less than 5%. On average, each individual genome carried ∼3.3 million SNPs and ∼492,000 indels/block substitutions, including approximately 179 variants that were predicted to cause loss of function of the gene products. Moreover, each individual genome carried an average of 44 such loss-of-function variants in a homozygous state, which would completely “knock out” the corresponding genes. Across all the 44 genomes, a total of 182 genes were “knocked-out” in at least one individual genome, among which 46 genes were “knocked out” in over 30% of our samples, suggesting that a number of genes are commonly “knocked-out” in general populations. Gene ontology analysis suggested that these commonly “knocked-out” genes are enriched in biological process related to antigen processing and immune response. Our results contribute towards a comprehensive characterization of human genomic variation, especially for less-common and rare variants, and provide an invaluable resource for future genetic studies of human variation and diseases.     

Genetic Stratigraphy of Key Demographic Events in Arabia

At the crossroads between Africa and Eurasia, Arabia is necessarily a melting pot, its peoples enriched by successive gene flow over the generations. Estimating the timing and impact of these multiple migrations are important steps in reconstructing the key demographic events in the human history. However, current methods based on genome-wide information identify admixture events inefficiently, tending to estimate only the more recent ages, as here in the case of admixture events across the Red Sea (∼8–37 generations for African input into Arabia, and 30–90 generations for “back-to-Africa” migrations). An mtDNA-based founder analysis, corroborated by detailed analysis of the whole-mtDNA genome, affords an alternative means by which to identify, date and quantify multiple migration events at greater time depths, across the full range of modern human history, albeit for the maternal line of descent only. In Arabia, this approach enables us to infer several major pulses of dispersal between the Near East and Arabia, most likely via the Gulf corridor. Although some relict lineages survive in Arabia from the time of the out-of-Africa dispersal, 60 ka, the major episodes in the peopling of the Peninsula took place from north to south in the Late Glacial and, to a lesser extent, the immediate post-glacial/Neolithic. Exchanges across the Red Sea were mainly due to the Arab slave trade and maritime dominance (from ∼2.5 ka to very recent times), but had already begun by the early Holocene, fuelled by the establishment of maritime networks since ∼8 ka. The main “back-to-Africa” migrations, again undetected by genome-wide dating analyses, occurred in the Late Glacial period for introductions into eastern Africa, whilst the Neolithic was more significant for migrations towards North Africa.     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny Chromosomes

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.     

Evolution of modern birds revealed by mitogenomics: timing the radiation and origin of major orders

Mitochondrial (mt) genes and genomes are among the major sources of data for evolutionary studies in birds. This places mitogenomic studies in birds at the core of intense debates in avian evolutionary biology. Indeed, complete mt genomes are actively been used to unveil the phylogenetic relationships among major orders, whereas single genes (e.g., cytochrome c oxidase I [COX1]) are considered standard for species identification and defining species boundaries (DNA barcoding). In this investigation, we study the time of origin and evolutionary relationships among Neoaves orders using complete mt genomes. First, we were able to solve polytomies previously observed at the deep nodes of the Neoaves phylogeny by analyzing 80 mt genomes, including 17 new sequences reported in this investigation. As an example, we found evidence indicating that columbiforms and charadriforms are sister groups. Overall, our analyses indicate that by improving the taxonomic sampling, complete mt genomes can solve the evolutionary relationships among major bird groups. Second, we used our phylogenetic hypotheses to estimate the time of origin of major avian orders as a way to test if their diversification took place prior to the Cretaceous/Tertiary (K/T) boundary. Such timetrees were estimated using several molecular dating approaches and conservative calibration points. Whereas we found time estimates slightly younger than those reported by others, most of the major orders originated prior to the K/T boundary. Finally, we used our timetrees to estimate the rate of evolution of each mt gene. We found great variation on the mutation rates among mt genes and within different bird groups. COX1 was the gene with less variation among Neoaves orders and the one with the least amount of rate heterogeneity across lineages. Such findings support the choice of COX 1 among mt genes as target for developing DNA barcoding approaches in birds.     

Old trees, new trees--is there any progress?

The amount of comparative data for phylogenetic analyses is constantly increasing. Data come from different directions such as morphology, molecular genetics, developmental biology and paleontology. With the increasing diversity of data and of analytical tools, the number of competing hypotheses on phylogenetic relationships rises, too. The choice of the phylogenetic tree as a basis for the interpretation of new data is important, because different trees will support different evolutionary interpretations of the data investigated. I argue here that, although many problematic aspects exist, there are several phylogenetic relationships that are supported by the majority of analyses and may be regarded as something like a robust backbone. This accounts, for example, for the monophyly of Metazoa, Bilateria, Deuterostomia, Protostomia (= Gastroneuralia), Gnathifera, Spiralia, Trochozoa and Arthropoda and probably also for the branching order of diploblastic taxa (“Porifera”, Trichoplax adhaerens, Cnidaria and Ctenophora). Along this “backbone”, there are several problematic regions, where either monophyly is questionable and/or where taxa “rotate” in narrow regions of the tree. This is illustrated exemplified by the probable paraphyly of Porifera and the phylogenetic relationships of basal spiralian taxa. Two problems span wider regions of the tree: the position of Arthropoda either as the sister taxon of Annelida (= Articulata) or of Cycloneuralia (= Ecdysozoa) and the position of tentaculate taxa either as sister taxa of Deuterostomia (= Radialia) or within the taxon Spiralia. The backbone makes it possible to develop a basic understanding of the evolution of genes, molecules and structures in metazoan animals.     

Subtle Population Genetic Structure in Yelloweye Rockfish (Sebastes ruberrimus) Is Consistent with a Major Oceanographic Division in British Columbia,Canada

The boundaries between oceanographic domains often function as dispersal barriers for many temperate marine species with a dispersive pelagic larval phase. Yelloweye rockfish (Sebastes ruberrimus, YR) are widely distributed across the northeastern Pacific Ocean, inhabiting coastal rocky reefs from the Aleutian Islands in Alaska through southern California. This species exhibits an extended pelagic larval duration and has the capacity for long distance larval transport. We assayed 2,862 YR individuals from 13 general areas in the northeast Pacific Ocean for allelic variation at nine microsatellite loci. Bayesian model-based clustering analyses grouped individuals from the Strait of Georgia (SG) into a distinct genetic cluster, while individuals from outer coastal water locations (OCLs) were partitioned equally across two genetic clusters, including the cluster associated with the SG fish. Pairwise F_ST values were consistently an order of magnitude higher for comparisons between the SG and OCLs than they were for all OCL-OCL comparisons (∼0.016 vs. ∼0.001). This same pattern was observed across two time points when individuals were binned into an “old” and “young” group according to birth year (old: ∼0.020 vs. 0.0003; young: ∼0.020 vs. ∼0.004). Additionally, mean allelic richness was markedly lower within the SG compared to the OCLs (8.00 vs. 10.54–11.77). These results indicate that the Strait of Georgia “deep-basin” estuary oceanographic domain acts as a dispersal barrier from the outer coastal waters via the Juan de Fuca Strait. Alternatively, selection against maladapted dispersers across this oceanographic transition may underlie the observed genetic differentiation between the Georgia basin and the outer coastal waters, and further work is needed to confirm the SG-OCL divide acts as a barrier to larval dispersal.