Search

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.

As microbiome science expands, academic centres scramble to fill many needs, from service provider to industry liaison. A newly created network aims to share strategies and accelerate knowledge transfer, and invites others to join the efforts.

Background The causes and consequences of genome size variation across Eukaryotes, which spans five orders of magnitude, have been hotly debated since before the advent of genome sequencing. Previous studies have mostly examined variation among larger taxonomic units (e.g., orders, or genera), while comparisons among closely related species are rare. Rotifers of the Brachionus plicatilis species complex exhibit a seven-fold variation in genome size and thus represent a unique opportunity to study such changes on a relatively short evolutionary timescale. Here, we sequenced and analysed the genomes of four species of this complex with nuclear DNA contents spanning 110–422 Mbp. To establish the likely mechanisms of genome size change, we analysed both sequencing read libraries and assemblies for signatures of polyploidy and repetitive element content. We also compared these genomes to that of B. calyciflorus, the closest relative with a sequenced genome (293 Mbp nuclear DNA content). Results Despite the very large differences in genome size, we saw no evidence of ploidy level changes across the B. plicatilis complex. However, repetitive element content explained a large portion of genome size variation (at least 54%). The species with the largest genome, B. asplanchnoidis, has a strikingly high 44% repetitive element content, while the smaller B. plicatilis genomes contain between 14 and 25% repetitive elements. According to our analyses, the B. calyciflorus genome contains 39% repetitive elements, which is substantially higher than previously reported (21%), and suggests that high repetitive element load could be widespread in monogonont rotifers. Conclusions Even though the genome sizes of these species are at the low end of the metazoan spectrum, their genomes contain substantial amounts of repetitive elements. Polyploidy does not appear to play a role in genome size variations in these species, and these variations can be mostly explained by changes in repetitive element content. This contradicts the naïve expectation that small genomes are streamlined, or less complex, and that large variations in nuclear DNA content between closely related species are due to polyploidy. Electronic supplementary material The online version of this article (10.1186/s12864-019-5859-y) contains supplementary material, which is available to authorized users.

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.

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.

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.

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.

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.

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

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.

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.

The Class Bdelloidea of the Phylum Rotifera is the largest metazoan taxon in which males, hermaphrodites, and meiosis are unknown. We conducted a molecular genetic test of this indication that bdelloid rotifers may have evolved without sexual reproduction or genetic exchange. The test is based on the expectation that after millions of years without these processes, genomes will no longer contain pairs of closely similar haplotypes and instead will contain highly divergent descendants of formerly allelic nucleotide sequences. We find that genomes of individual bdelloid rotifers, representing four different species, appear to lack pairs of closely similar sequences and contain representatives of two ancient lineages that began to diverge before the bdelloid radiation many millions of years ago when sexual reproduction and genetic exchange may have ceased.

Background Chemically mediated prezygotic barriers to reproduction likely play an important role in speciation. In facultatively sexual monogonont rotifers from the Brachionus plicatilis cryptic species complex, mate recognition of females by males is mediated by the Mate Recognition Protein (MRP), a globular glycoprotein on the surface of females, encoded by the mmr-b gene family. In this study, we sequenced mmr-b copies from 27 isolates representing 11 phylotypes of the B. plicatilis species complex, examined the mode of evolution and selection of mmr-b, and determined the relationship between mmr-b genetic distance and mate recognition among isolates. Results Isolates of the B. plicatilis species complex have 1–4 copies of mmr-b, each composed of 2–9 nearly identical tandem repeats. The repeats within a gene copy are generally more similar than are gene copies among phylotypes, suggesting concerted evolution. Compared to housekeeping genes from the same isolates, mmr-b has accumulated only half as many synonymous differences but twice as many non-synonymous differences. Most of the amino acid differences between repeats appear to occur on the outer face of the protein, and these often result in changes in predicted patterns of phosphorylation. However, we found no evidence of positive selection driving these differences. Isolates with the most divergent copies were unable to mate with other isolates and rarely self-crossed. Overall the degree of mate recognition was significantly correlated with the genetic distance of mmr-b. Conclusions Discrimination of compatible mates in the B. plicatilis species complex is determined by proteins encoded by closely related copies of a single gene, mmr-b. While concerted evolution of the tandem repeats in mmr-b may function to maintain identity, it can also lead to the rapid spread of a mutation through all copies in the genome and thus to reproductive isolation. The mmr-b gene is evolving rapidly, and novel alleles may be maintained and increase in frequency via asexual reproduction. Our analyses indicate that mate recognition, controlled by MMR-B, may drive reproductive isolation and allow saltational sympatric speciation within the B. plicatilis cryptic species complex, and that this process may be largely neutral.

Background. Rotifers are among the most common non-arthropod animals and are the most experimentally tractable members of the basal assemblage of metazoan phyla known as Gnathifera. The monogonont rotifer Brachionus plicatilis is a developing model system for ecotoxicology, aquatic ecology, cryptic speciation, and the evolution of sex, and is an important food source for finfish aquaculture. However, basic knowledge of the genome and transcriptome of any rotifer species has been lacking. Methodology/Principal Findings. We generated and partially sequenced a cDNA library from B. plicatilis and constructed a database of over 2300 expressed sequence tags corresponding to more than 450 transcripts. About 20% of the transcripts had no significant similarity to database sequences by BLAST; most of these contained open reading frames of significant length but few had recognized Pfam motifs. Sixteen transcripts accounted for 25% of the ESTs; four of these had no significant similarity to BLAST or Pfam databases. Putative up- and downstream untranslated regions are relatively short and AT rich. In contrast to bdelloid rotifers, there was no evidence of a conserved trans-spliced leader sequence among the transcripts and most genes were single-copy. Conclusions/Significance. Despite the small size of this EST project it revealed several important features of the rotifer transcriptome and of individual monogonont genes. Because there is little genomic data for Gnathifera, the transcripts we found with no known function may represent genes that are species-, class-, phylum- or even superphylum-specific; the fact that some are among the most highly expressed indicates their importance. The absence of trans-spliced leader exons in this monogonont species contrasts with their abundance in bdelloid rotifers and indicates that the presence of this phenomenon can vary at the subphylum level. Our EST database provides a relatively large quantity of transcript-level data for B. plicatilis, and more generally of rotifers and other gnathiferan phyla, and can be browsed and searched at gmod.mbl.edu., PLoS ONE, 2(8), e671; 2007

Background Understanding gene expression changes over lifespan in diverse animal species will lead to insights to conserved processes in the biology of aging and allow development of interventions to improve health. Rotifers are small aquatic invertebrates that have been used in aging studies for nearly 100 years and are now re-emerging as a modern model system. To provide a baseline to evaluate genetic responses to interventions that change health throughout lifespan and a framework for new hypotheses about the molecular genetic mechanisms of aging, we examined the transcriptome of an asexual female lineage of the rotifer Brachionus manjavacas at five life stages: eggs, neonates, and early-, late-, and post-reproductive adults. Results There are widespread shifts in gene expression over the lifespan of B. manjavacas; the largest change occurs between neonates and early reproductive adults and is characterized by down-regulation of developmental genes and up-regulation of genes involved in reproduction. The expression profile of post-reproductive adults was distinct from that of other life stages. While few genes were significantly differentially expressed in the late- to post-reproductive transition, gene set enrichment analysis revealed multiple down-regulated pathways in metabolism, maintenance and repair, and proteostasis, united by genes involved in mitochondrial function and oxidative phosphorylation. Conclusions This study provides the first examination of changes in gene expression over lifespan in rotifers. We detected differential expression of many genes with human orthologs that are absent in Drosophila and C. elegans, highlighting the potential of the rotifer model in aging studies. Our findings suggest that small but coordinated changes in expression of many genes in pathways that integrate diverse functions drive the aging process. The observation of simultaneous declines in expression of genes in multiple pathways may have consequences for health and longevity not detected by single- or multi-gene knockdown in otherwise healthy animals. Investigation of subtle but genome-wide change in these pathways during aging is an important area for future study. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3540-x) contains supplementary material, which is available to authorized users.

Background Sexual reproduction is a widely studied biological process because it is critically important to the genetics, evolution, and ecology of eukaryotes. Despite decades of study on this topic, no comprehensive explanation has been accepted that explains the evolutionary forces underlying its prevalence and persistence in nature. Monogonont rotifers offer a useful system for experimental studies relating to the evolution of sexual reproduction due to their rapid reproductive rate and close relationship to the putatively ancient asexual bdelloid rotifers. However, little is known about the molecular underpinnings of sex in any rotifer species. Results We generated mRNA-seq libraries for obligate parthenogenetic (OP) and cyclical parthenogenetic (CP) strains of the monogonont rotifer, Brachionus calyciflorus, to identify genes specific to both modes of reproduction. Our differential expression analysis identified receptors with putative roles in signaling pathways responsible for the transition from asexual to sexual reproduction. Differential expression of a specific copy of the duplicated cell cycle regulatory gene CDC20 and specific copies of histone H2A suggest that such duplications may underlie the phenotypic plasticity required for reproductive mode switch in monogononts. We further identified differential expression of genes involved in the formation of resting eggs, a process linked exclusively to sex in this species. Finally, we identified transcripts from the bdelloid rotifer Adineta ricciae that have significant sequence similarity to genes with higher expression in CP strains of B. calyciflorus. Conclusions Our analysis of global gene expression differences between facultatively sexual and exclusively asexual populations of B. calyciflorus provides insights into the molecular nature of sexual reproduction in rotifers. Furthermore, our results offer insight into the evolution of obligate asexuality in bdelloid rotifers and provide indicators important for the use of monogononts as a model system for investigating the evolution of sexual reproduction.