Comparative phylogeography of two Northern Rocky Mountain endemics: the widespread Anguispira kochi occidentalis and the narrow-range Anguispira nimapuna (Gastropoda: Discidae).

The Northern Rocky Mountain ecosystem supports rich biological diversity with many endemic and rare species. Extant endemics display two biogeographic patterns: widespread species with fragmented populations, and narrow-range endemics. These distributions are shown by the congeneric snails Anguispira kochi occidentalis and Anguispira nimapuna. These two taxa are disjunct from the remaining species of the genus, which achieves its greatest diversity in eastern North America. Given the disjunct nature of A. k. occidentalis and A. nimapuna, we here present a mtDNA phylogeny of the genus that includes both eastern and western species to assess the phylogenetic position of A. k. occidentalis and A. nimapuna. We then reconstruct the demographic history of A. k. occidentalis and A. nimapuna by analysing current patterns of genetic variation and interpreting the results considering the historical biogeography of the region. Both A. k. occidentalis and A. nimapuna represent unique taxa that are genetically and geographically distinct from their congeners. The current distribution and genetic structure of A. k. occidentalis has been shaped by both historical isolation in refugia and more recent northward shifts, whereas A. nimapuna is represented by two populations with shallow divergence in an area of long-term habitat stability.

Phylogenomic analyses recover a clade of large-bodied decapodiform cephalopods.

Phylogenetic relationships among the squids and cuttlefishes (Cephalopoda:Decapodiformes) have resisted clarification for decades, despite multiple analyses of morphological, molecular and combined data sets. More recently, analyses of complete mitochondrial genomes and hundreds of nuclear loci have yielded similarly ambiguous results. In this study, we re-evaluate hypotheses of decapodiform relationships by increasing taxonomic breadth and utilizing higher-quality genome and transcriptome data for several taxa. We also employ analytical approaches to (1) identify contamination in transcriptome data, (2) better assess model adequacy, and (3) account for potential biases. Using this larger data set, we consistently recover a clade comprising Myopsida (closed-eye squid), Sepiida (cuttlefishes), and Oegopsida (open-eye squid) that is sister to a Sepiolida (bobtail and bottletail squid) clade. Idiosepiida (pygmy squid) is consistently recovered as the sister group to all sampled decapodiform lineages. Further, a weighted Shimodaira-Hasegawa test applied to one of our larger data matrices rejects all alternatives to these ordinal-level relationships. At present, available nuclear genome-scale data support nested clades of relatively large-bodied decapodiform cephalopods to the exclusion of pygmy squids, but improved taxon sampling and additional genomic data will be needed to test these novel hypotheses rigorously.

Spiralian Genomes Reveal Gene Family Expansions Associated with Adaptation to Freshwater.

The colonization of freshwater habitats by marine-adapted organisms represents a major transition that has only occurred a few times in the evolution of animals. Only around half of the extant animal phyla have representatives in both marine and freshwater environments and even within those phyla some major clades are restricted to marine environments. Moving from marine to freshwater environments can create severe osmotic and ionic stresses and the mechanisms that animals have used to adapt to those stresses are still not well understood. In this study, we downloaded amino acid sequence data from 11 spiralian animal genomes (four freshwater taxa representing four different phyla as well as 7 marine taxa) and identified a number of gene family expansions that have occurred exclusively in the freshwater lineages. Further investigation of these gene families and the timing and nature of their expansions will illuminate one of the major evolutionary transitions in the history of life on Earth.

The grass squid Pickfordiateuthis pulchella is a paedomorphic loliginid.

The wide disparity in adult body size observed both within and among animal taxa has long attracted widespread interest, with several general rules having been proposed to explain trends in body size evolution. Adult body size disparity among the cephalopod mollusks is remarkable, with adult body sizes ranging from a few centimeters to several meters. Some of the smallest cephalopods are found within Pickfordiateuthis, a group comprising three described species of squid found in the western Atlantic and tropical eastern Pacific. Pickfordiateuthis pulchella, the type species of the genus, was initially proposed to be closely related to the loliginid squids (Loliginidae), with subsequent descriptions of additional species supporting a placement within Loliginidae. Pickfordiateuthis is remarkable in that all species reach sexual maturity at about one-fifth to one-tenth the size seen in most loliginid species. To date, no phylogenetic analyses have included representatives of Pickfordiateuthis. To infer the phylogenetic position of Pickfordiateuthis and explore its implications for body size evolution, we collected specimens of Pickfordiateuthis pulchella from Brazilian waters and sequenced regions of two loci-the mitochondrial large ribosomal subunit (rrnL a.k.a. 16S) gene and the nuclear gene rhodopsin. Maximum likelihood and Bayesian analyses of these sequences support a placement of Pickfordiateuthis pulchella as sister to a clade comprising the Western Hemisphere loliginid genera Doryteuthis and Lolliguncula. Analyses of body size evolution within Loliginidae suggest that a shift to a smaller body size optimum occurred along the lineage leading to P. pulchella, with some evidence of shifts toward larger sizes in the ancestors of Loligo and Sepioteuthis; these inferences seem to be robust to phylogenetic uncertainty and incomplete taxon sampling. The small size and juvenile-like morphological traits seen in adult Pickfordiateuthis (e.g., sepiolid-like fins and biserial sucker arrangement in the tentacles) may be due to paedomorphosis.

Phylogenomic Analysis of a Putative Missing Link Sparks Reinterpretation of Leech Evolution.

Leeches (Hirudinida) comprise a charismatic, yet often maligned group of organisms. Despite their ecological, economic, and medical importance, a general consensus on the phylogenetic relationships of major hirudinidan lineages is lacking. This absence of a consistent, robust phylogeny of early-diverging lineages has hindered our understanding of the underlying processes that enabled evolutionary diversification of this clade. Here, we used an anchored hybrid enrichment-based phylogenomic approach, capturing hundreds of loci to investigate phylogenetic relationships among major hirudinidan lineages and their closest living relatives. Our results suggest that a dramatic reinterpretation of early leech evolution is warranted. We recovered Branchiobdellida as sister to a clade that includes all major lineages of hirudinidans, but found Acanthobdella to be nested within Oceanobdelliformes. These results cast doubt on the utility of Acanthobdella as a "missing link" used to explain the origin of blood-feeding in hirudineans. Further, our results support a deep divergence between predominantly marine and freshwater lineages, while not supporting the reciprocal monophyly of jawed and proboscis-bearing leeches. To sum up, our phylogenomic resolution of early-diverging leeches provides a necessary foundation for illuminating the evolution of host-symbiont associations and key adaptations that have allowed leeches to colonize a wide diversity of habitats worldwide.

Na+ /K+ -ATPase gene duplications in clitellate annelids are associated with freshwater colonization.

Major habitat transitions, such as those from marine to freshwater habitats or from aquatic to terrestrial habitats, have occurred infrequently in animal evolution and may represent a barrier to diversification. Identifying genomic events associated with these transitions can help us better understand mechanisms that allow animals to cross these barriers and diversify in new habitats. Study of the Capitella telata and Helobdella robusta genomes allows examination of one such habitat transition (marine to freshwater) in Annelida. Initial examination of these genomes indicated that the freshwater leech H. robusta contains many more copies (12) of the sodium-potassium pump alpha-subunit (Na+ /K+ -ATPase) gene than does the marine polychaete C. telata (2). The sodium-potassium pump plays a key role in maintenance of cellular ionic balance and osmoregulation, and Na+ /K+ -ATPase duplications may have helped annelids invade and diversify in freshwater habitats. To assess whether the timing of Na+ /K+ -ATPase duplications coincided with the marine-to-freshwater transition in Clitellata, we used transcriptomic data from 18 annelid taxa, along with the two genomes, to infer a species phylogeny and identified Na+ /K+ -ATPase gene transcripts in order to infer the timing of gene duplication events using tree-based methods. The inferred timing of Na+ /K+ -ATPase duplication events is consistent with the timing of the initial marine-to-freshwater transition early in the history of clitellate annelids, supporting the hypothesis that gene duplications may have played a role in the annelid diversification into freshwater habitats.

Clarifying the phylogenetic relationships and taxonomy of Stenonema, Stenacron and Maccaffertium, three common eastern North American mayfly genera.

Stenonema, Stenacron, and Maccaffertium are three closely related genera of mayflies (Ephemeroptera:Heptageniidae) commonly found across North America. Due to their primarily aquatic life history and sensitivity to aquatic pollutants, these mayflies are often used as water quality indicators. However, there is little morphological variation within these genera, leading to difficulties in identification and rampant taxonomic confusion, limiting their utility as bioindicators. In an attempt to resolve the phylogenetic relationships of Stenonema, Stenacron, and Maccaffertium, and to clarify their higher-level classifications, we sequenced regions of two mitochondrial genes (cytochrome oxidase subunit 1 (cox1) and 16S ribosomal RNA (rrnl)) and two nuclear genes (Wingless (Wg) and histone H3) from 60 individuals representing most of the described species in these genera and included data from representatives of three heptageniid genera (Kageronia, Macdunnoa and Pseudiron) proposed in previous studies to be closely related to our focal taxa as well as two more distantly related heptageniid genera (Epeorus and Heptagenia) to root the phylogenies. Maximum likelihood and Bayesian analysis were conducted on single-gene and concatenated multi-gene data sets and species tree methods were utilized to resolve relationships. These analyses resolved Stenacron as a monophyletic group sister to a clade comprising Macdunnoa, Maccaffertium and Stenonema. Maccaffertium was found to be paraphyletic, with Stenonema femoratum resolved within Maccaffertium as sister to M. mexicanum. Many relationships remained unresolved or varied across analyses, making revision of the classification based on phylogenetic considerations challenging. To minimize confusion while naming clades and acknowledging uncertainty in our phylogenetic conclusions, we redefine Stenonema to include Maccaffertium and propose three subgenera-Stenonema, Maccaffertium and Lewisa- for key well-supported clades.

Assessing the utility of transcriptome data for inferring phylogenetic relationships among coleoid cephalopods.

Historically, deep-level relationships within the molluscan class Cephalopoda (squids, cuttlefishes, octopods and their relatives) have remained elusive due in part to the considerable morphological diversity of extant taxa, a limited fossil record for species that lack a calcareous shell and difficulties in sampling open ocean taxa. Many conflicts identified by morphologists in the early 1900s remain unresolved today in spite of advances in morphological, molecular and analytical methods. In this study we assess the utility of transcriptome data for resolving cephalopod phylogeny, with special focus on the orders of Decapodiformes (open-eye squids, bobtail squids, cuttlefishes and relatives). To do so, we took new and previously published transcriptome data and used a unique cephalopod core ortholog set to generate a dataset that was subjected to an array of filtering and analytical methods to assess the impacts of: taxon sampling, ortholog number, compositional and rate heterogeneity and incongruence across loci. Analyses indicated that datasets that maximized taxonomic coverage but included fewer orthologs were less stable than datasets that sacrificed taxon sampling to increase the number of orthologs. Clades recovered irrespective of dataset, filtering or analytical method included Octopodiformes (Vampyroteuthis infernalis + octopods), Decapodiformes (squids, cuttlefishes and their relatives), and orders Oegopsida (open-eyed squids) and Myopsida (e.g., loliginid squids). Ordinal-level relationships within Decapodiformes were the most susceptible to dataset perturbation, further emphasizing the challenges associated with uncovering relationships at deep nodes in the cephalopod tree of life.

Phylogenomic analyses of Crassiclitellata support major Northern and Southern Hemisphere clades and a Pangaean origin for earthworms.

BACKGROUND: Earthworms (Crassiclitellata) are a diverse group of annelids of substantial ecological and economic importance. Earthworms are primarily terrestrial infaunal animals, and as such are probably rather poor natural dispersers. Therefore, the near global distribution of earthworms reflects an old and likely complex evolutionary history. Despite a long-standing interest in Crassiclitellata, relationships among and within major clades remain unresolved. METHODS: In this study, we evaluate crassiclitellate phylogenetic relationships using 38 new transcriptomes in combination with publicly available transcriptome data. Our data include representatives of nearly all extant earthworm families and a representative of Moniligastridae, another terrestrial annelid group thought to be closely related to Crassiclitellata. We use a series of differentially filtered data matrices and analyses to examine the effects of data partitioning, missing data, compositional and branch-length heterogeneity, and outgroup inclusion. RESULTS AND DISCUSSION: We recover a consistent, strongly supported ingroup topology irrespective of differences in methodology. The topology supports two major earthworm clades, each of which consists of a Northern Hemisphere subclade and a Southern Hemisphere subclade. Divergence time analysis results are concordant with the hypothesis that these north-south splits are the result of the breakup of the supercontinent Pangaea. CONCLUSIONS: These results support several recently proposed revisions to the classical understanding of earthworm phylogeny, reveal two major clades that seem to reflect Pangaean distributions, and raise new questions about earthworm evolutionary relationships.

Patterns in wetland microbial community composition and functional gene repertoire associated with methane emissions.

UNLABELLED: Wetland restoration on peat islands previously drained for agriculture has potential to reverse land subsidence and sequester atmospheric carbon dioxide as peat accretes. However, the emission of methane could potentially offset the greenhouse gas benefits of captured carbon. As microbial communities play a key role in governing wetland greenhouse gas fluxes, we are interested in how microbial community composition and functions are associated with wetland hydrology, biogeochemistry, and methane emission, which is critical to modeling the microbial component in wetland methane fluxes and to managing restoration projects for maximal carbon sequestration. Here, we couple sequence-based methods with biogeochemical and greenhouse gas measurements to interrogate microbial communities from a pilot-scale restored wetland in the Sacramento-San Joaquin Delta of California, revealing considerable spatial heterogeneity even within this relatively small site. A number of microbial populations and functions showed strong correlations with electron acceptor availability and methane production; some also showed a preference for association with plant roots. Marker gene phylogenies revealed a diversity of major methane-producing and -consuming populations and suggested novel diversity within methanotrophs. Methanogenic archaea were observed in all samples, as were nitrate-, sulfate-, and metal-reducing bacteria, indicating that no single terminal electron acceptor was preferred despite differences in energetic favorability and suggesting spatial microheterogeneity and microniches. Notably, methanogens were negatively correlated with nitrate-, sulfate-, and metal-reducing bacteria and were most abundant at sampling sites with high peat accretion and low electron acceptor availability, where methane production was highest. IMPORTANCE: Wetlands are the largest nonanthropogenic source of atmospheric methane but also a key global carbon reservoir. Characterizing belowground microbial communities that mediate carbon cycling in wetlands is critical to accurately predicting their responses to changes in land management and climate. Here, we studied a restored wetland and revealed substantial spatial heterogeneity in biogeochemistry, methane production, and microbial communities, largely associated with the wetland hydraulic design. We observed patterns in microbial community composition and functions correlated with biogeochemistry and methane production, including diverse microorganisms involved in methane production and consumption. We found that methanogenesis gene abundance is inversely correlated with genes from pathways exploiting other electron acceptors, yet the ubiquitous presence of genes from all these pathways suggests that diverse electron acceptors contribute to the energetic balance of the ecosystem. These investigations represent an important step toward effective management of wetlands to reduce methane flux to the atmosphere and enhance belowground carbon storage.

Ancestral state reconstruction for Dendroctonus bark beetles: evolution of a tree killer.

While most bark beetles attack only dead or weakened trees, many species in the genus Dendroctonus have the ability to kill healthy conifers through mass attack of the host tree, and can exhibit devastating outbreaks. Other species in this group are able to successfully colonize trees in small numbers without killing the host. We reconstruct the evolution of these ecological and life history traits, first classifying the extant Dendroctonus species by attack type (mass or few), outbreaks (yes or no), host genus (Pinus and others), location of attacks on the tree (bole, base, etc.), whether the host is killed (yes or no), and if the larvae are gregarious or have individual galleries (yes or no). We then estimated a molecular phylogeny for a data set of cytochrome oxidase I sequences sampled from nearly all Dendroctonus species, and used this phylogeny to reconstruct the ancestral state at various nodes on the tree, employing maximum parsimony, maximum likelihood, and Bayesian methods. Our reconstructions suggest that extant Dendroctonus species likely evolved from an ancestor that killed host pines through mass attack of the bole, had individual larvae, and exhibited outbreaks. The ability to colonize a host tree in small numbers (as well as gregarious larvae and attacks at the tree base) apparently evolved later, possibly as two separate events in different clades. It is likely that tree mortality and outbreaks have been continuing features of the interaction between conifers and Dendroctonus bark beetles.

Mitochondrial DNA analysis of the bark beetle predator Thanasimus dubius F. (Coleoptera: Cleridae) reveals regional genetic differentiation.

The checkered beetle, Thanasimus dubius F., is an important predator of scolytid bark beetles that attack conifers. Relatively few studies exist that have addressed the population genetics of predatory beetles, especially those with potential as biological control agents. This study was conducted to investigate the population genetics of T. dubius across a large part of its range in the eastern United States. A 464-base pair portion of the mitochondrial cytochrome c oxidase subunit I was sequenced for 85 individuals resulting in 60 haplotypes. Analysis of molecular variance was conducted on the resulting haplotypes for all populations and as a hierarchical analysis between populations defined as broad-scale northern and southern groups. Results indicate a significant overall Phi ST = 0.220 (P < 0.001) for all populations with the hierarchical analysis revealing that this significant Phi ST is due to structuring of the populations between northern and southern regions (Phi CT = 0.388, P < 0.009). The observed genetic structure is possibly due to the discontinuous distribution of pine trees, which act as hosts for the prey of T. dubius, which has occurred historically in the central region of the United States that has been covered by prairie.

Discovery of the photosynthetic relatives of the "Maltese mushroom" Cynomorium.

BACKGROUND: Although recent molecular phylogenetic studies have identified the photosynthetic relatives of several enigmatic holoparasitic angiosperms, uncertainty remains for the last parasitic plant order, Balanophorales, often considered to include two families, Balanophoraceae and Cynomoriaceae. The nonphotosynthetic (holoparasitic) flowering plant Cynomorium coccineum has long been known to the Muslim world as "tarthuth" and to Europeans as the "Maltese mushroom"; C. songaricum is known in Chinese medicine as "suo yang." Interest in these plants is increasing and they are being extensively collected from wild populations for use in herbal medicines. RESULTS: Here we report molecular phylogenetic analyses of nuclear ribosomal DNA and mitochondrial matR sequence data that strongly support the independent origin of Balanophoraceae and Cynomoriaceae. Analyses of single gene and combined gene data sets place Cynomorium in Saxifragales, possibly near Crassulaceae (stonecrop family). Balanophoraceae appear related to Santalales (sandalwood order), a position previously suggested from morphological characters that are often assumed to be convergent. CONCLUSION: Our work shows that Cynomorium and Balanophoraceae are not closely related as indicated in all past and present classifications. Thus, morphological features, such as inflorescences bearing numerous highly reduced flowers, are convergent and were attained independently by these two holoparasite lineages. Given the widespread harvest of wild Cynomorium species for herbal medicines, we here raise conservation concerns and suggest that further molecular phylogenetic work is needed to identify its photosynthetic relatives. These relatives, which will be easier to cultivate, should then be examined for phytochemical activity purported to be present in the more sensitive Cynomorium.

Phylogenetic inference in Rafflesiales: the influence of rate heterogeneity and horizontal gene transfer.

BACKGROUND: The phylogenetic relationships among the holoparasites of Rafflesiales have remained enigmatic for over a century. Recent molecular phylogenetic studies using the mitochondrial matR gene placed Rafflesia, Rhizanthes and Sapria (Rafflesiaceae s. str.) in the angiosperm order Malpighiales and Mitrastema (Mitrastemonaceae) in Ericales. These phylogenetic studies did not, however, sample two additional groups traditionally classified within Rafflesiales (Apodantheaceae and Cytinaceae). Here we provide molecular phylogenetic evidence using DNA sequence data from mitochondrial and nuclear genes for representatives of all genera in Rafflesiales. RESULTS: Our analyses indicate that the phylogenetic affinities of the large-flowered clade and Mitrastema, ascertained using mitochondrial matR, are congruent with results from nuclear SSU rDNA when these data are analyzed using maximum likelihood and Bayesian methods. The relationship of Cytinaceae to Malvales was recovered in all analyses. Relationships between Apodanthaceae and photosynthetic angiosperms varied depending upon the data partition: Malvales (3-gene), Cucurbitales (matR) or Fabales (atp1). The latter incongruencies suggest that horizontal gene transfer (HGT) may be affecting the mitochondrial gene topologies. The lack of association between Mitrastema and Ericales using atp1 is suggestive of HGT, but greater sampling within eudicots is needed to test this hypothesis further. CONCLUSIONS: Rafflesiales are not monophyletic but composed of three or four independent lineages (families): Rafflesiaceae, Mitrastemonaceae, Apodanthaceae and Cytinaceae. Long-branch attraction appears to be misleading parsimony analyses of nuclear small-subunit rDNA data, but model-based methods (maximum likelihood and Bayesian analyses) recover a topology that is congruent with the mitochondrial matR gene tree, thus providing compelling evidence for organismal relationships. Horizontal gene transfer appears to be influencing only some taxa and some mitochondrial genes, thus indicating that the process is acting at the single gene (not whole genome) level.

Should we be worried about long-branch attraction in real data sets? Investigations using metazoan 18S rDNA.

Although long-branch attraction (LBA) is frequently cited as the cause of anomalous phylogenetic groupings, few examples of LBA involving real sequence data are known. We have found several cases of probable LBA by analyzing subsamples from an alignment of 18S rDNA sequences for 133 metazoans. In one example, maximum parsimony analysis of sequences from two rotifers, a ctenophore, and a polychaete annelid resulted in strong support for a tree grouping two "long-branch taxa" (a rotifer and the ctenophore). Maximum-likelihood analysis of the same sequences yielded strong support for a more biologically reasonable "rotifer monophyly" tree. Attempts to break up long branches for problematic subsamples through increased taxon sampling reduced, but did not eliminate, LBA problems. Exhaustive analyses of all quartets for a subset of 50 sequences were performed in order to compare the performance of maximum likelihood, equal-weights parsimony, and two additional variants of parsimony; these methods do differ substantially in their rates of failure to recover trees consistent with well established, but highly unresolved phylogenies. Power analyses using simulations suggest that some incorrect inferences by maximum parsimony are due to statistical inconsistency and that when estimates of central branch lengths for certain quartets are very low, maximum-likelihood analyses have difficulty recovering accepted phylogenies even with large amounts of data. These examples demonstrate that LBA problems can occur in real data sets, and they provide an opportunity to investigate causes of incorrect inferences.

Bilaterian phylogeny based on analyses of a region of the sodium-potassium ATPase beta-subunit gene.

Molecular investigations of deep-level relationships within and among the animal phyla have been hampered by a lack of slowly evolving genes that are amenable to study by molecular systematists. To provide new data for use in deep-level metazoan phylogenetic studies, primers were developed to amplify a 1.3-kb region of the alpha subunit of the nuclear-encoded sodium-potassium ATPase gene from 31 bilaterians representing several phyla. Maximum parsimony, maximum likelihood, and Bayesian analyses of these sequences (combined with ATPase sequences for 23 taxa downloaded from GenBank) yield congruent trees that corroborate recent findings based on analyses of other data sets (e.g., the 18S ribosomal RNA gene). The ATPase-based trees support monophyly for several clades (including Lophotrochozoa, a form of Ecdysozoa, Vertebrata, Mollusca, Bivalvia, Gastropoda, Arachnida, Hexapoda, Coleoptera, and Diptera) but do not support monophyly for Deuterostomia, Arthropoda, or Nemertea. Parametric bootstrapping tests reject monophyly for Arthropoda and Nemertea but are unable to reject deuterostome monophyly. Overall, the sodium-potassium ATPase alpha-subunit gene appears to be useful for deep-level studies of metazoan phylogeny.