Search

Background Microbial transmission from parent to offspring is hypothesized to be widespread in vertebrates. However, evidence for this is limited as many evolutionarily important clades remain unexamined. There is currently no data on the microbiota associated with any Chondrichthyan species during embryonic development, despite the global distribution, ecological importance, and phylogenetic position of this clade. In this study, we take the first steps towards filling this gap by investigating the microbiota associated with embryonic development in the little skate, Leucoraja erinacea, a common North Atlantic species and popular system for chondrichthyan biology. Methods To assess the potential for bacterial transmission in an oviparous chondrichthyan, we used 16S rRNA amplicon sequencing to characterize the microbial communities associated with the skin, gill, and egg capsule of the little skate, at six points during ontogeny. Community composition was analyzed using the QIIME2 pipeline and microbial continuity between stages was tracked using FEAST. Results We identify site-specific and stage-specific microbiota dominated by the bacterial phyla Proteobacteria and Bacteroidetes. This composition is similar to, but distinct from, that of previously published data on the adult microbiota of other chondrichthyan species. Our data reveal that the skate egg capsule harbors a highly diverse bacterial community–particularly on the internal surface of the capsule–and facilitates intergenerational microbial transfer to the offspring. Embryonic skin and external gill tissues host similar bacterial communities; the skin and gill communities later diverge as the internal gills and skin denticles develop. Conclusions Our study is the first exploration of the chondrichthyan microbiota throughout ontogeny and provides the first evidence of vertical transmission in this group.

The metazoan 70-kDa heat shock protein (HSP70) family contains several members localized in different subcellular compartments. The cytosolic members have been classified into inducible HSP70s and constitutive heat shock cognates (HSC70s), but their distinction and evolutionary relationship remain unclear because of occasional reports of “constitutive HSP70s” and the lack of cross-phylum comparisons. Here we provide novel insights into the evolution of these important molecular chaperones. Phylogenetic analyses of 125 full-length HSP70s from a broad range of phyla revealed an ancient duplication that gave rise to two lineages from which all metazoan cytosolic HSP70s descend. One lineage (A) contains a relatively small number of genes from many invertebrate phyla, none of which have been shown to be constitutively expressed (i.e., either inducible or unknown). The other lineage (B) included both inducible and constitutive genes from diverse phyla. Species-specific duplications are present in both lineages, and Lineage B contains well-supported phylum-specific clades for Platyhelminthes, Rotifera, Nematoda, Porifera/Cnidaria, and Chordata. Some genes in Lineage B have likely independently acquired inducibility, which may explain the sporadic distribution of “HSP70” or “HSC70” in previous phylogenetic analyses. Consistent with the diversification history within each group, inducible members show lower purifying selection pressure compared to constitutive members. These results illustrate the evolutionary history of the HSP70 family, encouraging us to propose a new nomenclature: “HSP70 + subcellular localization + linage + copy number in the organism + inducible or constitutive, if known.” e.g., HSP70cA1i for cytosolic Lineage A, copy 1, inducible.

To understand the chromatin changes underlying differential gene expression during induced differentiation of human leukemic HL-60/S4 cells, we conducted RNA-Seq analysis on quadruplicate cultures of undifferentiated, granulocytic- and macrophage-differentiated cell forms. More than half of mapped genes exhibited altered transcript levels in the differentiated cell forms. In general, more genes showed increased mRNA levels in the granulocytic form and in the macrophage form, than showed decreased levels. The majority of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were significantly enriched in genes that exhibited differential transcript levels after either RA or TPA treatment. Changes in transcript levels for groups of genes with characteristic protein phenotypes, such as genes encoding cytoplasmic granular proteins, nuclear envelope and cytoskeletal proteins, cell adhesion proteins, and proteins involved in the cell cycle and apoptosis illustrate the profound differences among the various cell states. In addition to the transcriptome analyses, companion karyotyping by M-FISH of undifferentiated HL-60/S4 cells revealed a plethora of chromosome alterations, compared with normal human cells. The present mRNA profiling provides important information related to nuclear shape changes (e.g., granulocyte lobulation), deformability of the nuclear envelope and linkage between the nuclear envelope and cytoskeleton during induced myeloid chromatin differentiation.

In the Anthropocene, in which we now live, climate change is impacting most life on Earth. Microorganisms support the existence of all higher trophic life forms. To understand how humans and other life forms on Earth (including those we are yet to discover) can withstand anthropogenic climate change, it is vital to incorporate knowledge of the microbial ‘unseen majority’. We must learn not just how microorganisms affect climate change (including production and consumption of greenhouse gases) but also how they will be affected by climate change and other human activities. This Consensus Statement documents the central role and global importance of microorganisms in climate change biology. It also puts humanity on notice that the impact of climate change will depend heavily on responses of microorganisms, which are essential for achieving an environmentally sustainable future., Nature Reviews Microbiology, 17 (9), ISSN:1740-1526, ISSN:1740-1534

Understanding patterns and processes in biological diversity is a critical task given current and rapid environmental change. Such knowledge is even more essential when the taxa under consideration are important ecological and evolutionary models. One of these cases is the monogonont rotifer cryptic species complex Brachionus plicatilis, which is by far the most extensively studied group of rotifers, is widely used in aquaculture, and is known to host a large amount of unresolved diversity. Here we collate a dataset of previously available and newly generated sequences of COI and ITS1 for 1273 isolates of the B. plicatilis complex and apply three approaches in DNA taxonomy (i.e. ABGD, PTP, and GMYC) to identify and provide support for the existence of 15 species within the complex. We used these results to explore phylogenetic signal in morphometric and ecological traits, and to understand correlation among the traits using phylogenetic comparative models. Our results support niche conservatism for some traits (e.g. body length) and phylogenetic plasticity for others (e.g. genome size).

A monophyletic origin of endoparasitic thorny-headed worms (Acanthocephala) and wheel-animals (Rotifera) is widely accepted. However, the phylogeny inside the clade, be it called Syndermata or Rotifera, has lacked validation by mitochondrial (mt) data. Herein, we present the first mt genome of the key taxon Seison and report conflicting results of phylogenetic analyses: while mt sequence-based topologies showed monophyletic Lemniscea (Bdelloidea + Acanthocephala), gene order analyses supported monophyly of Pararotatoria (Seisonidea + Acanthocephala) and Hemirotifera (Bdelloidea + Pararotatoria). Sequence-based analyses obviously suffered from substitution saturation, compositional bias, and branch length heterogeneity; however, we observed no compromising effects in gene order analyses. Moreover, gene order-based topologies were robust to changes in coding (genes vs. gene pairs, two-state vs. multistate, aligned vs. non-aligned), tree reconstruction methods, and the treatment of the two monogonont mt genomes. Thus, mt gene order verifies seisonids as sister to acanthocephalans within monophyletic Hemirotifera, while deviating results of sequence-based analyses reflect artificial signal. This conclusion implies that the complex life cycle of extant acanthocephalans evolved from a free-living state, as retained by most monogononts and bdelloids, via an epizoic state with a simple life cycle, as shown by seisonids. Hence, Acanthocephala represent a rare example where ancestral transitional stages have counterparts amongst the closest relatives.

Background Bdelloid rotifers are the oldest, most diverse and successful animal taxon for which males, hermaphrodites, and traditional meiosis are unknown. Their degenerate tetraploid genome, with 2–4 copies of most loci, includes thousands of genes acquired from all domains of life by horizontal transfer. Many bdelloid species thrive in ephemerally aquatic habitats by surviving desiccation at any life stage with no loss of fecundity or lifespan. Their unique genomic diversity and the intense selective pressure of desiccation provide an exceptional opportunity to study the evolution of diversity and novelty in genes involved in DNA repair. Results We used genomic data and RNA-Seq of the desiccation process in the bdelloid Adineta vaga to characterize DNA damage reversal, translesion synthesis, and the major DNA repair pathways: base, nucleotide, and alternate excision repair, mismatch repair (MMR), and double strand break repair by homologous recombination (HR) and classical non-homologous end joining (NHEJ). We identify multiple horizontally transferred DNA damage response genes otherwise unknown in animals (AlkD, Fpg, LigK UVDE), and the presence of genes often considered vertebrate specific, particularly in the NHEJ complex and X family polymerases. While 75–100% of genes involved in MMR and HR are present in 0–2 copies, genes involved in NHEJ, which are present in only a single copy in nearly all other animals, are retained in 3–8 copies. We present structural predictions and expression evidence of neo- or sub-functionalization of multiple copy genes involved in NHEJ and other repair processes. Conclusion The horizontally-acquired genes and duplicated genes in BER and NHEJ suggest resilience to oxidative damage is conferred in part by increased DNA damage recognition and efficient end repair capabilities. The pattern of gene loss and retention in MMR and HR may facilitate recombination and gene conversion between divergent sequences, thus providing at least some of the benefits of sex. The unique retention and divergence of duplicates genes in NHEJ may be facilitated by the lack of efficient selection in the absence of meiotic recombination and independent assortment, and may contribute to the evolutionary success of bdelloids. Electronic supplementary material The online version of this article (10.1186/s12862-018-1288-9) contains supplementary material, which is available to authorized users.

Lifespan extension under low temperature is well conserved across both endothermic and exothermic taxa, but the mechanism underlying this change in aging is poorly understood. Low temperature is thought to decrease metabolic rate, thus slowing the accumulation of cellular damage from reactive oxygen species, although recent evidence suggests involvement of specific cold-sensing biochemical pathways. We tested the effect of low temperature on aging in 11 strains of Brachionus rotifers, with the hypothesis that if the mechanism of lifespan extension is purely thermodynamic, all strains should have a similar increase in lifespan. We found differences in change in median lifespan ranging from a 6% decrease to a 100% increase, as well as differences in maximum and relative lifespan extension and in mortality rate. Low temperature delays reproductive senescence in most strains, suggesting an extension of healthspan, even in strains with little to no change in lifespan. The combination of low temperature and caloric restriction in one strain resulted in an additive lifespan increase, indicating these interventions may work via non- or partially-overlapping pathways. The known low temperature sensor TRPA1 is present in the rotifer genome, but chemical TRPA1 agonists did not affect lifespan, suggesting that this gene may be involved in low temperature sensation but not in chemoreception in rotifers. The congeneric variability in response to low temperature suggests that the mechanism of low temperature lifespan extension is an active genetic process rather than a passive thermodynamic one and is dependent upon genotype.

Despite decades of research, the genetic causes of aging remain elusive. The natural phenotypic and genetic diversity in the Brachionus plicatilis species complex of rotifers can address this challenge. Rotifers are small aquatic invertebrates that have been used in aging studies for nearly 100 years but have only recently emerged as a modern model system. Ease of culturing, small size, and short generation time offer experimental advantages. Caloric restriction (CR), reduced temperature, hypoxia, c-Jun N-terminal kinase or target of rapamycin inhibition, and glycerol supplementation increase life span and in some cases health span, suggesting the physiological and genetic basis of these interventions are conserved between rotifers, established models, and humans. The life cycle of Brachionus produces three types of genetically identical females that have different life span, health span, and fecundity responses to CR interventions, allowing investigation of epigenetic and other maternal effects. There is also variation among closely related members of the B. plicatilis species complex in their response to these interventions, ranging from increases to decreases in life span, health span, and fecundity. Some strains display the classic trade-off between life span and fecundity while others show an increase in both, with implications for the role of the disposable soma theory of aging in different environments. Comparing the response of these natural variants can be a powerful tool to identify genetic pathways modulating health span.

Cell differentiation is associated with changes in chromatin organization and gene expression. In this study, we examine chromatin structure following differentiation of the human myeloid leukemia cell line (HL-60/S4) into granulocytes with retinoic acid (RA) or into macrophage with phorbol ester (TPA). We performed ChIP-seq of histone H3 and its modifications, analyzing changes in nucleosome occupancy, nucleosome repeat length, eu-/heterochromatin redistribution and properties of epichromatin (surface chromatin adjacent to the nuclear envelope). Nucleosome positions changed genome-wide, exhibiting a specific class of alterations involving nucleosome loss in extended (∼1kb) regions, pronounced in enhancers and promoters. Genes that lost nucleosomes at their promoters showed a tendency to be upregulated. On the other hand, nucleosome gain did not show simple effects on transcript levels. The average genome-wide nucleosome repeat length (NRL) did not change significantly with differentiation. However, we detected an approximate 10 bp NRL decrease around the haematopoietic transcription factor (TF) PU.1 and the architectural protein CTCF, suggesting an effect on NRL proximal to TF binding sites. Nucleosome occupancy changed in regions associated with active promoters in differentiated cells, compared with untreated HL-60/S4 cells. Epichromatin regions revealed an increased GC content and high nucleosome density compared with surrounding chromatin. Epichromatin showed depletion of major histone modifications and revealed enrichment with PML body-associated genes. In general, chromatin changes during HL-60/S4 differentiation appeared to be more localized to regulatory regions, compared with genome-wide changes among diverse cell types studied elsewhere.

While many studies have focused on the detrimental effects of advanced maternal age and harmful prenatal environments on progeny, little is known about the role of beneficial non-Mendelian maternal inheritance on aging. Here, we report the effects of maternal age and maternal caloric restriction (CR) on the life span and health span of offspring for a clonal culture of the monogonont rotifer Brachionus manjavacas. Mothers on regimens of chronic CR (CCR) or intermittent fasting (IF) had increased life span compared with mothers fed ad libitum (AL). With increasing maternal age, life span and fecundity of female offspring of AL-fed mothers decreased significantly and life span of male offspring was unchanged, whereas body size of both male and female offspring increased. Maternal CR partially rescued these effects, increasing the mean life span of AL-fed female offspring but not male offspring and increasing the fecundity of AL-fed female offspring compared with offspring of mothers of the same age. Both maternal CR regimens decreased male offspring body size, but only maternal IF decreased body size of female offspring, whereas maternal CCR caused a slight increase. Understanding the genetic and biochemical basis of these different maternal effects on aging may guide effective interventions to improve health span and life span.

A long-standing question in evolutionary biology is how sexual reproduction has persisted in eukaryotic lineages. As cyclical parthenogens, monogonont rotifers are a powerful model for examining this question, yet the molecular nature of sexual reproduction in this lineage is currently understudied. To examine genes involved in meiosis, we generated partial genome assemblies for 2 distantly related monogonont species, Brachionus calyciflorus and B. manja vacas. Here we present an inventory of 89 meiotic genes, of which 80 homologs were identified and annotated from these assemblies. Using phylogenetic analysis, we show that several meiotic genes have undergone relatively recent duplication events that appear to be specific to the monogonont lineage. Further, we compare the expression of “meiosis-specific” genes involved in recombination and all annotated copies of the cell cycle regulatory gene CDC20 between obligate parthenogenetic (OP) and cyclical parthenogenetic (CP) strains of B. caly ciflorus. We show that “meiosis-specific” genes are expressed in both CP and OP strains, whereas the expression of one of the CDC20 genes is specific to cyclical parthenogenesis. The data presented here provide insights into mechanisms of cyclical parthenogenesis and establish expectations for studies of obligate asexual relatives of monogononts, the bdelloid rotifer lineage.

Rotifers are microscopic cosmopolitan zooplankton used as models in ecotoxicological and aging studies due to their several advantages such as short lifespan, ease of culture, and parthenogenesis that enables clonal culture. However, caution is required when measuring their survival time as it is affected by maternal age and maternal feeding conditions. Here we provide a protocol for powerful and reproducible measurement of the survival time in Brachionus rotifers following a careful synchronization of culture conditions over several generations. Empirically, poor synchronization results in early mortality and a gradual decrease in survival rate, thus resulting in weak statistical power. Indeed, under such conditions, calorie restriction (CR) failed to significantly extend the lifespan of B. plicatilis although CR-induced longevity has been demonstrated with well-synchronized rotifer samples in past and present studies. This protocol is probably useful for other invertebrate models, including the fruitfly Drosophila melanogaster and the nematode Caenorhabditis elegans, because maternal age effects have also been reported in these species.

Aim To test whether within-species and among-species patterns of abundance and latitudinal range in marine bacteria resemble those found for macro-organisms, and whether these patterns differ along latitudinal clines. Location Global pelagic marine environments. Methods Taxon-specific sequence abundance and location were retrieved from the open-access V6-rRNA pyrotag sequence data base VAMPS (http://vamps.mbl.edu/), which holds a massive collection of marine bacterial community data sets from the International Census of Marine Microbes sampling effort of global ocean water masses. Data were randomly subsampled to correct for spatial bias and for differences in sampling effort. Results We show that bacterial latitudinal ranges are narrower than expected by chance. When present in both Northern and Southern hemispheres, taxa occupy restricted ranges at similar latitudes on both sides of the equator. A significant and positive relationship exists between sequence abundance and latitudinal range, although this pattern contains a large amount of variance. Abundant taxa in the tropics and in the Northern Hemisphere generally have smaller ranges than those in the Southern Hemisphere. We show that the mean latitudinal range of bacterial taxa increases with latitude, supporting the existence of a Rapoport effect in marine bacterioplankton. Finally, we show that bacterioplankton communities contain a higher proportion of abundant taxa as they approach the poles. Main conclusions Macroecological patterns such as the abundance–range relationship, in general, extend to marine bacteria. However, differences in the shape of these relationships between bacteria and macro-organisms call into question whether the processes and their relative importance in shaping global marine bacteria and macro-organism distributions are the same.

Verrucomicrobia is a bacterial phylum that is commonly detected in soil, but little is known about the distribution and diversity of this phylum in the marine environment. To address this, we analyzed the marine microbial community composition in 506 samples from the International Census of Marine Microbes as well as 11 coastal samples taken from the California Current. These samples from both the water column and sediments covered a wide range of environmental conditions. Verrucomicrobia were present in 98% of the analyzed samples, and thus appeared nearly ubiquitous in the ocean. Based on the occurrence of amplified 16S ribosomal RNA sequences, Verrucomicrobia constituted on average 2% of the water column and 1.4% of the sediment bacterial communities. The diversity of Verrucomicrobia displayed a biogeography at multiple taxonomic levels and thus, specific lineages appeared to have clear habitat preference. We found that subdivision 1 and 4 generally dominated marine bacterial communities, whereas subdivision 2 was more frequent in low salinity waters. Within the subdivisions, Verrucomicrobia community composition were significantly different in the water column compared with sediment as well as within the water column along gradients of salinity, temperature, nitrate, depth and overall water column depth. Although we still know little about the ecophysiology of Verrucomicrobia lineages, the ubiquity of this phylum suggests that it may be important for the biogeochemical cycle of carbon in the ocean.

The mate recognition protein (MRP) gene is a member of a family of extracellular matrix protein genes, called MRP Motif Repeat (MMR) genes, with no known homologs. Two sets of MMR genes, designated MMR-A and MMR-B, were found in Brachionus manjavacas. MMR-B has previously been shown to encode the MRP in the Brachionus plicatilis species complex. MMR family genes share the same basic structure: a signal peptide sequence, followed by nearly identical 276 bp (MMR-A) or 261 bp (MMR-B) repeats, with a truncated final repeat. Each repeat of the predicted MMR-A and -B proteins is expected to have a secondary structure of 5 α-helices, ranging in length from 11 to 20 amino acids, separated by coils of 1–3 amino acids. Hydrophobic and hydrophilic amino acids are predicted to be partitioned to opposite sides of each α-helix, suggesting that MMR proteins are globular with a hydrophobic core. MMR-A and MMR-B proteins vary in their post-translational modifications, resulting in differences in size and charge, and likely causing differences in the physical properties of the proteins on the surface of the female, and their ability to be recognized by a receptor on a male rotifer. The identity of MMR gene repeats is theorized to be maintained by concerted evolution, through a process of unequal crossing over and/or gene conversion, with new mutations likely to be lost. Rarely, however, the same process of concerted evolution can rapidly spread a mutation across all of the repeats. When a mutation results in conformational changes in the protein detectable by males, it could lead to reproductive isolation and thereby to speciation. Thus, changes in MRP could be a driving force in the high degree of species diversity seen within the B. plicatilis cryptic species complex.

Many studies have investigated physiological roles of the membrane associated progesterone receptor (MAPR), but little is known of its evolution. Marked variations in response to exogenous progesterone have been reported for four brachionid rotifer species, suggesting differences in progesterone signaling and reception. Here we report sequence variation for the MAPR gene in the Brachionus plicatilis species complex. Phylogenetic analysis of this receptor is compared with relatedness based on cytochrome c oxidase subunit 1 sequences. Nonsynonymous to synonymous site substitution rate ratios, amino acid divergence, and variations in predicted phosphorylation sites are examined to assess evolution of the MAPR among brachionid clades.

Summary Deep sequencing of PCR amplicon libraries facilitates the detection of low-abundance populations in environmental DNA surveys of complex microbial communities. At the same time, deep sequencing can lead to overestimates of microbial diversity through the generation of low-frequency, error-prone reads. Even with sequencing error rates below 0.005 per nucleotide position, the common method of generating operational taxonomic units (OTUs) by multiple sequence alignment and complete-linkage clustering significantly increases the number of predicted OTUs and inflates richness estimates. We show that a 2% single-linkage preclustering methodology followed by an average-linkage clustering based on pairwise alignments more accurately predicts expected OTUs in both single and pooled template preparations of known taxonomic composition. This new clustering method can reduce the OTU richness in environmental samples by as much as 30‐60% but does not reduce the fraction of OTUs in long-tailed rank abundance curves that defines the rare biosphere.

Meiotic sex evolved early in the history of eukaryotes. Giardia duodenalis (syn. Giardia lamblia, Giardia intestinalis), a parasitic protist belonging to an early diverging lineage of eukaryotes, shows no cytological or physiological evidence of meiotic or sexual processes. Recent molecular analyses challenge the idea that G. duodenalis is a strictly clonal organism by providing evidence of recombination between homologous chromosomes within one subgroup (Assemblage A) of this species as well as genetic transfer from one subgroup to another (Assemblage A-B). Because recombination is not well documented and because it is not known whether the observed inter-assemblage transfer represents true reciprocal genetic exchange or a non-sexual process, we analyzed genic sequences from all major subgroups (Assemblages A-G) of this species. For all assemblages, we detected molecular signatures consistent with meiotic sex or genetic exchange, including low levels of heterozygosity, as indicated by allelic sequence divergence within isolates, and intra- and inter-assemblage recombination. The identification of recombination between assemblages suggests a shared gene pool and calls into question whether it is appropriate to divide the genetically distinct assemblages of G. duodenalis into a species complex.

Giardia intestinalis is a microbial eukaryotic parasite that causes diarrheal disease in humans and other vertebrates worldwide. The negative effect on quality of life and economics caused by G. intestinalis may be increased by its potential status as a zoonosis, or a disease that can be trans- mitted from animals to humans. The zoonotic potential of G. intestinalis has been implied for over 2 decades, with human-infecting genotypes (belonging to the 2 major subgroups, Assemblages A and B) occurring in wildlife and domesticated animals. There are recent reports of G. intestinalis in shell- fish, seals, sea lions and whales, suggesting that marine animals are also potential reservoirs of human disease. However, the prevalence, genetic diversity and effect of G. intestinalis in marine environments and the role that marine animals play in transmission of this parasite to humans are rel- atively unexplored. Here, we provide the first thorough molecular characterization of G. intestinalis in marine vertebrates. Using a multi-locus sequencing approach, we identify human-infecting G. intestinalis haplotypes of both Assemblages A and B in the fecal material of dolphins, porpoises, seals, herring gulls Larus argentatus, common eiders Somateria mollissima and a thresher shark Alopias vulpinus. Our results indicate that G. intestinalis is prevalent in marine ecosystems, and a wide range of marine hosts capable of harboring zoonotic forms of this parasite exist. The presence of G. intestinalis in marine ecosystems raises concerns about how this disease might be transmitted among different host species.

The monogonont rotifer Brachionus plicatilis is an emerging model system for a diverse array of questions in limnological ecosystem dynamics, the evolution of sexual recombination, cryptic speciation, and the phylogeny of basal metazoans. We sequenced the complete mitochondrial genome of B. plicatilis sensu strictu NH1L and found that it is composed of 2 circular chromosomes, designated mtDNA-I (11,153 bp) and mtDNA-II (12,672 bp). Hybridization to DNA isolated from mitochondria demonstrated that mtDNA-I is present at 4 times the copy number of mtDNA-II. The only nucleotide similarity between the 2 chromosomes is a 4.9-kbp region of 99.5% identity including a transfer RNA (tRNA) gene and an extensive noncoding region that contains putative D-loop and control sequence. The mtDNA-I chromosome encodes 4 proteins (ATP6, COB, NAD1, and NAD2), 13 tRNAs, and the large and small subunit ribosomal RNAs; mtDNA-II encodes 8 proteins (COX1-3, NAD3-6, and NAD4L) and 9 tRNAs. Gene order is not conserved between B. plicatilis and its closest relative with a sequenced mitochondrial genome, the acanthocephalan Leptorhynchoides thecatus, or other sequenced mitochondrial genomes. Polymerase chain reaction assays and Southern hybridization to DNA from 18 strains of Brachionus suggest that the 2-chromosome structure has been stable for millions of years. The novel organization of the B. plicatilis mitochondrial genome into 2 nearly equal chromosomes of 4-fold different copy number may provide insight into the evolution of metazoan mitochondria and the phylogenetics of rotifers and other basal animal phyla.

The evolution of marine microbes over billions of years predicts that the composition of microbial communities should be much greater than the published estimates of a few thousand distinct kinds of microbes per liter of seawater. By adopting a massively parallel tag sequencing strategy, we show that bacterial communities of deep water masses of the North Atlantic and diffuse flow hydrothermal vents are one to two orders of magnitude more complex than previously reported for any microbial environment. A relatively small number of different populations dominate all samples, but thousands of low-abundance populations account for most of the observed phylogenetic diversity. This >rare biosphere> is very ancient and may represent a nearly inexhaustible source of genomic innovation. Members of the rare biosphere are highly divergent from each other and, at different times in earth's history, may have had a profound impact on shaping planetary processes. © 2006 by The National Academy of Sciences of the USA.

Rotifera and Acanthocephala are generally regarded as separate phyla sharing a basal position among triploblast protostomes. This paper presents the first molecular phylogenetic examination of the relationship of Acanthocephala to all three rotifer classes, Seisonidea, Monogononta, and Bdelloidea. Inclusion of Acanthocephala within Rotifera, probably as a sister-taxon to a clade composed of Bdelloidea and Monogononta (the Eurotatoria), is strongly supported by both parsimony and distance methods, using a region of the nuclear coding gene hsp82. Previous molecular evidence for the inclusion of Acanthocephala in the Rotifera suggested that Acanthocephala is a sister-taxon of Bdelloidea, forming the clade Lemniscea. No support is found for this clade, and evidence is presented that the monogonont rotifer used in those analyses, Brachionus plicatilis, may be evolving in an anomalous manner.

Rotifers of the asexual class Bdelloidea are unusual in possessing two or more divergent copies of every gene that has been examined. Phylogenetic analysis of the heat-shock gene hsp82 and the TATA-box-binding protein gene tbp in multiple bdelloid species suggested that for each gene, each copy belonged to one of two lineages that began to diverge before the bdelloid radiation. Such gene trees are consistent with the two lineages having descended from former alleles that began to diverge after meiotic segregation ceased or from subgenomes of an alloploid ancestor of the bdelloids. However, the original analyses of bdelloid gene-copy divergence used only a single outgroup species and were based on parsimony and neighbor joining. We have now used maximum likelihood and Bayesian inference methods and, for hsp82 , multiple outgroups in an attempt to produce more robust gene trees. Here we report that the available data do not unambiguously discriminate between gene trees that root the origin of hsp82 and tbp copy divergence before the bdelloid radiation and those which indicate that the gene copies began to diverge within bdelloid families. The remarkable presence of multiple diverged gene copies in individual genomes is nevertheless consistent with the loss of sex in an ancient ancestor of bdelloids.

The approximately 360 species of the Class Bdelloidea (Phylum Rotifera) constitute the largest, most diverse and oldest eukaryotic group for which obligate asexuality has morphological, cytological and molecular support. DNA sequence studies indicate that the class constitutes a single clonal lineage, 50-100 million years old. Here we report cytofluorometric measurements of oocyte nuclear DNA content and GC proportion for six bdelloid species, representing the three major families of Bdelloidea. We find wide differences between species in nuclear DNA content (0.5-2.4 pg) and proportion of GC (37-57%). The proportion of GC in the nuclear DNA of all six species is greater than that in their coding regions examined to date, and tends to be greater in species with greater nuclear DNA content. Interspecific differences in nuclear DNA content do not correlate with chromosome numbers, which range from 10 to 13, and are therefore not due to differences in ploidy. In addition, since bdelloids seemingly lack functional reverse transcriptase genes, the observed interspecific differences in nuclear DNA content and GC proportion do not appear to have resulted from retrotransposon activity. The values of oocyte nuclear DNA content obtained by cytofluorometric measurement are close to the values of total genome size we previously obtained from measurements of DNA hybridization intensity, indicating that bdelloid genomes do not undergo significant somatic chromatin diminution and, unlike the genomes of sexually reproducing diploids or polyploids, do not consist of two or more nearly identical haplotypes.

The human mouth is an excellent system to study the dynamics of microbial communities and their interactions with their host. We employed oligotyping to analyze, with single-nucleotide resolution, oral microbial 16S ribosomal RNA (rRNA) gene sequence data from a time course sampled from the tongue of two individuals, and we interpret our results in the context of oligotypes that we previously identified in the oral data from the Human Microbiome Project. Our previous work established that many of these oligotypes had dramatically different distributions between individuals and across oral habitats, suggesting that they represented functionally different organisms. Here we demonstrate the presence of a consistent tongue microbiome but with rapidly fluctuating proportions of the characteristic taxa. In some cases closely related oligotypes representing strains or variants within a single species displayed fluctuating relative abundances over time, while in other cases an initially dominant oligotype was replaced by another oligotype of the same species. We use this high temporal and taxonomic level of resolution to detect correlated changes in oligotype abundance that could indicate which taxa likely interact synergistically or occupy similar habitats, and which likely interact antagonistically or prefer distinct habitats. For example, we found a strong correlation in abundance over time between two oligotypes from different families of Gamma Proteobacteria, suggesting a close functional or ecological relationship between them. In summary, the tongue is colonized by a microbial community of moderate complexity whose proportional abundance fluctuates widely on time scales of days. The drivers and functional consequences of these community dynamics are not known, but we expect they will prove tractable to future, targeted studies employing taxonomically resolved analysis of high-throughput sequencing data sampled at appropriate temporal intervals and spatial scales.

Loss of sexual reproduction is considered an evolutionary dead end for metazoans, but bdelloid rotifers challenge this view as they appear to have persisted asexually for millions of years1. Neither male sex organs nor meiosis have ever been observed in these microscopic animals: oocytes are formed through mitotic divisions, with no reduction of chromosome number and no indication of chromosome pairing2. However, current evidence does not exclude that they may engage in sex on rare, cryptic occasions. Here we report the genome of a bdelloid rotifer, Adineta vaga (Davis, 1873)3, and show that its structure is incompatible with conventional meiosis. At gene scale, the genome of A. vaga is tetraploid and comprises both anciently duplicated segments and less divergent allelic regions. However, in contrast to sexual species, the allelic regions are rearranged and sometimes even found on the same chromosome. Such structure does not allow meiotic pairing; instead, we find abundant evidence of gene conversion, which may limit the accumulation of deleterious mutations in the absence of meiosis. Gene families involved in resistance to oxidation, carbohydrate metabolism and defence against transposons are significantly expanded, which may explain why transposable elements cover only 3% of the assembled sequence. Furthermore, 8% of the genes are likely to be of non-metazoan origin and were probably acquired horizontally. This apparent convergence between bdelloids and prokaryotes sheds new light on the evolutionary significance of sex., National Science Foundation (U.S.) (grant MCB-0821956), National Science Foundation (U.S.) (NSF grant MCB-1121334), German Research Foundation (grant HA 5163/2-1), Ministry of Education and Science of the Russian Federation (grant 11.G34.31.0008), National Science Foundation (U.S.) (NSF CAREER grant number 0644282), National Science Foundation (U.S.) (NSF grant MCB-0923676), Belgian National Foundation for Scientific Research (FRFC grant 2.4.655.09.F), University of Namur (start-up grant)

Caloric restriction (CR) is cited as the most robust means of increasing lifespan across a range of taxa, yet there is a high degree of variability in the response to CR, both within and between species. To examine the intraspecific evolutionary conservation of lifespan extension by CR, we tested the effects of chronic caloric restriction (CCR) at multiple food levels and of intermittent fasting (IF) in twelve isolates from the Brachionus plicatilis species complex of monogonont rotifers. While CCR generally increased or did not change lifespan and total fecundity, IF caused increased, unchanged, or decreased lifespan, depending upon the isolate, and decreased total fecundity in all but one isolate. Lifespan under ad libitum (AL) feeding varied among isolates and predicted the lifespan response to CR: longer-lived isolates under AL were less likely to have a significant increase in lifespan under CCR and were more likely to have a significantly shortened lifespan under IF. Lifespan under AL conditions and the response to CR were not correlated with hydroperiodicity of native habitat or with time in culture. Lack of trade-off between lifespan and fecundity under CCR, and differences in lifespan and fecundity under CCR and IF, even when average food intake was similar, suggest that longevity changes are not always directly determined by energy intake and that CCR and IF regimens extend lifespan through diverse genetic mechanisms.

We measured life span and fecundity of three reproductive modes in a clone of the monogonont rotifer Brachionus manjavacas subjected to chronic caloric restriction (CCR) over a range of food concentrations or to intermittent fasting (IF). IF increased life span 50%-70% for all three modes, whereas CCR increased life span of asexual females derived from sexually or asexually produced eggs, but not that of sexual females. The main effect of CR on both asexual modes was to delay death at young ages, rather than to prevent death at middle ages or to greatly extend maximum life span; in contrast CR in sexual females greatly increased the life span of a few long-lived individuals. Lifetime fecundity did not decrease with CCR, suggesting a lack of resource allocation trade-off between somatic maintenance and reproduction. Multiple outcomes for a clonal lineage indicate that different responses are established through epigenetic programming, whereas differences in life-span allocations suggest that multiple genetic mechanisms mediate life-span extension.

Spatial patterns of marine Synechococcus diversity across ocean domains have been reported on extensively. However, much less is known of seasonal and multiannual patterns of change in Synechococcus community composition. Here we report on the genotypic diversity of Synechococcus populations in the Gulf of Aqaba, Northern Red Sea, over seven annual cycles of deep mixing and stabile stratification, using ntcA as a phylogenetic marker. Synechococcus clone libraries were dominated by clade II and XII genotypes and a total of eight different clades were identified. Inclusion of ntcA sequences from the Global Ocean Sampling database in our analyses identified members of clade XII from beyond the Gulf of Aqaba, extending its known distribution. Most of the Synechococcus diversity was attributed to members of clade II during the spring bloom, while clade III contributed significantly to diversity during summer stratification. Clade XII diversity was most prevalent in fall and winter. Clade abundances were estimated from pyrosequencing of the V6 hypervariable region of 16S rRNA. Members of clade II dominated Synechococcus communities throughout the year, whereas the less frequent genotypes showed a pattern of seasonal succession. Based on the prevailing nutritional conditions we observed that clade I members thrive at higher nutrient concentrations during winter mixing. Clades V, VI and X became apparent during the transition periods between mixing and stratification. Clade III became prominent during sumeer stratification. We propose that members of clades V, VI, and X, and clade III are Synechococcus ecotypes that are adapted to intermediate and low nutrient levels respectively. This is the first time that molecular analyses have correlated population dynamics of Synechococcus genotypes with temporal fluctuations in nutrient regimes. Since these Synechococcus genotypes are routinely observed in the Gulf of Aqaba we suggest that seasonal fluctuations in nutrient levels create temporal niches that sustain their coexistence.

Background Populations of Atlantic killifish (Fundulus heteroclitus) have evolved resistance to the embryotoxic effects of polychlorinated biphenyls (PCBs) and other halogenated and nonhalogenated aromatic hydrocarbons that act through an aryl hydrocarbon receptor (AHR)-dependent signaling pathway. The resistance is accompanied by reduced sensitivity to induction of cytochrome P450 1A (CYP1A), a widely used biomarker of aromatic hydrocarbon exposure and effect, but whether the reduced sensitivity is specific to CYP1A or reflects a genome-wide reduction in responsiveness to all AHR-mediated changes in gene expression is unknown. We compared gene expression profiles and the response to 3,3',4,4',5-pentachlorobiphenyl (PCB-126) exposure in embryos (5 and 10 dpf) and larvae (15 dpf) from F. heteroclitus populations inhabiting the New Bedford Harbor, Massachusetts (NBH) Superfund site (PCB-resistant) and a reference site, Scorton Creek, Massachusetts (SC; PCB-sensitive). Results Analysis using a 7,000-gene cDNA array revealed striking differences in responsiveness to PCB-126 between the populations; the differences occur at all three stages examined. There was a sizeable set of PCB-responsive genes in the sensitive SC population, a much smaller set of PCB-responsive genes in NBH fish, and few similarities in PCB-responsive genes between the two populations. Most of the array results were confirmed, and additional PCB-regulated genes identified, by RNA-Seq (deep pyrosequencing). Conclusions The results suggest that NBH fish possess a gene regulatory defect that is not specific to one target gene such as CYP1A but rather lies in a regulatory pathway that controls the transcriptional response of multiple genes to PCB exposure. The results are consistent with genome-wide disruption of AHR-dependent signaling in NBH fish.

BackgroundMarine microbial communities have been essential contributors to global biomass, nutrient cycling, and biodiversity since the early history of Earth, but so far their community distribution patterns remain unknown in most marine ecosystems.Methodology/principal findingsThe synthesis of 9.6 million bacterial V6-rRNA amplicons for 509 samples that span the global ocean's surface to the deep-sea floor shows that pelagic and benthic communities greatly differ, at all taxonomic levels, and share Conclusions/significanceThis first synthesis of global bacterial distribution across different ecosystems of the World's oceans shows remarkable horizontal and vertical large-scale patterns in bacterial communities. This opens interesting perspectives for the definition of biogeographical biomes for bacteria of ocean waters and the seabed.

Low-temperature hydrothermal vent fluids represent access points to diverse microbial communities living in oceanic crust. This study examined the distribution, relative abundance, and diversity of Epsilonproteobacteria in 14 low-temperature vent fluids from five volcanically active seamounts of the Mariana Arc using a 454 tag sequencing approach. Most vent fluids were enriched in cell concentrations compared with background seawater, and quantitative PCR results indicated that all fluids were dominated by bacteria. Operational taxonomic unit-based statistical tools applied to 454 data show that all vents from the northern end of the Mariana Arc grouped together, to the exclusion of southern arc seamounts, which were as distinct from one another as they were from northern seamounts. Statistical analysis also showed a significant relationship between seamount and individual vent groupings, suggesting that community membership may be linked to geographical isolation and not geochemical parameters. However, while there may be large-scale geographic differences, distance is not the distinguishing factor in the microbial community composition. At the local scale, most vents host a distinct population of Epsilonproteobacteria, regardless of seamount location. This suggests that there may be barriers to exchange and dispersal for these vent endemic microorganisms at hydrothermal seamounts of the Mariana Arc.

Giardia duodenalis is an intestinal parasite of many vertebrates. The presence of G. duodenalis in the marine environment due to anthropogenic and wildlife activity is well documented, including the contributions from untreated sewage and storm water, agricultural run-off and droppings from terrestrial animals. Recently, studies have detected this protistan parasite in the faeces of marine vertebrates such as whales, dolphins, seals and shore birds. To explore the population biology of G. duodenalis in marine life, we determined the prevalence of G. duodenalis in two species of seal (Halichoerus grypus, Phoca vitulina vitulina and Phoca vitulina richardsi) from the east and west coasts of the USA, sequenced two loci from G. duodenalis-positive samples to assess molecular diversity and examined G. duodenalis distribution amongst these seals and other marine vertebrates along the east coast. We found a significant difference in the presence of G. duodenalis between east and west coast seal species. Only the zoonotic lineages of G. duodenalis, Assemblages A and B and a novel lineage, which we designated as Assemblage H, were identified in marine vertebrates. Assemblages A and B are broadly distributed geographically and show a lack of host specificity. Only grey seal (Halichoerus grypus) samples and one gull sample (Larus argentatus) from a northern location of Cape Cod, Massachusetts, USA, showed the presence of Assemblage H haplotypes; only one other study of harbour seals from the Puget Sound region of Washington, USA previously recorded the presence of an Assemblage H haplotype. Assemblage H sequences form a monophyletic clade that appears as divergent from the other seven Assemblages of G. duodenalis as these assemblages are from each other. The discovery of a previously uncharacterised lineage of G. duodenalis suggests that this parasite has more genetic diversity and perhaps a larger host range than previously believed.

Lamin B receptor (LBR) is an integral membrane protein of the interphase nuclear envelope (NE). The N-terminal end resides in the nucleoplasm, binding to lamin B and heterochromatin, with the interactions disrupted during mitosis. The C-terminal end resides within the inner nuclear membrane, retreating with the ER away from condensing chromosomes during mitotic NE breakdown. Some of these properties are interpretable in terms of our current structural knowledge of LBR, but many of the structural features remain unknown. LBR apparently has an evolutionary history which brought together at least two ancient conserved structural domains (i.e., Tudor and sterol reductase). This convergence may have occurred with the emergence of the chordates and echinoderms. It is not clear what survival values have maintained LBR structure during evolution. But it seems likely that roles in post-mitotic nuclear reformation, interphase NE growth and compartmentalization of nuclear architecture might have provided some evolutionary advantage to preservation of the LBR gene. REVIEW

Background Mate choice is of central importance to most animals, influencing population structure, speciation, and ultimately the survival of a species. Mating behavior of male brachionid rotifers is triggered by the product of a chemosensory gene, a glycoprotein on the body surface of females called the mate recognition pheromone. The mate recognition pheromone has been biochemically characterized, but little was known about the gene(s). We describe the isolation and characterization of the mate recognition pheromone gene through protein purification, N-terminal amino acid sequence determination, identification of the mate recognition pheromone gene from a cDNA library, sequencing, and RNAi knockdown to confirm the functional role of the mate recognition pheromone gene in rotifer mating. Results A 29 kD protein capable of eliciting rotifer male circling was isolated by high-performance liquid chromatography. Two transcript types containing the N-terminal sequence were identified in a cDNA library; further characterization by screening a genomic library and by polymerase chain reaction revealed two genes belonging to each type. Each gene begins with a signal peptide region followed by nearly perfect repeats of an 87 to 92 codon motif with no codons between repeats and the final motif prematurely terminated by the stop codon. The two Type A genes contain four and seven repeats and the two Type B genes contain three and five repeats, respectively. Only the Type B gene with three repeats encodes a peptide with a molecular weight of 29 kD. Each repeat of the Type B gene products contains three asparagines as potential sites for N-glycosylation; there are no asparagines in the Type A genes. RNAi with Type A double-stranded RNA did not result in less circling than in the phosphate-buffered saline control, but transfection with Type B double-stranded RNA significantly reduced male circling by 17%. The very low divergence between repeat units, even at synonymous positions, suggests that the repeats are kept nearly identical through a process of concerted evolution. Information-rich molecules like surface glycoproteins are well adapted for chemical communication and aquatic animals may have evolved signaling systems based on these compounds, whereas insects use cuticular hydrocarbons. Conclusion Owing to its critical role in mating, the mate recognition pheromone gene will be a useful molecular marker for exploring the mechanisms and rates of selection and the evolution of reproductive isolation and speciation using rotifers as a model system. The phylogenetic variation in the mate recognition pheromone gene can now be studied in conjunction with the large amount of ecological and population genetic data being gathered for the Brachionus plicatilis species complex to understand better the evolutionary drivers of cryptic speciation.

PCR-based surveys of microbial communities commonly use regions of the small-subunit ribosomal RNA (SSU rRNA) gene to determine taxonomic membership and estimate total diversity. Here we show that the length of the target amplicon has a significant effect on assessments of microbial richness and community membership. Using operational taxonomic unit (OTU)- and taxonomy-based tools, we compared the V6 hypervariable region of the bacterial SSU rRNA gene of three amplicon libraries of c. 100, 400 and 1000 base pairs (bp) from each of two hydrothermal vent fluid samples. We found that the smallest amplicon libraries contained more unique sequences, higher diversity estimates and a different community structure than the other two libraries from each sample. We hypothesize that a combination of polymerase dissociation, cloning bias and mispriming due to secondary structure accounts for the differences. While this relationship is not linear, it is clear that the smallest amplicon libraries contained more different types of sequences, and accordingly, more diverse members of the community. Because divergent and lower abundant taxa can be more readily detected with smaller amplicons, they may provide better assessments of total community diversity and taxonomic membership than longer amplicons in molecular studies of microbial communities.

Despite considerable attention, what drives asexual populations to extinction and allows the maintenance of sexual reproduction remains a mystery. Bdelloid rotifers appear to be an exception to the general rule that multicellular eukaryotes that abandon sex are doomed to early extinction, and an understanding of their biology may reveal how they have avoided the fate of other asexual lineages. Here, we discuss aspects of the ecology, physiology, and genome evolution of bdelloid rotifers that may contribute to their evolutionary success. We propose that many of the unusual characteristics of bdelloids derive from their adaptation to desiccation-prone habitats but may also contribute to their long-term success in the absence of sexual reproduction. However, many unanswered questions about these characteristics and the basic biology of bdelloids remain. Population genetics and molecular ecology are promising approaches for the next generation of studies using this model system to provide answers to the questions of how and why sex is the dominant form of reproduction in plants and animals.

The fourth author's name was listed incorrectly in the article's citation. The correct citation is: Huse SM, Dethlefsen L, Huber JA, Mark Welch D, Relman DA, et al. (2008) Exploring Microbial Diversity and Taxonomy Using SSU rRNA Hypervariable Tag Sequencing. PLoS Genet 4(11): e1000255. doi:10.1371/journal.pgen.1000255

Rotifers of class Bdelloidea have evolved for millions of years apparently without sexual reproduction. We have sequenced 45- to 70-kb regions surrounding the four copies of the hsp82 gene of the bdelloid rotifer Philodina roseola , each of which is on a separate chromosome. The four regions comprise two colinear gene-rich pairs with gene content, order, and orientation conserved within each pair. Only a minority of genes are common to both pairs, also in the same orientation and order, but separated by gene-rich segments present in only one or the other pair. The pattern is consistent with degenerate tetraploidy with numerous segmental deletions, some in one pair of colinear chromosomes and some in the other. Divergence in 1,000-bp windows varies along an alignment of a colinear pair, from zero to as much as 20% in a pattern consistent with gene conversion associated with recombinational repair of DNA double-strand breaks. Although pairs of colinear chromosomes are a characteristic of sexually reproducing diploids and polyploids, a quite different explanation for their presence in bdelloids is suggested by the recent finding that bdelloid rotifers can recover and resume reproduction after suffering hundreds of radiation-induced DNA double-strand breaks per oocyte nucleus. Because bdelloid primary oocytes are in G 1 and therefore lack sister chromatids, we propose that bdelloid colinear chromosome pairs are maintained as templates for the repair of DNA double-strand breaks caused by the frequent desiccation and rehydration characteristic of bdelloid habitats.

The analytical power of environmental DNA sequences for modeling microbial ecosystems depends on accurate assessments of population structure, including diversity (richness) and relative abundance (evenness). We investigated both aspects of population structure for microbial communities at two neighboring hydrothermal vents by examining the sequences of more than 900,000 microbial small-subunit ribosomal RNA amplicons. The two vent communities have different population structures that reflect local geochemical regimes. Descriptions of archaeal diversity were nearly exhaustive, but despite collecting an unparalleled number of sequences, statistical analyses indicated additional bacterial diversity at every taxonomic level. We predict that hundreds of thousands of sequences will be necessary to capture the vast diversity of microbial communities, and that different patterns of evenness for both high- and low-abundance taxa may be important in defining microbial ecosystem dynamics.

Rotifers are a key component of many freshwater ecosystems, but surveys of the composition of rotifer communities are limited by the labor-intensiveness of sample processing, particularly of non-planktonic taxa, and by the shortage of investigators qualified to identify a broad range of rotifer species. Additional problems are posed by species that must be identified from living specimens, and by members of cryptic species complexes. As DNA sequencing becomes easier and cheaper, it has become practical to obtain representative DNA sequences from identified rotifer species for use in genome-based surveys to determine which rotifers are present in a new sample, avoiding the difficulties of traditional surveys. Here we discuss two genome-based tools used in surveys of microbial communities: serial analysis of gene tags (SAGT) and microarray hybridization. SAGT is a method for inexpensively obtaining characteristic short DNA sequences from a sample that can both identify taxa for which the tag sequence is known and signal the presence of additional uncharacterized species. Microarray hybridization allows detection of DNA sequences in the sample that are identical or similar to sequences present on the microarray. We also report the construction and hybridization of a small microarray of rotifer sequences, demonstrating that this method can discriminate among bdelloid families, and is likely to make much finer discriminations if appropriate sequences are present on the microarray. These techniques are most powerful when combined with traditional systematics in collaborative efforts, which may be fostered through the data base of rotifer biology, WheelBase (http://jbpc.mbl.edu/wheelbase).

Rotifera is composed of groups with unusual ultrastructural, physiological, and reproductive characters. Our ability to understand the evolution of these features is complicated by the fact that the phylogenetic relationships among the three traditional rotifer groups (Seisonidea, Monogononta, and Bdelloidea) and Acanthocephala remain unresolved. Here, I present maximum likelihood and Bayesian analyses of rotifer-acanthocephalan relationships using both the protein-coding gene hsp82 and a combined data set of hsp82 and ribosomal small subunit (SSU) DNA sequences, using nucleotide and codon based models of evolution. Statistical analysis of the phylogenetic support for any of the likely relationships among rotifer groups suggests that more than a combined hsp82 + SSU data set will be needed to resolve rotifer-acanthocephalan phylogeny with any degree of certainty.

DNA sequencing has shown individual bdelloid rotifer genomes to contain two or more diverged copies of every gene examined and has revealed no closely similar copies. These and other findings are consistent with long-term asexual evolution of bdelloids. It is not entirely ruled out, however, that bdelloid genomes consist of previously undetected pairs of sequences so similar as to be identical over the regions sequenced, as might result if bdelloids were highly inbred sexual diploids or polyploids. Here, we employ fluorescent in situ hybridization with cosmid probes to determine the copy number and chromosomal distribution of the heat shock gene hsp82 and adjacent sequences in the bdelloid Philodina roseola . We conclude that the four copies identified by sequencing are the only ones present and that each is on a separate chromosome. Bdelloids therefore are not highly homozygous sexually reproducing diploids or polyploids.

The class Bdelloidea of the phylum Rotifera is the largest well studied eukaryotic taxon in which males and meiosis are unknown, and the only one for which these indications of ancient asexuality are supported by cytological and molecular genetic evidence. We estimated the rates of synonymous and nonsynonymous substitutions in the hsp82 heat shock gene in bdelloids and in facultatively sexual rotifers of the class Monogononta, employing distance based and maximum likelihood methods. Relative-rate tests, using acanthocephalan rotifers as an outgroup, showed slightly higher rates of nonsynonymous substitution and slightly lower rates of synonymous substitution in bdelloids as compared with monogononts. The opposite trend, however, was seen in intraclass pairwise comparisons. If, as it seems, bdelloids have evolved asexually, an equality of bdelloid and monogonont substitution rates would suggest that the maintenance of sexual reproduction in monogononts is not attributable to an effect of sexual reproduction in limiting the load of deleterious nucleotide substitutions.

Here we report on the occurrence and position of introns found in three genes of rotifers. A region of the gene for the TATA-box binding protein was examined in three species of Bdelloidea and one of Monogononta. There are two introns in both copies of this gene present in each of the three bdelloids examined — one at a position where introns occur in other eukaryotes and the other at a novel position; the monogonont has no introns in the region examined. A region of the gene encoding the 82 kD heat shock protein was examined in 10 species, with every rotifer class represented. Introns were found in only two species, both bdelloids: one of the species has an intron in all three copies of the gene; the other has an intron in only one of the three copies. Both introns occur at novel positions. The gene for triosephosphate isomerase was examined in one bdelloid. Both copies of the gene in this species contain introns, all at conserved positions: one copy contains five introns, the other copy three. These observations demonstrate the presence of introns in bdelloid rotifers, some in conserved positions, others apparently newly arisen during bdelloid evolution.

The past decade has seen the application of DNA sequence data to phylogenetic investigations of Rotifera, both expanding and challenging our understanding of the evolution of the phylum. Evidence that Acanthocephala, long regarded as a separate but closely related phylum, is a highly derived class of Rotifera demonstrates the potential of molecular analyses to suggest relationships not obvious from morphological analysis. Phylogenies based on the sequence of the gene for the small ribosomal RNA suggest that rotifers and acanthocephalans are associated with Platyhelminthes and Gastrotricha, perhaps in a clade with Gnathostomula and Cycliophora; at present, this group lacks a clear morphological synapomorphy. A more complete resolution of the molecular phylogeny of Rotifera will require surveying multiple genes and several species from each clade under investigation.