# Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer Adineta vaga Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation

Boris Hespeels, Xiang Li, Jean-François Flot, Lise-Marie Pigneur, Jeremy Malaisse, Corinne da Silva, Karine van Doninck

► To cite this version:

Boris Hespeels, Xiang Li, Jean-François Flot, Lise-Marie Pigneur, Jeremy Malaisse, et al. Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer Adineta vaga Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation. PLoS ONE, 2015, 10 (7), pp.e0131313. 10.1371/journal.pone.0131313. cea-04573313

HAL Id: cea-04573313

https://cea.hal.science/cea-04573313

Submitted on 13 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer Adineta vaga Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation

Boris Hespeels, Xiang Li, Jean-François Flot, Lise-Marie Pigneur, Jeremy Malaisse, Corinne da Silva, Karine van Doninck

► To cite this version:

Boris Hespeels, Xiang Li, Jean-François Flot, Lise-Marie Pigneur, Jeremy Malaisse, et al. Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer Adineta vaga Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation. PLoS ONE, 2015, 10 (7), pp.e0131313. 10.1371/journal.pone.0131313. cea-04573313

HAL Id: cea-04573313

https://cea.hal.science/cea-04573313v1

Submitted on 13 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

CrossMark

click for updates

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Hespeels B, Li X, Flot J-F, Pigneur L-M, Malaisse J, Da Silva C, et al. (2015) Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer *Adineta vaga* Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation. PLoS ONE 10(7): e0131313. doi:10.1371/journal.pone.0131313

**Academic Editor:** Richard A Wilson, University of Nebraska-Lincoln, UNITED STATES

**Received:** February 19, 2015

**Accepted:** May 31, 2015

**Published:** July 10, 2015

**Copyright:** © 2015 Hespeels et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by UNamur through an assistant PhD grant allocated to BH and a starting grant allocated to KVD and by the Belgian National Fund for Scientific Research (Fonds National pour la Recherche Scientifique - FNRS) through a PhD FRIA (Fonds pour la formation à la Recherche dans l'Industrie et dans l'Agriculture) thesis grant allocated to XL and an FRFC (Fonds de la Recherche Fondamentale Collective) grant

RESEARCH ARTICLE

# Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer *Adineta vaga* Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation

Boris Hespeels¹, Xiang Li¹, Jean-François Flot², Lise-Marie Pigneur¹, Jeremy Malaisse³, Corinne Da Silva⁴, Karine Van Doninck¹,⁵*

1. LEGE laboratory, URBE, Department of Biology, University of Namur, Namur, Belgium, 2. Department of Genetics, Evolution and Environment, University College London, London, United Kingdom, 3. Research Unit for Molecular Physiology, Cell and Tissue Laboratory, Namur, Belgium, 4. CEA-Institut de Génomique, GENOSCOPE, Centre National de Séquençage, Evry, France, 5. naXys, Namur Center for Complex Systems, University of Namur, Namur, Belgium

* karine.vandoninck@unamur.be

## Abstract

The disaccharide sugar trehalose is essential for desiccation resistance in most metazoans that survive dryness; however, neither trehalose nor the enzymes involved in its metabolism have ever been detected in bdelloid rotifers despite their extreme resistance to desiccation. Here we screened the genome of the bdelloid rotifer *Adineta vaga* for genes involved in trehalose metabolism. We discovered a total of four putative trehalose-6-phosphate synthase (TPS) and seven putative trehalase (TRE) gene copies in the genome of this ameiotic organism; however, no trehalose-6-phosphate phosphatase (TPP) gene or domain was detected. The four TPS copies of *A. vaga* appear more closely related to plant and fungi proteins, as well as to some protists, whereas the seven TRE copies fall in bacterial clades. Therefore, *A. vaga* likely acquired its trehalose biosynthesis and hydrolysis genes by horizontal gene transfers. Nearly all residues important for substrate binding in the predicted TPS domains are highly conserved, supporting the hypothesis that several copies of the genes might be functional. Besides, RNAseq library screening showed that trehalase genes were highly expressed compared to TPS genes, explaining probably why trehalose had not been detected in previous studies of bdelloids. A strong overexpression of their TPS genes was observed when bdelloids enter desiccation, suggesting a possible signaling role of trehalose-6-phosphate or trehalose in this process.
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Hespeels B, Li X, Flot J-F, Pigneur L-M, Malaisse J, Da Silva C, et al. (2015) Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer *Adineta vaga* Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation. PLoS ONE 10(7): e0131313. doi:10.1371/journal.pone.0131313

**Academic Editor:** Richard A Wilson, University of Nebraska-Lincoln, UNITED STATES

**Received:** February 19, 2015

**Accepted:** May 31, 2015

**Published:** July 10, 2015

**Copyright:** © 2015 Hespeels et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by UNamur through an assistant PhD grant allocated to BH and a starting grant allocated to KVD and by the Belgian National Fund for Scientific Research (Fonds National pour la Recherche Scientifique - FNRS) through a PhD FRIA (Fonds pour la formation à la Recherche dans l'Industrie et dans l'Agriculture) thesis grant allocated to XL and an FRFC (Fonds de la Recherche Fondamentale Collective) grant

RESEARCH ARTICLE

# Against All Odds: Trehalose-6-Phosphate Synthase and Trehalase Genes in the Bdelloid Rotifer *Adineta vaga* Were Acquired by Horizontal Gene Transfer and Are Upregulated during Desiccation

Boris Hespeels¹, Xiang Li¹, Jean-François Flot², Lise-Marie Pigneur¹, Jeremy Malaisse³, Corinne Da Silva⁴, Karine Van Doninck¹,⁵*

1. LEGE laboratory, URBE, Department of Biology, University of Namur, Namur, Belgium, 2. Department of Genetics, Evolution and Environment, University College London, London, United Kingdom, 3. Research Unit for Molecular Physiology, Cell and Tissue Laboratory, Namur, Belgium, 4. CEA-Institut de Génomique, GENOSCOPE, Centre National de Séquençage, Evry, France, 5. naXys, Namur Center for Complex Systems, University of Namur, Namur, Belgium

* karine.vandoninck@unamur.be

## Abstract

The disaccharide sugar trehalose is essential for desiccation resistance in most metazoans that survive dryness; however, neither trehalose nor the enzymes involved in its metabolism have ever been detected in bdelloid rotifers despite their extreme resistance to desiccation. Here we screened the genome of the bdelloid rotifer *Adineta vaga* for genes involved in trehalose metabolism. We discovered a total of four putative trehalose-6-phosphate synthase (TPS) and seven putative trehalase (TRE) gene copies in the genome of this ameiotic organism; however, no trehalose-6-phosphate phosphatase (TPP) gene or domain was detected. The four TPS copies of *A. vaga* appear more closely related to plant and fungi proteins, as well as to some protists, whereas the seven TRE copies fall in bacterial clades. Therefore, *A. vaga* likely acquired its trehalose biosynthesis and hydrolysis genes by horizontal gene transfers. Nearly all residues important for substrate binding in the predicted TPS domains are highly conserved, supporting the hypothesis that several copies of the genes might be functional. Besides, RNAseq library screening showed that trehalase genes were highly expressed compared to TPS genes, explaining probably why trehalose had not been detected in previous studies of bdelloids. A strong overexpression of their TPS genes was observed when bdelloids enter desiccation, suggesting a possible signaling role of trehalose-6-phosphate or trehalose in this process.
2.4.655.09.F and a CDR research grant DAMAGE (19597258) allocated to KVD.

**Competing Interests:** The authors have declared that no competing interests exist.

Introduction

Rotifers are microscopic invertebrates characterized by a ciliated head structure and a jaw-like grinding organ, the mastax. Within the phylum Rotifera, bdelloids are of particular interest because of their survival and diversification for tens of million years without males, reproducing asexually [1,2]. Recently, our team revealed a particular genome structure in the bdelloid lineage *Adineta vaga* confirming this ameiotic evolution: the allelic regions were rearranged and sometimes found on the same chromosome preventing pairing of homologues chromosomes and proper segregation of alleles. Our genomic results and a recent transcriptomic study by Boschetti et al. [3] on the closely related bdelloid species *Adineta ricciae* also revealed that around 8% of the genes in the genome are of apparent non-metazoan origin and were therefore probably acquired horizontally. In *A. ricciae*, these genes of foreign origin appear to be expressed on a scale unprecedented in animals. Moreover, there is evidence for ancient HGTs in the genus *Adineta* when comparing both species [4] suggesting that HGTs have contributed significantly to the adaptation of bdelloid rotifers and may represent an advantage during their ameiotic evolution.

The 400 described morphospecies of bdelloid rotifers are ubiquitous inhabitants of nearly all freshwater habitats on Earth, with a preference for limno-terrestrial environments such as mosses, lichens and temporary freshwater pools [5]. Their ability to thrive in such habitats lies in their extreme resistance to desiccation at any stage of their life cycle. When habitats dry out, bdelloids enter a state of suspended life in the form of a “tun” (when their body becomes ovoid and contracted), then regain activity when water becomes available. This phenomenon, known as “anhydrobiosis” [6], was discovered in bdelloids over 300 years ago by Van Leeuwenhoek [7]. Anhydrobiosis is not common among metazoans and is often restricted to particular developmental stages such as the resting eggs of the monogonont rotifers [8] and *Daphnia* [9], the cysts of some crustaceans such as *Artemia* [10], and the larvae of a few insect species, including chironomids [11]. However, there are three well-known invertebrate phyla able to withstand desiccation at any stage in their life-cycle, including the adult form: tardigrades, nematodes and bdelloid rotifers [12]. These taxa represent unique model systems to investigate the desiccation phenomenon.

Desiccation survival involves complex changes and adaptations at morphological, physiological and molecular levels that either protect cellular components from desiccation-induced damage or promote their repair following rehydration. A number of candidate genes involved in desiccation resistance have been identified by genomic, transcriptomic and proteomic analyses in several species. There appears to be a set of molecular adaptations common to all “anhydroorganisms”, which has led Potts et al. [13] to coin the term “desiccome” to describe those “universal” molecules involved in desiccation tolerance. One of the critical components of the desiccome are non-reducing disaccharides found to accumulate during desiccation, such as trehalose (in bacteria, fungi, some resurrection plants and metazoans) and sucrose (in most plants, pollen and seeds) [14–16]. These sugars have been proposed to play a role in osmotic adjustment, to stabilize biomolecules and membranes, and to act as a replacement for water [17,18].

In Archaebacteria an accumulation of trehalose occurs in response to stress suggesting that the protective role of trehalose during cell dehydration might be an ancient adaptation that is evolutionary preserved [19,20]. Accumulation of high levels of trehalose was also observed in desiccated cysts of *Artemia fransiscana* and the desiccation tolerant chironomid larvae up to 15% and 20% respectively of dry body mass [10,21]. Accumulation of trehalose or overexpression of TPS genes upon desiccation was also reported in numerous desiccation-resistant nematodes [22,23]. In tardigrades, trehalose was found in all investigated species but the level was
2.4.655.09.F and a CDR research grant DAMAGE (19597258) allocated to KVD.

**Competing Interests:** The authors have declared that no competing interests exist.

Introduction

Rotifers are microscopic invertebrates characterized by a ciliated head structure and a jaw-like grinding organ, the mastax. Within the phylum Rotifera, bdelloids are of particular interest because of their survival and diversification for tens of million years without males, reproducing asexually [1,2]. Recently, our team revealed a particular genome structure in the bdelloid lineage *Adineta vaga* confirming this ameiotic evolution: the allelic regions were rearranged and sometimes found on the same chromosome preventing pairing of homologues chromosomes and proper segregation of alleles. Our genomic results and a recent transcriptomic study by Boschetti et al. [3] on the closely related bdelloid species *Adineta ricciae* also revealed that around 8% of the genes in the genome are of apparent non-metazoan origin and were therefore probably acquired horizontally. In *A. ricciae*, these genes of foreign origin appear to be expressed on a scale unprecedented in animals. Moreover, there is evidence for ancient HGTs in the genus *Adineta* when comparing both species [4] suggesting that HGTs have contributed significantly to the adaptation of bdelloid rotifers and may represent an advantage during their ameiotic evolution.

The 400 described morphospecies of bdelloid rotifers are ubiquitous inhabitants of nearly all freshwater habitats on Earth, with a preference for limno-terrestrial environments such as mosses, lichens and temporary freshwater pools [5]. Their ability to thrive in such habitats lies in their extreme resistance to desiccation at any stage of their life cycle. When habitats dry out, bdelloids enter a state of suspended life in the form of a “tun” (when their body becomes ovoid and contracted), then regain activity when water becomes available. This phenomenon, known as “anhydrobiosis” [6], was discovered in bdelloids over 300 years ago by Van Leeuwenhoek [7]. Anhydrobiosis is not common among metazoans and is often restricted to particular developmental stages such as the resting eggs of the monogonont rotifers [8] and *Daphnia* [9], the cysts of some crustaceans such as *Artemia* [10], and the larvae of a few insect species, including chironomids [11]. However, there are three well-known invertebrate phyla able to withstand desiccation at any stage in their life-cycle, including the adult form: tardigrades, nematodes and bdelloid rotifers [12]. These taxa represent unique model systems to investigate the desiccation phenomenon.

Desiccation survival involves complex changes and adaptations at morphological, physiological and molecular levels that either protect cellular components from desiccation-induced damage or promote their repair following rehydration. A number of candidate genes involved in desiccation resistance have been identified by genomic, transcriptomic and proteomic analyses in several species. There appears to be a set of molecular adaptations common to all “anhydroorganisms”, which has led Potts et al. [13] to coin the term “desiccome” to describe those “universal” molecules involved in desiccation tolerance. One of the critical components of the desiccome are non-reducing disaccharides found to accumulate during desiccation, such as trehalose (in bacteria, fungi, some resurrection plants and metazoans) and sucrose (in most plants, pollen and seeds) [14–16]. These sugars have been proposed to play a role in osmotic adjustment, to stabilize biomolecules and membranes, and to act as a replacement for water [17,18].

In Archaebacteria an accumulation of trehalose occurs in response to stress suggesting that the protective role of trehalose during cell dehydration might be an ancient adaptation that is evolutionary preserved [19,20]. Accumulation of high levels of trehalose was also observed in desiccated cysts of *Artemia fransiscana* and the desiccation tolerant chironomid larvae up to 15% and 20% respectively of dry body mass [10,21]. Accumulation of trehalose or overexpression of TPS genes upon desiccation was also reported in numerous desiccation-resistant nematodes [22,23]. In tardigrades, trehalose was found in all investigated species but the level was
close to detection limit in the species *Milnesium tardigradum* and there appears to be divergence in the responses to desiccation, suggesting that trehalose is not be the unique key molecule for desiccation resistance in these organisms [24,25].

Despite the presence of trehalose anabolytic and catabolytic enzymes in the whole plant kingdom, high amounts of trehalose (up to 12% of dry weight) were only detected in few desiccation-tolerant plants [20,26]. Surprisingly, attempts to improve drought tolerance by engineering genetically modified plants accumulating trehalose resulted in plants that were less tolerant to drought than wild-type strains [27]. There are therefore evidences that trehalose produced by plants does not primarily act as a bioprotectant. Instead, sucrose is hypothesized to take over the role of trehalose as a protective sugar during desiccation in plants [26,28]. Even in a highly desiccation-tolerant plant such as *Selaginella lepidophylla*, a recent publication suggested that the high level of trehalose detected is synthetized by plant endophytes and not by the plant itself [20]. The widespread occurrence and evolutionary conservation of genes involved in trehalose metabolism in the plant kingdom suggests their possible implication in something else than mere desiccation or drought resistance by trehalose accumulation. Indeed, in the last decade it has been highlighted repeatedly that the important roles of TPS and trehalose-6-phosphate are in regulating plant metabolism, growth, development and abiotic stress response [29–32].

Among the five biosynthetic pathways for trehalose recognized to date, the TPS/TPP pathway is the most common: it is found in the three domains of life and is the only pathway utilized by plants and metazoans (summarized in [33]). In this pathway, the synthesis of trehalose from glucose is catalyzed by the enzymes TPS and trehalose phosphatase (TPP) (Figure A in S1 File). The degradation of trehalose, on the other hand, most commonly involves the enzyme trehalase (TRE), which catalyzes the breakdown of trehalose into two D-glucose molecules (Figure A in S1 File) [33].

Surprisingly, disaccharides such as trehalose have never been detected in the bdelloid rotifers while most of the species appear desiccation resistant. Carbohydrate analysis by gas chromatography (GC) of samples extracted from desiccated *Philodina roseola*, *A. vaga* and *Macrotrachela quadricornifera* yielded no evidence of trehalose or any other disaccharide, and all attempts to amplify bdelloid TPS genes failed [34–36]. A recent analysis of a cDNA library of the bdelloid rotifer *Adineta ricciae* enriched for genes upregulated following dehydration for 24h at 98% relative humidity (RH) did not detect any genes involved in the biosynthesis of non-reducing disaccharides such as trehalose [37]. By contrast, in the resting eggs of the facultative sexual rotifers of the class Monogononta, the sister clade of bdelloids, trehalose was detected using GC [38]. Furthermore, EST sequencing revealed that sexual as well as amictic females of the monogonont species *Brachionus plicatilis* expressed TPS [8,38]. However, in contrast to the desiccated cysts of *Artemia fransiscana*, a low amount of trehalose (0.35% trehalose/dry weight) was detected in desiccated resting eggs of monogononts suggesting that trehalose may be involved in osmotic regulation in these organisms rather than in resistance to desiccation [36]. It is nevertheless surprising that the bdelloid rotifers lack trehalose metabolism while this sugar is present, even at low concentrations, in all desiccation resistant metazoans.

The recently sequenced genome of *A. vaga* comprises the most diverse repertoire of carbohydrate-active enzymes (CAZymes) reported among metazoans so far, including 623 glycoside hydrolases (GHs, involved in the hydrolysis of sugar bonds), and 412 glycosyltransferases (GTs, responsible for building sugar bonds) [3]. Inspired by this rich repertoire, the genome of *A. vaga* was screened here for candidate genes involved in trehalose biosynthesis and degradation and the phylogenetic origin of the TPS and TRE enzymes of bdelloid rotifers was verified. Moreover, cDNA libraries and quantitative real-time PCRs (qPCRs) were conducted at
close to detection limit in the species *Milnesium tardigradum* and there appears to be divergence in the responses to desiccation, suggesting that trehalose is not be the unique key molecule for desiccation resistance in these organisms [24,25].

Despite the presence of trehalose anabolytic and catabolytic enzymes in the whole plant kingdom, high amounts of trehalose (up to 12% of dry weight) were only detected in few desiccation-tolerant plants [20,26]. Surprisingly, attempts to improve drought tolerance by engineering genetically modified plants accumulating trehalose resulted in plants that were less tolerant to drought than wild-type strains [27]. There are therefore evidences that trehalose produced by plants does not primarily act as a bioprotectant. Instead, sucrose is hypothesized to take over the role of trehalose as a protective sugar during desiccation in plants [26,28]. Even in a highly desiccation-tolerant plant such as *Selaginella lepidophylla*, a recent publication suggested that the high level of trehalose detected is synthetized by plant endophytes and not by the plant itself [20]. The widespread occurrence and evolutionary conservation of genes involved in trehalose metabolism in the plant kingdom suggests their possible implication in something else than mere desiccation or drought resistance by trehalose accumulation. Indeed, in the last decade it has been highlighted repeatedly that the important roles of TPS and trehalose-6-phosphate are in regulating plant metabolism, growth, development and abiotic stress response [29–32].

Among the five biosynthetic pathways for trehalose recognized to date, the TPS/TPP pathway is the most common: it is found in the three domains of life and is the only pathway utilized by plants and metazoans (summarized in [33]). In this pathway, the synthesis of trehalose from glucose is catalyzed by the enzymes TPS and trehalose phosphatase (TPP) (Figure A in S1 File). The degradation of trehalose, on the other hand, most commonly involves the enzyme trehalase (TRE), which catalyzes the breakdown of trehalose into two D-glucose molecules (Figure A in S1 File) [33].

Surprisingly, disaccharides such as trehalose have never been detected in the bdelloid rotifers while most of the species appear desiccation resistant. Carbohydrate analysis by gas chromatography (GC) of samples extracted from desiccated *Philodina roseola*, *A. vaga* and *Macrotrachela quadricornifera* yielded no evidence of trehalose or any other disaccharide, and all attempts to amplify bdelloid TPS genes failed [34–36]. A recent analysis of a cDNA library of the bdelloid rotifer *Adineta ricciae* enriched for genes upregulated following dehydration for 24h at 98% relative humidity (RH) did not detect any genes involved in the biosynthesis of non-reducing disaccharides such as trehalose [37]. By contrast, in the resting eggs of the facultative sexual rotifers of the class Monogononta, the sister clade of bdelloids, trehalose was detected using GC [38]. Furthermore, EST sequencing revealed that sexual as well as amictic females of the monogonont species *Brachionus plicatilis* expressed TPS [8,38]. However, in contrast to the desiccated cysts of *Artemia fransiscana*, a low amount of trehalose (0.35% trehalose/dry weight) was detected in desiccated resting eggs of monogononts suggesting that trehalose may be involved in osmotic regulation in these organisms rather than in resistance to desiccation [36]. It is nevertheless surprising that the bdelloid rotifers lack trehalose metabolism while this sugar is present, even at low concentrations, in all desiccation resistant metazoans.

The recently sequenced genome of *A. vaga* comprises the most diverse repertoire of carbohydrate-active enzymes (CAZymes) reported among metazoans so far, including 623 glycoside hydrolases (GHs, involved in the hydrolysis of sugar bonds), and 412 glycosyltransferases (GTs, responsible for building sugar bonds) [3]. Inspired by this rich repertoire, the genome of *A. vaga* was screened here for candidate genes involved in trehalose biosynthesis and degradation and the phylogenetic origin of the TPS and TRE enzymes of bdelloid rotifers was verified. Moreover, cDNA libraries and quantitative real-time PCRs (qPCRs) were conducted at
different time points before, during and after desiccation of *A. vaga* to determine the expression profile of the detected genes involved in trehalose metabolism. While we did not study the activity of the detected trehalose enzymes, our results start to lift the veil on the origin of the trehalose genes of *A. vaga* and their metabolism.

## Materials and Methods

### Bdelloid rotifer cultures

Experiments were performed using isogenic *A. vaga* clones issued from a single individual from the laboratory of Matthew Meselson at Harvard University [39]. The cultures were maintained hydrated at 16°C in 150 × 20 mm Petri dishes supplemented with natural spring water (Spa) and fed with *Escherichia coli* MG1655.

### Annotation of candidate genes and sequence divergence calculations

Representative sequences of proteins involved in trehalose biosynthesis and hydrolysis pathways in eubacteria, archaea, plants, fungi and animals (Figure A in S1 File) were downloaded from the NCBI database and aligned against the genome of *A. vaga* using TBLASTN [40] as implemented in BioEdit [41] with a E-value threshold of 10<sup>-10</sup> (as recommended by the program MCScanX to detect homologues within and between genomes [42]). Intron-exon boundaries were identified by comparison with homologous amino acid sequences using GENEWISE2 [43] or, in a few cases, were identified manually based on the GT/AG rules for intron boundaries [44]. Translated CDS were used as queries for BLASTP [45] searches against the CAZyme repertoire found in the *A. vaga* genome [39] and the nr GenBank database. For each putative protein, the alien index (AI) [46] was determined as in [39]. Genes present on the genomic contigs of *A. vaga* were predicted according to the automated annotation of *A. vaga* [39] supplemented with manual curation and functional annotation as in Hur et al. [47].

The well-characterized sequence of *E. coli* TPS (OtsA), the 3D structure of which had been experimentally determined [48,49], was aligned with the TPS sequences found in *A. vaga* in order to check for conservation of amino acid residues involved in the binding of glucose 6-phosphate (acceptor) and UDP-glucose (donor). The alignment of the proteins was performed using MUSCLE [50] and manually edited with MEGA 6.06 [51] to remove non-conserved regions. Based on the previous protein alignment, a 3D-structure modeling of AvTpsA was performed using SWISS-MODEL (in Alignment Mode) [52–54] based on the known tridimensional structure of the OtsA chain B homotetramer in complex with G6P and UDP (1GZ5 [48]). A 3D alignment with the original 3D structure of OtsA was obtained using Pymol v1.3. In order to look for 3D conservation of active residues, an additional 3D structure of TPS1 was modeled using Phyre2 [55]. This second model was designed based on the 3D structure of OtsA chain A in complex with UDP-2-deoxy-2-fluoroglucose [49].

The sequence divergence between the identified candidate genes of *A. vaga* was calculated at the gene and protein sequence level using Geneious 5.4.6 (Biomatters).

### Sequencing *A. vaga* candidate *tps* and *tre* genes

Two overlapping sets of PCR primers for each *tps* gene copy and one set of primer for each *tre* gene copy identified in the genome of *A. vaga* were designed using Geneious 5.4.6 and Batch-Primer3 v1.0. Each 25 µL PCR contained 1X GoTaq reaction buffer (1.5 mM MgCl₂), 0.2 mM of each dNTP, 0.5 µM of each primer, 0.5 U of GoTaq DNA Polymerase (Promega), and 1 µL (ca. 10–200 ng) of genomic DNA. The amplification profile included an initial denaturation step at 94°C for 4 min, followed by 30 cycles of 45 s denaturation at 94°C, 45 s annealing at the
different time points before, during and after desiccation of *A. vaga* to determine the expression profile of the detected genes involved in trehalose metabolism. While we did not study the activity of the detected trehalose enzymes, our results start to lift the veil on the origin of the trehalose genes of *A. vaga* and their metabolism.

## Materials and Methods

### Bdelloid rotifer cultures

Experiments were performed using isogenic *A. vaga* clones issued from a single individual from the laboratory of Matthew Meselson at Harvard University [39]. The cultures were maintained hydrated at 16°C in 150 × 20 mm Petri dishes supplemented with natural spring water (Spa) and fed with *Escherichia coli* MG1655.

### Annotation of candidate genes and sequence divergence calculations

Representative sequences of proteins involved in trehalose biosynthesis and hydrolysis pathways in eubacteria, archaea, plants, fungi and animals (Figure A in S1 File) were downloaded from the NCBI database and aligned against the genome of *A. vaga* using TBLASTN [40] as implemented in BioEdit [41] with a E-value threshold of 10<sup>-10</sup> (as recommended by the program MCScanX to detect homologues within and between genomes [42]). Intron-exon boundaries were identified by comparison with homologous amino acid sequences using GENEWISE2 [43] or, in a few cases, were identified manually based on the GT/AG rules for intron boundaries [44]. Translated CDS were used as queries for BLASTP [45] searches against the CAZyme repertoire found in the *A. vaga* genome [39] and the nr GenBank database. For each putative protein, the alien index (AI) [46] was determined as in [39]. Genes present on the genomic contigs of *A. vaga* were predicted according to the automated annotation of *A. vaga* [39] supplemented with manual curation and functional annotation as in Hur et al. [47].

The well-characterized sequence of *E. coli* TPS (OtsA), the 3D structure of which had been experimentally determined [48,49], was aligned with the TPS sequences found in *A. vaga* in order to check for conservation of amino acid residues involved in the binding of glucose 6-phosphate (acceptor) and UDP-glucose (donor). The alignment of the proteins was performed using MUSCLE [50] and manually edited with MEGA 6.06 [51] to remove non-conserved regions. Based on the previous protein alignment, a 3D-structure modeling of AvTpsA was performed using SWISS-MODEL (in Alignment Mode) [52–54] based on the known tridimensional structure of the OtsA chain B homotetramer in complex with G6P and UDP (1GZ5 [48]). A 3D alignment with the original 3D structure of OtsA was obtained using Pymol v1.3. In order to look for 3D conservation of active residues, an additional 3D structure of TPS1 was modeled using Phyre2 [55]. This second model was designed based on the 3D structure of OtsA chain A in complex with UDP-2-deoxy-2-fluoroglucose [49].

The sequence divergence between the identified candidate genes of *A. vaga* was calculated at the gene and protein sequence level using Geneious 5.4.6 (Biomatters).

### Sequencing *A. vaga* candidate *tps* and *tre* genes

Two overlapping sets of PCR primers for each *tps* gene copy and one set of primer for each *tre* gene copy identified in the genome of *A. vaga* were designed using Geneious 5.4.6 and Batch-Primer3 v1.0. Each 25 µL PCR contained 1X GoTaq reaction buffer (1.5 mM MgCl₂), 0.2 mM of each dNTP, 0.5 µM of each primer, 0.5 U of GoTaq DNA Polymerase (Promega), and 1 µL (ca. 10–200 ng) of genomic DNA. The amplification profile included an initial denaturation step at 94°C for 4 min, followed by 30 cycles of 45 s denaturation at 94°C, 45 s annealing at the
appropriate temperature for each primer pair, and 50 s elongation at 72°C; and a final elonga-
tion step of 10 min at 72°C. The amplified fragments were Sanger-sequenced in both directions
at Genoscreen (Lille, France) to confirm the sequences of the identified *A. vaga* *tps* and *tre*
genes. Primers used for amplifying the trehalose-6-phosphate synthase (TPS) and trehalase
(TRE) genes are listed in Table A in S2 File.

Phylogenetic analysis

Since HGT is apparent in bdelloid rotifers, the phylogenetic trees were build based on an
exhaustive search of TPS and TRE genes found in the different domains of life. The selection of
all TPS genes given in Avonce et al. [33] were supplemented with metazoan sequences from
the Uniprot database. Moreover, sequences homologous to *A. vaga* genes were selected from
the best BLASTP hits against the non-redundant protein database of NCBI. For the trehalase
phylogenies, a sample of homologous TRE genes was selected from the Swiss-Prot database
using BLASTP, then enriched with the best BLASTP hits of *A. vaga* genes against the non-
redundant protein database of NCBI. Unpublished monogonont trehalase sequences were pro-
vided by Dr. David Mark Welch (MBL, Woods Hole).

The complete list of sequences used in our study is provided in Table B in S2 File. All TPS
and TRE amino acid sequences were aligned using MUSCLE [50] in MEGA 6.06 [51]. The
alignments were checked manually, then the best-fitting evolutionary model for each align-
ment was evaluated in PROTTEST 3.4 [56] using the Bayesian Information Criterion [57]
and the Akaike Information Criterion [58] (and concurred in selecting the WAG+I+G+F and
LG+I+G models for the TPS and TRE phylogenies, respectively).

Neighbor joining trees were constructed in MEGA 6.06 [51] using the JTT model and 1000
bootstraps. Maximum-likelihood trees were constructed in raxml GUI v1.3 [59,60] with 1000
rapid bootstraps according to the same evolutionary model as PROTTEST results. A Bayesian
phylogeny was inferred using MrBayes v3.1.1 [61] as implemented in Topali v2.5 [62]. Two
runs of 2.5 and 3 million generations were conducted for TPS and TRE respectively, sampling
trees every 10 generation and discarding the first 25% as burn-in. The PSRF (Potential scale
reduction factor) was in both case equal to 1, meaning that excellent convergence had been
achieved. We also used the probability-consensus pruning method (under maximum likeli-
hood) implemented in MetaPIGA v3.1 [63]. The program selected the best-fitting protein sub-
stitution model (WAG+G+F and GTR20+G for TPS and TRE respectively) implementing the
Bayesian Information Criterion. To estimate the posterior probability distribution of possible
trees, replicated metaGA searches were performed and stopped when a series of mean relative
error values among 10 consecutive consensus trees remained below 5% (between 100–10,000
replicates). Trees were edited using FigTree v1.4.0 (http://tree.bio.ed.ac.uk/software/figtree/).

Desiccation assays

Bdelloid rotifers of *A. vaga* were desiccated following the protocol published in [64]: briefly,
healthy *A. vaga* cultures were washed with filtered Spa water, individuals were collected using a
cell scraper, transferred to a 50-ml Falcon tube and centrifuged. Concentrated bdelloid individ-
uals were transferred with 1,2 ml supernatant to the center of Petri dishes containing 3% LMP
agarose (Invitrogen Carlsbad). Plates with similar amount of bdelloids were placed in a climatic
chamber (WEKK 0028) and submitted to the following desiccation protocol: first a linear
decrease of relative humidity from 70% to 55% for 17 h, then a linear decrease of relative
humidity from 55% to 41% for 1 h, and maintenance of relative humidity at 41% during the 14
days of desiccation. *A. vaga* individuals were rehydrated by adding 30 mL filtered Spa water to
each plate. During all the experiment, temperature was kept at 23°C. RNA for qPCR analysis
appropriate temperature for each primer pair, and 50 s elongation at 72°C; and a final elonga-
tion step of 10 min at 72°C. The amplified fragments were Sanger-sequenced in both directions
at Genoscreen (Lille, France) to confirm the sequences of the identified *A. vaga* *tps* and *tre*
genes. Primers used for amplifying the trehalose-6-phosphate synthase (TPS) and trehalase
(TRE) genes are listed in Table A in S2 File.

Phylogenetic analysis

Since HGT is apparent in bdelloid rotifers, the phylogenetic trees were build based on an
exhaustive search of TPS and TRE genes found in the different domains of life. The selection of
all TPS genes given in Avonce et al. [33] were supplemented with metazoan sequences from
the Uniprot database. Moreover, sequences homologous to *A. vaga* genes were selected from
the best BLASTP hits against the non-redundant protein database of NCBI. For the trehalase
phylogenies, a sample of homologous TRE genes was selected from the Swiss-Prot database
using BLASTP, then enriched with the best BLASTP hits of *A. vaga* genes against the non-
redundant protein database of NCBI. Unpublished monogonont trehalase sequences were pro-
vided by Dr. David Mark Welch (MBL, Woods Hole).

The complete list of sequences used in our study is provided in Table B in S2 File. All TPS
and TRE amino acid sequences were aligned using MUSCLE [50] in MEGA 6.06 [51]. The
alignments were checked manually, then the best-fitting evolutionary model for each align-
ment was evaluated in PROTTEST 3.4 [56] using the Bayesian Information Criterion [57]
and the Akaike Information Criterion [58] (and concurred in selecting the WAG+I+G+F and
LG+I+G models for the TPS and TRE phylogenies, respectively).

Neighbor joining trees were constructed in MEGA 6.06 [51] using the JTT model and 1000
bootstraps. Maximum-likelihood trees were constructed in raxml GUI v1.3 [59,60] with 1000
rapid bootstraps according to the same evolutionary model as PROTTEST results. A Bayesian
phylogeny was inferred using MrBayes v3.1.1 [61] as implemented in Topali v2.5 [62]. Two
runs of 2.5 and 3 million generations were conducted for TPS and TRE respectively, sampling
trees every 10 generation and discarding the first 25% as burn-in. The PSRF (Potential scale
reduction factor) was in both case equal to 1, meaning that excellent convergence had been
achieved. We also used the probability-consensus pruning method (under maximum likeli-
hood) implemented in MetaPIGA v3.1 [63]. The program selected the best-fitting protein sub-
stitution model (WAG+G+F and GTR20+G for TPS and TRE respectively) implementing the
Bayesian Information Criterion. To estimate the posterior probability distribution of possible
trees, replicated metaGA searches were performed and stopped when a series of mean relative
error values among 10 consecutive consensus trees remained below 5% (between 100–10,000
replicates). Trees were edited using FigTree v1.4.0 (http://tree.bio.ed.ac.uk/software/figtree/).

Desiccation assays

Bdelloid rotifers of *A. vaga* were desiccated following the protocol published in [64]: briefly,
healthy *A. vaga* cultures were washed with filtered Spa water, individuals were collected using a
cell scraper, transferred to a 50-ml Falcon tube and centrifuged. Concentrated bdelloid individ-
uals were transferred with 1,2 ml supernatant to the center of Petri dishes containing 3% LMP
agarose (Invitrogen Carlsbad). Plates with similar amount of bdelloids were placed in a climatic
chamber (WEKK 0028) and submitted to the following desiccation protocol: first a linear
decrease of relative humidity from 70% to 55% for 17 h, then a linear decrease of relative
humidity from 55% to 41% for 1 h, and maintenance of relative humidity at 41% during the 14
days of desiccation. *A. vaga* individuals were rehydrated by adding 30 mL filtered Spa water to
each plate. During all the experiment, temperature was kept at 23°C. RNA for qPCR analysis
and/or cDNA preparation was extracted at six different time points before, during or after desiccation (Figure B in S1 File).

## RNA extraction and cDNA library construction

Total RNAs were extracted using the RNAqueous-4PCR Kit (Ambion, Austin) at the six time points (Figure B in S1 File). Total RNA was enriched in mRNA based on its polyA tail, chemically fragmented and converted into single-stranded cDNA using random hexamer priming, generating double-stranded cDNA. Next, paired-end libraries were prepared following the Illumina 222s protocol (Illumina DNA sample kit): briefly, fragments were end-repaired, 3'-adenylated, and ligated with Illumina adapters. DNA fragments ranging in size from 300 to 600bp (including the adapters) were PCR-amplified using adapter-specific primers. Libraries were purified then quantified using a Qubit Fluorometer (Life technologies), and library profiles were evaluated using an Agilent 2100 bioanalyzer. A paired-end flow cell of 101-bp reads was sequenced for each library on an Illumina HiSeq2000 platform. Raw data were normalized per kilobase of gene and per billion of mapped reads (RPKB) to allow comparison between conditions. The cDNA reads were mapped against the putative open reading frames (ORFs) of the identified TPS and TRE genes using Gaze v2 [65] to quantify their expression at the different time points.

## Quantitative PCR

The expression levels of the TPS and TRE genes at each time point (Figure B in S1 File) were confirmed by conducting qPCRs. RNA extraction for the experiment followed the protocol detailed above. Each cDNA synthesis reaction was performed using 390 ng of total RNA, 200 U of M-MLV Reverse Transcriptase (Promega), 80 µM dNTP mix (Promega) and 15 ng/µL random hexamers (Invitrogen) in 20 µL of 1X M-MLV reaction buffer (Promega M531A) for 1 h at 42°C, then exposed for 5 min at 95°C and kept on ice. FastStart Universal SYBR Green Master (Rox) (Roche, Basel, Switzerland) was used for qPCR. The oligonucleotide primer sequences (300 nmol/L; IDT) for TPS and TRE genes are listed in Table C in S2 File. All amplification products were checked using BLAST searches against the published genome of *A. vaga*, gel visualization, dissociation curves and sequencing.

After investigating the stability of the expression of eight possible housekeeping genes (18S, 28S, 53A, L40, L32, α tubulin, β tubulin and actin) across 12 Illumina RNA-seq libraries obtained at different time points (see [39]) the L40 ribosomal subunit was chosen for the qPCR analyses because it had the most stable expression pattern (stdev/average = 0.54). Each time point was investigated in triplicate (to account for environmental variation) and each qPCR reaction was performed in duplicate (to account for experimental variation). Relative quantifications were performed following the ΔΔCt method [66] supplemented by analyses of variance (ANOVA1).

## Results

### Identification of the *Adineta vaga* genes involved in trehalose metabolism

**The trehalose biosynthesis pathway.** The TBLASTN search results for *A. vaga* homologues of trehalose phosphate synthase (TPS) from a wide range of species (bacteria, fungi, plants and metazoans) yielded significant E-values (less than $10^{-10}$) for four predicted genes belonging to three different scaffolds (av51, av255 and av681) in the *A. vaga* genomic dataset. Indeed, *A. vaga* was previously reported to be an ancient degenerate tetraploid [39], hence finding
and/or cDNA preparation was extracted at six different time points before, during or after desiccation (Figure B in S1 File).

## RNA extraction and cDNA library construction

Total RNAs were extracted using the RNAqueous-4PCR Kit (Ambion, Austin) at the six time points (Figure B in S1 File). Total RNA was enriched in mRNA based on its polyA tail, chemically fragmented and converted into single-stranded cDNA using random hexamer priming, generating double-stranded cDNA. Next, paired-end libraries were prepared following the Illumina 222s protocol (Illumina DNA sample kit): briefly, fragments were end-repaired, 3'-adenylated, and ligated with Illumina adapters. DNA fragments ranging in size from 300 to 600bp (including the adapters) were PCR-amplified using adapter-specific primers. Libraries were purified then quantified using a Qubit Fluorometer (Life technologies), and library profiles were evaluated using an Agilent 2100 bioanalyzer. A paired-end flow cell of 101-bp reads was sequenced for each library on an Illumina HiSeq2000 platform. Raw data were normalized per kilobase of gene and per billion of mapped reads (RPKB) to allow comparison between conditions. The cDNA reads were mapped against the putative open reading frames (ORFs) of the identified TPS and TRE genes using Gaze v2 [65] to quantify their expression at the different time points.

## Quantitative PCR

The expression levels of the TPS and TRE genes at each time point (Figure B in S1 File) were confirmed by conducting qPCRs. RNA extraction for the experiment followed the protocol detailed above. Each cDNA synthesis reaction was performed using 390 ng of total RNA, 200 U of M-MLV Reverse Transcriptase (Promega), 80 µM dNTP mix (Promega) and 15 ng/µL random hexamers (Invitrogen) in 20 µL of 1X M-MLV reaction buffer (Promega M531A) for 1 h at 42°C, then exposed for 5 min at 95°C and kept on ice. FastStart Universal SYBR Green Master (Rox) (Roche, Basel, Switzerland) was used for qPCR. The oligonucleotide primer sequences (300 nmol/L; IDT) for TPS and TRE genes are listed in Table C in S2 File. All amplification products were checked using BLAST searches against the published genome of *A. vaga*, gel visualization, dissociation curves and sequencing.

After investigating the stability of the expression of eight possible housekeeping genes (18S, 28S, 53A, L40, L32, α tubulin, β tubulin and actin) across 12 Illumina RNA-seq libraries obtained at different time points (see [39]) the L40 ribosomal subunit was chosen for the qPCR analyses because it had the most stable expression pattern (stdev/average = 0.54). Each time point was investigated in triplicate (to account for environmental variation) and each qPCR reaction was performed in duplicate (to account for experimental variation). Relative quantifications were performed following the ΔΔCt method [66] supplemented by analyses of variance (ANOVA1).

## Results

### Identification of the *Adineta vaga* genes involved in trehalose metabolism

**The trehalose biosynthesis pathway.** The TBLASTN search results for *A. vaga* homologues of trehalose phosphate synthase (TPS) from a wide range of species (bacteria, fungi, plants and metazoans) yielded significant E-values (less than $10^{-10}$) for four predicted genes belonging to three different scaffolds (av51, av255 and av681) in the *A. vaga* genomic dataset. Indeed, *A. vaga* was previously reported to be an ancient degenerate tetraploid [39], hence finding
more than two copies of a given gene was not surprising. Two TPS genes were found on a palindromic, highly colinear allelic region [39] of scaffold av51, and were annotated *AvTpsA* and *AvTpsA'* (Fig 1A). The surrounding region of the TPS copies on scaffolds av255 and av681 were also highly colinear and therefore probably allelic [39]: these two copies were annotated as *AvTpsB* and *AvTpsB'* (Fig 1A). One copy (*AvTpsB*) was characterized by the insertion of an A at position 223 (in comparison with *AvTpsB'* ) generating a stop codon in the translated protein. In *AvTpsB'* , the substitution of a C by a T at position 766 also resulted in a stop codon and a truncated protein sequence. We scrutinized the corresponding cDNA sequencing dataset and found no evidence for RNA editing (see below), suggesting that these two copies might be pseudogenes and we therefore prefixed their names with the Greek letter Ψ (Fig 1A). In downstream phylogenetic analyses, we restored both genes to their presumable full length by deleting the A in position 223 of Ψ*AvTpsB* and by replacing the T in position 766 of Ψ*AvTpsB'* by a C.

The putative allelic relationships between the copies, inferred from colinearity, were evalu-
ated by calculating the pairwise sequence divergence: *AvTpsA* and *AvTpsA'* were nearly identi-
cal and shared 98.2% identity at the nucleotide level and 98.3% at the protein level whereas
Ψ*AvTpsB* and Ψ*AvTpsB'* shared 96.9% identity at the nucleotide level and 94.5% at the protein
level. The two pairs (A and B) containing those genes were 71.3% identical at the nucleotide
level and 57.0% identical at the protein level. Although synteny was well conserved within each
pair of TPS genes, no synteny was detected between the two pairs except at the *Tps* position
(Fig 1A), suggesting that the four contigs containing the TPS genes of *A. vaga* do not form a
quartet of ohnologues such as reported previously in *A. vaga* for hsp82 or histone H3 [47,67].
All the four *Tps* gene copies had been previously annotated in the repertoire of CAZymes of
the *A. vaga* genome [39].

Within the *A. vaga* genomic dataset no significant hit to the trehalose phosphate phospha-
tase (TPP) gene present in bacteria, fungi, plants and metazoans (Figure A in S1 File) was
detected. The TPP catalytic site of the gene is characterized by three short, highly conserved
motifs found in archaea, eubacteria, plants and opisthokonts [33,68], and it has been suggested
that TPP domains became fused to the TPS protein in some bacteria and archaea as well as in
most eukaryotes [33,69,70]. A close examination of the C-terminal regions of the four pre-
dicted TPS proteins of *A. vaga* did however not reveal any putative TPP domain.

In order to check whether the amino acid residues important for UDP and G6P binding
were conserved in the *A. vaga* TPS proteins, their sequences were aligned with the OtsA
sequence of *Escherichia coli* and with the TPS sequences of *Arabidopsis thaliana* and *Plasmo-
diophora brassicae* (respectively *AtTPS1* and *PbTPS1* [71]) that are within a same clade as *A.
vaga* in our TPS phylogenies (see below). Both *AtTPS1* and *PbTPS* are known to be fully func-
tional. Alignment of the four predicted TPS proteins of *A. vaga* with the OtsA protein sequence
of *E. coli* revealed that 16 out of 19 residues important for UDP and G6P binding were con-
served in *AvTpsA* and *AvTpsA'*, whereas 15 and 13 out of 19 of these important residues were
conserved in *ΨAvTpsB* and *ΨAvTpsB'* respectively (Figure C in S1 File). By comparison, 16
out of 19 residues are conserved between *A. thaliana* and *E. coli* [69] and 17 out of 19 are con-
served between *P. brassicae* and *E. coli*. Interestingly, 18 out of 19 residues were conserved
between *A. vaga AvTpsA-A'* and *A. thaliana* TPS1 and 17 out of 19 with *P. brassicae* TPS1
(Figure C in S1 File). The predicted 3D models of *AvTpsA* were close to published structure of
the OtsA of *E. coli*, with a root-mean-square equal to 0.166 Å and 0.234 Å for the 3D structures
predicted by SwissModel and Phyre2 respectively (Fig 2A and 2B). Comparing the residues
involved in UDP (Fig 2B) and G6P binding (Fig 2C) between the crystal structure of OtsA and
the predicted protein structure of *AvTpsA*, *AtTPS1* and *PbTPS1* also suggested that the
more than two copies of a given gene was not surprising. Two TPS genes were found on a palindromic, highly colinear allelic region [39] of scaffold av51, and were annotated *AvTpsA* and *AvTpsA'* (Fig 1A). The surrounding region of the TPS copies on scaffolds av255 and av681 were also highly colinear and therefore probably allelic [39]: these two copies were annotated as *AvTpsB* and *AvTpsB'* (Fig 1A). One copy (*AvTpsB*) was characterized by the insertion of an A at position 223 (in comparison with *AvTpsB'* ) generating a stop codon in the translated protein. In *AvTpsB'* , the substitution of a C by a T at position 766 also resulted in a stop codon and a truncated protein sequence. We scrutinized the corresponding cDNA sequencing dataset and found no evidence for RNA editing (see below), suggesting that these two copies might be pseudogenes and we therefore prefixed their names with the Greek letter Ψ (Fig 1A). In downstream phylogenetic analyses, we restored both genes to their presumable full length by deleting the A in position 223 of Ψ*AvTpsB* and by replacing the T in position 766 of Ψ*AvTpsB'* by a C.

The putative allelic relationships between the copies, inferred from colinearity, were evalu-
ated by calculating the pairwise sequence divergence: *AvTpsA* and *AvTpsA'* were nearly identi-
cal and shared 98.2% identity at the nucleotide level and 98.3% at the protein level whereas
Ψ*AvTpsB* and Ψ*AvTpsB'* shared 96.9% identity at the nucleotide level and 94.5% at the protein
level. The two pairs (A and B) containing those genes were 71.3% identical at the nucleotide
level and 57.0% identical at the protein level. Although synteny was well conserved within each
pair of TPS genes, no synteny was detected between the two pairs except at the *Tps* position
(Fig 1A), suggesting that the four contigs containing the TPS genes of *A. vaga* do not form a
quartet of ohnologues such as reported previously in *A. vaga* for hsp82 or histone H3 [47,67].
All the four *Tps* gene copies had been previously annotated in the repertoire of CAZymes of
the *A. vaga* genome [39].

Within the *A. vaga* genomic dataset no significant hit to the trehalose phosphate phospha-
tase (TPP) gene present in bacteria, fungi, plants and metazoans (Figure A in S1 File) was
detected. The TPP catalytic site of the gene is characterized by three short, highly conserved
motifs found in archaea, eubacteria, plants and opisthokonts [33,68], and it has been suggested
that TPP domains became fused to the TPS protein in some bacteria and archaea as well as in
most eukaryotes [33,69,70]. A close examination of the C-terminal regions of the four pre-
dicted TPS proteins of *A. vaga* did however not reveal any putative TPP domain.

In order to check whether the amino acid residues important for UDP and G6P binding
were conserved in the *A. vaga* TPS proteins, their sequences were aligned with the OtsA
sequence of *Escherichia coli* and with the TPS sequences of *Arabidopsis thaliana* and *Plasmo-
diophora brassicae* (respectively *AtTPS1* and *PbTPS1* [71]) that are within a same clade as *A.
vaga* in our TPS phylogenies (see below). Both *AtTPS1* and *PbTPS* are known to be fully func-
tional. Alignment of the four predicted TPS proteins of *A. vaga* with the OtsA protein sequence
of *E. coli* revealed that 16 out of 19 residues important for UDP and G6P binding were con-
served in *AvTpsA* and *AvTpsA'*, whereas 15 and 13 out of 19 of these important residues were
conserved in *ΨAvTpsB* and *ΨAvTpsB'* respectively (Figure C in S1 File). By comparison, 16
out of 19 residues are conserved between *A. thaliana* and *E. coli* [69] and 17 out of 19 are con-
served between *P. brassicae* and *E. coli*. Interestingly, 18 out of 19 residues were conserved
between *A. vaga AvTpsA-A'* and *A. thaliana* TPS1 and 17 out of 19 with *P. brassicae* TPS1
(Figure C in S1 File). The predicted 3D models of *AvTpsA* were close to published structure of
the OtsA of *E. coli*, with a root-mean-square equal to 0.166 Å and 0.234 Å for the 3D structures
predicted by SwissModel and Phyre2 respectively (Fig 2A and 2B). Comparing the residues
involved in UDP (Fig 2B) and G6P binding (Fig 2C) between the crystal structure of OtsA and
the predicted protein structure of *AvTpsA*, *AtTPS1* and *PbTPS1* also suggested that the
a.

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<table><tr><td>1, 7, 10, 12, 13, 24</td><td>Hypothetical proteins</td><td>15</td><td>Phage-related minor tail protein (COG5218)</td></tr><tr><td>2</td><td>Metallophosphatase (MPP) superfamily</td><td>15', 16</td><td>Topoisomerase II-associated protein (PAT1)</td></tr><tr><td>3, 20</td><td>Trehalose-6-phosphate synthase (TPS)</td><td>17</td><td>Cadherin repeat-like superfamily</td></tr><tr><td>4, 4'</td><td>BNR (bacterial neuraminidase repeat) repeat-like domain</td><td>18</td><td>Iron-binding zinc finger CDGSH-type superfamily</td></tr><tr><td>5, 5'</td><td>FAD-dependent oxidoreductase superfamily</td><td>19</td><td>Aspartate aminotransferase (AAT) superfamily (fold type I)</td></tr><tr><td>6</td><td>Limulus factor C, Coch-5b2, and Lg11 (LCCL) superfamily</td><td>21, 25</td><td>Major facilitator superfamily (MFS)</td></tr><tr><td>8</td><td>Solute carrier families 5 and 6-like (SLC5-6-like sbd) superfamily</td><td>22</td><td>JmC domain(Cupin) superfamily</td></tr><tr><td>9</td><td>C-terminal, alpha helical domain of the glutathione S-transferase family (GST_C) superfamily</td><td>23</td><td>RNA recognition motif (RRM) superfamily</td></tr><tr><td>9'</td><td>Elongation factor 1 beta (EF1B) superfamily</td><td>26</td><td>Tetratricopeptide repeat (TPR)</td></tr><tr><td>11</td><td>Latrophilin/CL-1-like GPS domain</td><td>27</td><td>Ubiquitin-like domain containing C-terminal domain (CTD) phosphatase 1</td></tr><tr><td>11'</td><td>Lipase/floxygenase domain (PLAT/LH2 domain) superfamily</td><td>27'</td><td>Halocid dehalogenase (HAD)-like superfamily</td></tr><tr><td>14</td><td>RING- finger (Really Interesting New Gene) superfamily</td><td>28</td><td>Cytochrome b5 (Cyt-b5-like heme/steroid binding superfamily</td></tr><tr><td>15</td><td>Protein-interacting V-domain of mammalian Alix and related domains (V_Alix_like) superfamily</td><td>28'</td><td>Membrane fatty acid desaturase (Membrane_FADS)-like superfamily</td></tr></table>

b.

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<table><tr><td>1</td><td>Tubby C 2 (Tub-2) superfamily</td><td>18, 30</td><td>NADB Rossmann superfamily</td></tr><tr><td>2, 5</td><td>RNA recognition motif (RRM) superfamily</td><td>20, 28</td><td>Transient-receptor-potential calcium channel superfamily</td></tr><tr><td>3</td><td>Retinal pigment epithelial membrane protein (RPE65) superfamily</td><td>22</td><td>Serine Proteinase Inhibitors (serpins) superfamily</td></tr><tr><td>4, 6', 7, 8, 10, 21, 27, 31, 32, 35, 37, 42, 43, 44, 45, 48</td><td>Hypothetical proteins</td><td>23</td><td>Ubiquitin-like protein</td></tr><tr><td>6</td><td>cGMP-specific phosphodiesterases, adenyl cyclases and Fh1A (GAF) domain</td><td>24, 25, 33</td><td>NHL repeat containing-like protein</td></tr><tr><td>6'</td><td>H2D motif (HDc) domain</td><td>26</td><td>BED zinc finger (zfBed) superfamily</td></tr><tr><td>9, 19</td><td>Trehalase</td><td>29</td><td>Tollil-1 Receptor (TIR) superfamily</td></tr><tr><td>11</td><td>DnaJ/Hsp40 (heat shock protein 40) protein</td><td>34</td><td>Oligopeptide transporter (OPT) superfamily</td></tr><tr><td>12</td><td>RING- finger (Really Interesting New Gene) superfamily</td><td>36</td><td>ADP-ribose/glycyrolase (ADP-ribose/GH) superfamily</td></tr><tr><td>13</td><td>Josephin domain</td><td>38, 39</td><td>Low-density lipoprotein (LDL) receptor superfamily</td></tr><tr><td>14</td><td>Vault protein inter-alpha-trypsin (VIT) domain</td><td>40</td><td>Fatty acid hydroxylase superfamily</td></tr><tr><td>14'</td><td>Von Willebrand factor type A (vWA) domain</td><td>41</td><td>Calcium-binding epidermal growth factor (EGF-C)-like superfamily</td></tr><tr><td>15</td><td>Tetratricopeptide repeat (TPR)</td><td>47</td><td>Hedgehog/Integrin (Hint) superfamily</td></tr><tr><td>16</td><td>Integrase (rve) superfamily</td><td>48</td><td>Mitochondrial carrier superfamily</td></tr><tr><td>17</td><td>Trypsin-like serine protease (Tryp_Spc) superfamily</td><td></td><td></td></tr></table>
a.

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<table><thead><tr><th></th><th></th><th></th></tr></thead><tbody><tr><td>1, 7, 10, 12, 13, 24</td><td>Hypothetical proteins</td><td>15<sup>1</sup></td></tr><tr><td>2</td><td>Metallophosphatase (MPP) superfamily</td><td>15<sup>1</sup>, 16</td></tr><tr><td>3, 20</td><td>Trehalose-6-phosphate synthase (TPS)</td><td>17</td></tr><tr><td>4, 4'</td><td>BNR (bacterial neuraminidase repeat) repeat-like domain</td><td>18</td></tr><tr><td>5, 5'</td><td>FAD-dependent oxidoreductase superfamily</td><td>19</td></tr><tr><td>6</td><td>Limulus factor C, Coch-562, and Lg11 (LCCL) superfamily</td><td>21, 25</td></tr><tr><td>8</td><td>Solute carrier families 5 and 6-like (SLC5-6-like sbd) superfamily</td><td>22</td></tr><tr><td>9</td><td>C-terminal, alpha helical domain of the glutathione S-transferase family (GST_C) superfamily</td><td>23</td></tr><tr><td>9'</td><td>Elongation factor 1 beta (EF1B) superfamily</td><td>26</td></tr><tr><td>11</td><td>Latrophilin/CL-1-like GPS domain</td><td>27</td></tr><tr><td>11'</td><td>Lipase/floxygenase domain (PLAT/LH2 domain) superfamily</td><td>27'</td></tr><tr><td>14</td><td>RING- finger (Really Interesting New Gene) superfamily</td><td>28</td></tr><tr><td>15</td><td>Protein-interacting V-domain of mammalian Alix and related domains (V_Alix_like) superfamily</td><td>28<sup>1</sup></td></tr></tbody></table>

b.

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<table><thead><tr><th></th><th></th><th></th></tr></thead><tbody><tr><td>1</td><td>Tubby C 2 (Tub-2) superfamily</td><td>18, 30</td></tr><tr><td>2, 5</td><td>RNA recognition motif (RRM) superfamily</td><td>20, 28</td></tr><tr><td>3</td><td>Retinal pigment epithelial membrane protein (RPE65) superfamily</td><td>22</td></tr><tr><td>4, 6', 7, 8, 10, 21, 27, 31, 32, 35, 37, 42, 43, 44, 45, 48</td><td>Hypothetical proteins</td><td>23</td></tr><tr><td>6</td><td>cGMP-specific phosphodiesterases, adenyl cyclases and Fh1A (GAF) domain</td><td>24, 25, 33</td></tr><tr><td>6'</td><td>H3 histone (HDc) domain</td><td>26</td></tr><tr><td>9, 19</td><td>Trehalase</td><td>29</td></tr><tr><td>11</td><td>DnaJ/Hsp40 (heat shock protein 40) protein</td><td>34</td></tr><tr><td>12</td><td>RING- finger (Really Interesting New Gene) superfamily</td><td>36</td></tr><tr><td>13</td><td>Josephin domain</td><td>38, 39</td></tr><tr><td>14</td><td>Vault protein inter-alpha-trypsin (VIT) domain</td><td>40</td></tr><tr><td>14'</td><td>Von Willebrand factor type A (vWA) domain</td><td>41</td></tr><tr><td>15</td><td>Tetratricopeptide repeat (TPR)</td><td>47</td></tr><tr><td>16</td><td>Integrase (rve) superfamily</td><td>48</td></tr><tr><td>17</td><td>Trypsin-like serine protease (Tryp_Spc) superfamily</td><td></td></tr></tbody></table>
**Fig 1. Annotation of the scaffolds containing (a) TPS and (b) TRE genes in *Adineta vaga*.** Color-filled rectangles represent genes, numbers refer to the list at the bottom of the figure. Genes were colored in green or red if AI<45 or AI≥45, respectively. The genes annotated automatically in the *A. vaga* genome are surrounded with continuous lines, whereas discontinued lines surround additional genes detected during manual curation. The gene orientation is indicated by the position above (forward direction) or below (reverse direction) the line.

doi:10.1371/journal.pone.0131313.g001

topology of these residues was conserved (despite few amino acid substitutions between the
tested species).

No homologues of the genes involved in the TreP, TreT and TreY pathways of trehalose
biosynthesis (Figure A in S1 File) were detected in the genome of *A. vaga*.

**The trehalose hydrolysis pathway.** The TBLASTN alignments of the TRE proteins from a wide range of species (bacteria, fungi, plants and metazoa) against the *A. vaga* genome yielded six significant hits (E-values < 10<sup>-10</sup>) to five different scaffolds (av458, av6, av367, av920, av904) (Fig 1B and Table 1). Two TRE genes were found on scaffold av458 in a palindromic, highly colinear allelic region [3], and were therefore annotated *AvTreA* and *AvTreA'*. The regions surrounding two other TRE copies on scaffold av367 and av920 were highly colinear and therefore likely allelic [39]: these two allelic copies were annotated as *AvTreB* and *AvTreB'*. The TRE gene found on scaffold av6 and the adjacent genes exhibited a coverage twice the genome average suggesting that two allelic TRE copies had been fused: we labeled these two copies *AvTreC* and *AvTreC'*. An additional copy that had not been annotated in the published *A. vaga* genome was found in the terminal region of scaffold av904 and labeled *AvTreD*. The genomic region homologous to scaffold av904, scaffold av275, did not contain any detectable trehalase gene copy (Fig 1B), making a total of seven TRE gene copies in the published *A. vaga* genome. Since *AvTreD* had a stop codon located at the 886–888 bp of its coding region (and its cDNA sequence did not reveal evidence for RNA editing), we considered it to be a pseudogene and therefore prefixed its name with the Greek letter Ψ. The putative allelic relations between TRE copies were confirmed by analyzing the pairwise divergences between their sequences (see Table D in S2 File). In the regions surrounding the TRE genes on scaffolds av458 and av6, synteny was well preserved (Fig 1B), suggesting that *AvTreA-A'* and *AvTreC-C'* were ohnologues and form a quartet of four homologous regions (as in Hur et al and Van Doninck et al). The three other copies did not present such organization. The two other existing classes of trehalose hydrolytic enzymes (phosphotrehalase and trehalose phosphorylase), which had never been observed in animals so far, were not detected in the *A. vaga* genome either.

Alien index, introns and GC content

All the TPS and TRE genes identified had an alien index AI>45 (Table 1), an indication of their probable non-metazoan origin. At most one intron was detected in these genes (Table 1), which is significantly fewer than reported for probable core metazoan genes of *A. vaga* (AI≤-45, mean intron number 7.9 [39]). Seven out of the 11 TPS and TRE genes had a GC content (%) more than one percent below the genome average (33.3% [39]) (Table 1), another signature of plausible horizontal origin (see Flot et al. 2013).

Phylogenetic placement

**Phylogenetic tree of TPS gene sequences.** Tree topologies obtained by RaxML, Mr Bayes, Neighbor-Joining and MetaPIGA confirmed that the TPS sequences were grouped into two major branches as previously described by Avonce et al. [33] (Fig 3): the first group contains sequences from fungi, plants and some protists, whereas the second group consists of
**Fig 1. Annotation of the scaffolds containing (a) TPS and (b) TRE genes in *Adineta vaga*.** Color-filled rectangles represent genes, numbers refer to the list at the bottom of the figure. Genes were colored in green or red if AI<45 or AI≥45, respectively. The genes annotated automatically in the *A. vaga* genome are surrounded with continuous lines, whereas discontinued lines surround additional genes detected during manual curation. The gene orientation is indicated by the position above (forward direction) or below (reverse direction) the line.

doi:10.1371/journal.pone.0131313.g001

topology of these residues was conserved (despite few amino acid substitutions between the
tested species).

No homologues of the genes involved in the TreP, TreT and TreY pathways of trehalose
biosynthesis (Figure A in S1 File) were detected in the genome of *A. vaga*.

**The trehalose hydrolysis pathway.** The TBLASTN alignments of the TRE proteins from a wide range of species (bacteria, fungi, plants and metazoa) against the *A. vaga* genome yielded six significant hits (E-values < 10<sup>-10</sup>) to five different scaffolds (av458, av6, av367, av920, av904) (Fig 1B and Table 1). Two TRE genes were found on scaffold av458 in a palindromic, highly colinear allelic region [3], and were therefore annotated *AvTreA* and *AvTreA'*. The regions surrounding two other TRE copies on scaffold av367 and av920 were highly colinear and therefore likely allelic [39]: these two allelic copies were annotated as *AvTreB* and *AvTreB'*. The TRE gene found on scaffold av6 and the adjacent genes exhibited a coverage twice the genome average suggesting that two allelic TRE copies had been fused: we labeled these two copies *AvTreC* and *AvTreC'*. An additional copy that had not been annotated in the published *A. vaga* genome was found in the terminal region of scaffold av904 and labeled *AvTreD*. The genomic region homologous to scaffold av904, scaffold av275, did not contain any detectable trehalase gene copy (Fig 1B), making a total of seven TRE gene copies in the published *A. vaga* genome. Since *AvTreD* had a stop codon located at the 886–888 bp of its coding region (and its cDNA sequence did not reveal evidence for RNA editing), we considered it to be a pseudogene and therefore prefixed its name with the Greek letter Ψ. The putative allelic relations between TRE copies were confirmed by analyzing the pairwise divergences between their sequences (see Table D in S2 File). In the regions surrounding the TRE genes on scaffolds av458 and av6, synteny was well preserved (Fig 1B), suggesting that *AvTreA-A'* and *AvTreC-C'* were ohnologues and form a quartet of four homologous regions (as in Hur et al and Van Doninck et al). The three other copies did not present such organization. The two other existing classes of trehalose hydrolytic enzymes (phosphotrehalase and trehalose phosphorylase), which had never been observed in animals so far, were not detected in the *A. vaga* genome either.

Alien index, introns and GC content

All the TPS and TRE genes identified had an alien index AI>45 (Table 1), an indication of their probable non-metazoan origin. At most one intron was detected in these genes (Table 1), which is significantly fewer than reported for probable core metazoan genes of *A. vaga* (AI≤-45, mean intron number 7.9 [39]). Seven out of the 11 TPS and TRE genes had a GC content (%) more than one percent below the genome average (33.3% [39]) (Table 1), another signature of plausible horizontal origin (see Flot et al. 2013).

Phylogenetic placement

**Phylogenetic tree of TPS gene sequences.** Tree topologies obtained by RaxML, Mr Bayes, Neighbor-Joining and MetaPIGA confirmed that the TPS sequences were grouped into two major branches as previously described by Avonce et al. [33] (Fig 3): the first group contains sequences from fungi, plants and some protists, whereas the second group consists of
<figure><img src="image_10.png" /><figcaption>b.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>c.</figcaption></figure>

<figure><img src="image_12.png" /></figure>
<figure><img src="image_10.png" /><figcaption>b.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>c.</figcaption></figure>

<figure><img src="image_12.png" /></figure>
**Fig 2.** (a.) 3D alignment of the OtsA chain B crystal structure obtained for *E. coli* (white) with the two protein structures of TPS1 modeled using SWISS-MODEL (light blue) and PHYRE2 (orange). UDP (green) and G6P (red) are localized as observed in the OtsA crystal structure. (b.) Close-up on the conserved sites involved in the binding of G6P (red). The residues of the *E. coli* OtsA are shown in white and those of the AvTpsA predicted using SWISS-MODEL in light blue. SwissModel TPS1 *A. thaliana* (light green) and SwissModel predicted TPS1 *P. brassicacea* (gold) are also displayed in order to compare their topologies. (c.) Close-up on the conserved sites related to the binding of UDP (green). The residues are labeled based on the 3D structure of AvTpsA predicted using SWISS-MODEL.

doi:10.1371/journal.pone.0131313.g002

metazoans and bacteria. The TPS sequences of *A. vaga* were part of the clade of plant and functional TPS proteins (strongly supported by bootstrap branch supports of 100-97-84 for RaxML, NJ, and MetaPIGA respectively, and by a posterior probability of 1 using Mr Bayes). The two *A. vaga* AvTpsA and AvTpsA' sequences clustered with the protist TPS sequence of *P. brassicae* (bootstrap supports of 95-76-62, posterior probability of 1) and formed a sister clade to the plant Class I TPS proteins (bootstrap supports of 100-97-87, posterior probability of 1) confirming their non-metazoan origin. Most of these sequences have biochemically proven TPS enzymatic activity, except AtTPS2-AtTPS4 [71]. The position of ΨAvTpsB and ΨAvTpsB' within the tree was not resolved as the phylogenetic trees were not congruent or branch support was insufficient. In contrast, the TPS sequences of monogonont rotifers (Brachionus manjavacas and Brachionus calyciflorus) grouped outside the plant and fungal clade and were associated with other metazoans (bootstrap support values 100-100-99, posterior probability 1). As previously described in [33], the nematode sequences did not cluster with those of other metazoans.

<table><caption>Table 1. Summary of top BLASTP hits for the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes of *Adineta vaga* against GenBank.</caption><thead><tr><th>Gene ID, name</th><th>Scaffold</th><th>Nb of introns</th><th>GC (%)</th><th>Al</th><th>Best hit E-value</th><th>Metazoan best hit E-value</th><th>Best hit, taxonomy</th><th>Best hit, description</th></tr></thead><tbody><tr><td>GSADVT00013141001, AvTpsA</td><td>Av51</td><td>1</td><td>33.2</td><td>165</td><td>1.00E-161</td><td>4.00E-90</td><td>Eukaryota; Protozoa; Rhizaria; Cercozoa</td><td>glycosyltransferase family 20 protein</td></tr><tr><td>GSADVT00013173001, AvTpsA'</td><td>Av51</td><td>1</td><td>33.5</td><td>174</td><td>4.00E-162</td><td>2.00E-86</td><td>Eukaryota; Protozoa; Rhizaria; Cercozoa</td><td>glycosyltransferase family 20 protein</td></tr><tr><td>GSADVT00043720001, ΨAvTpsB*</td><td>Av255</td><td>0</td><td>31.1</td><td>82</td><td>1.00E-107</td><td>9.00E-72</td><td>Eukaryota; Fungi; Chytridiomycota</td><td>glycosyltransferase family 20 protein</td></tr><tr><td>GSADVT00063692001, ΨAvTpsB'*</td><td>Av681</td><td>0</td><td>31.2</td><td>78</td><td>1.00E-102</td><td>8.00E-69</td><td>Eukaryota; Fungi; Dikarya</td><td>glycosyltransferase family 20 protein</td></tr><tr><td>GSADVT00056365001, AvTreA</td><td>Av458</td><td>0</td><td>31.7</td><td>140</td><td>1.00E-147</td><td>7.00E-87</td><td>Bacteria; Bacteroidetes; Sphingobacteria</td><td>alpha, alpha-trehalase</td></tr><tr><td>GSADVT00056372001, AvTreA'</td><td>Av458</td><td>0</td><td>31.7</td><td>140</td><td>1.00E-147</td><td>1.00E-86</td><td>Bacteria; Bacteroidetes; Sphingobacteria</td><td>alpha, alpha-trehalase</td></tr><tr><td>GSADVT00051563001, AvTreB</td><td>Av367</td><td>1</td><td>27</td><td>281</td><td>0.00E+00</td><td>1.00E-78</td><td>Bacteria; Bacteroidetes; Cytophagia</td><td>alpha, alpha-trehalase</td></tr><tr><td>GSADVT00021696001, AvTreB'</td><td>Av920</td><td>1</td><td>26.9</td><td>233</td><td>6.00E-180</td><td>1.00E-78</td><td>Bacteria; Bacteroidetes; Cytophagia</td><td>alpha, alpha-trehalase</td></tr><tr><td>GSADVT00001947001, AvTreC-C'</td><td>Av6</td><td>0</td><td>33.6</td><td>139</td><td>2.00E-149</td><td>4.00E-89</td><td>Bacteria; Bacteroidetes; Cytophagia</td><td>alpha, alpha-trehalase</td></tr><tr><td>N/A, ΨAvTreD*</td><td>Av904</td><td>1</td><td>25</td><td>226</td><td>1.00E-177</td><td>9.00E-80</td><td>Bacteria; Bacteroidetes; Cytophagia</td><td>alpha, alpha-trehalase</td></tr></tbody></table>

doi:10.1371/journal.pone.0131313.t001

<footnote>* putative pseudogene: the stop codon in the middle was corrected to predict a full-length protein</footnote>
**Fig 2.** (a.) 3D alignment of the OtsA chain B crystal structure obtained for *E. coli* (white) with the two protein structures of TPS1 modeled using SWISS-MODEL (light blue) and PHYRE2 (orange). UDP (green) and G6P (red) are localized as observed in the OtsA crystal structure. (b.) Close-up on the conserved sites involved in the binding of G6P (red). The residues of the *E. coli* OtsA are shown in white and those of the AvTpsA predicted using SWISS-MODEL in light blue. SwissModel TPS1 *A. thaliana* (light green) and SwissModel predicted TPS1 *P. brassicacea* (gold) are also displayed in order to compare their topologies. (c.) Close-up on the conserved sites related to the binding of UDP (green). The residues are labeled based on the 3D structure of AvTpsA predicted using SWISS-MODEL.

doi:10.1371/journal.pone.0131313.g002

metazoans and bacteria. The TPS sequences of *A. vaga* were part of the clade of plant and functional TPS proteins (strongly supported by bootstrap branch supports of 100-97-84 for RaxML, NJ, and MetaPIGA respectively, and by a posterior probability of 1 using Mr Bayes). The two *A. vaga* AvTpsA and AvTpsA' sequences clustered with the protist TPS sequence of *P. brassicae* (bootstrap supports of 95-76-62, posterior probability of 1) and formed a sister clade to the plant Class I TPS proteins (bootstrap supports of 100-97-87, posterior probability of 1) confirming their non-metazoan origin. Most of these sequences have biochemically proven TPS enzymatic activity, except AtTPS2-AtTPS4 [71]. The position of ΨAvTpsB and ΨAvTpsB' within the tree was not resolved as the phylogenetic trees were not congruent or branch support was insufficient. In contrast, the TPS sequences of monogonont rotifers (Brachionus manjavacas and Brachionus calyciflorus) grouped outside the plant and fungal clade and were associated with other metazoans (bootstrap support values 100-100-99, posterior probability 1). As previously described in [33], the nematode sequences did not cluster with those of other metazoans.

<table><caption>Table 1. Summary of top BLASTP hits for the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes of *Adineta vaga* against GenBank.</caption><thead><tr><th>Gene ID, name</th><th>Scaffold</th><th>Nb of introns</th><th>GC (%)</th><th>Al</th><th>Best hit E-value</th><th>Metazoan best hit E-value</th><th>Best hit, taxonomy</th><th>Best hit, description</th></tr></thead><tbody><tr><td>GSADVT00013141001, AvTpsA</td><td>Av51</td><td>1</td><td>33.2</td><td>165</td><td>1.00E-161</td><td>4.00E-90</td><td>Eukaryota; Protozoa; Rhizaria; Cercozoa</td><td>glycosyltransferase family 20 protein</td></tr><tr><td>GSADVT00013173001, AvTpsA'</td><td>Av51</td><td>1</td><td>33.5</td><td>174</td><td>4.00E-162</td><td>2.00E-86</td><td>Eukaryota; Protozoa; Rhizaria; Cercozoa</td><td>glycosyltransferase family 20 protein</td></tr><tr><td>GSADVT00043720001, ΨAvTpsB*</td><td>Av255</td><td>0</td><td>31.1</td><td>82</td><td>1.00E-107</td><td>9.00E-72</td><td>Eukaryota; Fungi; Chytridiomycota</td><td>glycosyltransferase family 20 protein</td></tr><tr><td>GSADVT00063692001, ΨAvTpsB'*</td><td>Av681</td><td>0</td><td>31.2</td><td>78</td><td>1.00E-102</td><td>8.00E-69</td><td>Eukaryota; Fungi; Dikarya</td><td>glycosyltransferase family 20 protein</td></tr><tr><td>GSADVT00056365001, AvTreA</td><td>Av458</td><td>0</td><td>31.7</td><td>140</td><td>1.00E-147</td><td>7.00E-87</td><td>Bacteria; Bacteroidetes; Sphingobacteria</td><td>alpha, alpha-trehalase</td></tr><tr><td>GSADVT00056372001, AvTreA'</td><td>Av458</td><td>0</td><td>31.7</td><td>140</td><td>1.00E-147</td><td>1.00E-86</td><td>Bacteria; Bacteroidetes; Sphingobacteria</td><td>alpha, alpha-trehalase</td></tr><tr><td>GSADVT00051563001, AvTreB</td><td>Av367</td><td>1</td><td>27</td><td>281</td><td>0.00E+00</td><td>1.00E-78</td><td>Bacteria; Bacteroidetes; Cytophagia</td><td>alpha, alpha-trehalase</td></tr><tr><td>GSADVT00021696001, AvTreB'</td><td>Av920</td><td>1</td><td>26.9</td><td>233</td><td>6.00E-180</td><td>1.00E-78</td><td>Bacteria; Bacteroidetes; Cytophagia</td><td>alpha, alpha-trehalase</td></tr><tr><td>GSADVT00001947001, AvTreC-C'</td><td>Av6</td><td>0</td><td>33.6</td><td>139</td><td>2.00E-149</td><td>4.00E-89</td><td>Bacteria; Bacteroidetes; Cytophagia</td><td>alpha, alpha-trehalase</td></tr><tr><td>N/A, ΨAvTreD*</td><td>Av904</td><td>1</td><td>25</td><td>226</td><td>1.00E-177</td><td>9.00E-80</td><td>Bacteria; Bacteroidetes; Cytophagia</td><td>alpha, alpha-trehalase</td></tr></tbody></table>

doi:10.1371/journal.pone.0131313.t001

<footnote>* putative pseudogene: the stop codon in the middle was corrected to predict a full-length protein</footnote>
<figure><img src="image_13.png" /></figure>
<figure><img src="image_13.png" /></figure>
**Fig 3.** Maximum-likelihood phylogenetic tree based on the amino acid sequences available for trehalose-6 phosphate synthase (TPS) domains. The tree displayed was generated using RaxML with 1000 rapid bootstrap according to the model WAG+I+G+. Numbers above the branches are bootstrap support percentages and posterior probability of 1 obtained with MrBayes (See M&M). When the node was not recovered in the Bayesian tree, the bootstrap value was replaced by "-". Bootstrap support values above 70% for MetaPIGA and NJ are indicated by lines respectively drawn above and below the ML bootstrap values and posterior probabilities*. Color code: green (plant), red (metazoans), black (bacteria and archeae), cyan (protists) and brown (fungi).

doi:10.1371/journal.pone.0131313.g003

## Phylogenetic tree of TRE gene sequences

Unrooted phylogenetic trees (Fig 4) constructed on the basis of the alignment of the amino acid sequences of known or putative TRE predicted the existence of three major clades: one clade of fungal trehalases (bootstrap supports of 100-100-100, posterior probability of 1); one clade of bacterial trehalases (bootstrap supports of 91-100-92, posterior probability of 1); and one clade grouping trehalase sequences of plants, metazoans and some protists (bootstrap supports of 100-100-98, posterior probability of 1). The bacterial clade was further divided into two subclades: one that mostly comprised sequences of the class Bacteroidetes (bootstrap supports of 84-89-85, posterior probability of 1), and a second one that was composed of sequences of the class Proteobacteria (bootstrap supports of 99-100-100, posterior probability of 1). All the seven *A. vaga* trehalase sequences fell within the Bacteroidetes clade, confirming their non-metazoan origin. However, it was not possible to resolve more precisely the position of the *A. vaga* TRE genes inside Bacteroidetes due to incongruence between trees or poor branch support. In contrast, the TRE sequences of the monogonont rotifer *Brachionus calyciflorus* were closely related to those of plants, metazoans and *Dictyostelium discoideum*. Moreover, a close relation between *D. discoideum* and *Brachionus calyciflorus* was observed (bootstrap supports of 88-84-69, posterior probability of 1).

## Expression profiles of the *tps* and *tre* genes

**RNAseq.** To obtain an overview of the expression of the TPS and TRE genes identified in *A. vaga*, we analyzed cDNA libraries representing hydrated *A. vaga* (control condition), early desiccated *A. vaga* (after 37h drying process) and *A. vaga* individuals rehydrated since 1.5 hours following 14 days of desiccation (Figure B in S1 File). Mapping the reads obtained from the cDNA libraries against the putative coding sequences of the four *tps* genes identified in *A. vaga* revealed few matches to *AvTpsA* and *AvTpsA'* (RPKB<1,000 for each gene in each time point) and no match to *ΨAvTpsB* and *ΨAvTpsB'*. For the TRE genes, all the six copies except the putative pseudogene *ΨAvTreD* were expressed at all time points with a high number of matching ESTs (10⁶–10⁸ RPKB, Figure D in S1 File). At the three time points, the number of reads matching the TRE genes was over thousand times greater than the number of reads matching the TPS genes, indicating a much higher expression level for trehalase compared to trehalose-6-phosphate synthase (Figure D in S1 File).

**qPCR.** Quantitative PCRs (qPCRs) were performed on the TPS and TRE genes of *A. vaga* to study in detail their expression kinetics during desiccation. We screened five time points during the drying process, desiccated state and rehydration process and compared them with hydrated (control) bdelloids (Figure B in S1 File). We found an increase of *AvTpsA* and *AvTpsA'* mRNA expression levels during drying that reached a peak of 25 fold at 37 h drying (newly desiccated state) in comparison with the control, then decreased during rehydration (Fig 5A). The difference in relative expression is significant at timepoint 37h of drying.

Similar to the expression pattern of *AvTpsA-A'*, the trehalase copies *AvTreA* and *AvTreA'* were upregulated during drying process with a peak at 37 h and a decrease down to control
**Fig 3.** Maximum-likelihood phylogenetic tree based on the amino acid sequences available for trehalose-6 phosphate synthase (TPS) domains. The tree displayed was generated using RaxML with 1000 rapid bootstrap according to the model WAG+I+G+. Numbers above the branches are bootstrap support percentages and posterior probability of 1 obtained with MrBayes (See M&M). When the node was not recovered in the Bayesian tree, the bootstrap value was replaced by "-". Bootstrap support values above 70% for MetaPIGA and NJ are indicated by lines respectively drawn above and below the ML bootstrap values and posterior probabilities*. Color code: green (plant), red (metazoans), black (bacteria and archeae), cyan (protists) and brown (fungi).

doi:10.1371/journal.pone.0131313.g003

## Phylogenetic tree of TRE gene sequences

Unrooted phylogenetic trees (Fig 4) constructed on the basis of the alignment of the amino acid sequences of known or putative TRE predicted the existence of three major clades: one clade of fungal trehalases (bootstrap supports of 100-100-100, posterior probability of 1); one clade of bacterial trehalases (bootstrap supports of 91-100-92, posterior probability of 1); and one clade grouping trehalase sequences of plants, metazoans and some protists (bootstrap supports of 100-100-98, posterior probability of 1). The bacterial clade was further divided into two subclades: one that mostly comprised sequences of the class Bacteroidetes (bootstrap supports of 84-89-85, posterior probability of 1), and a second one that was composed of sequences of the class Proteobacteria (bootstrap supports of 99-100-100, posterior probability of 1). All the seven *A. vaga* trehalase sequences fell within the Bacteroidetes clade, confirming their non-metazoan origin. However, it was not possible to resolve more precisely the position of the *A. vaga* TRE genes inside Bacteroidetes due to incongruence between trees or poor branch support. In contrast, the TRE sequences of the monogonont rotifer *Brachionus calyciflorus* were closely related to those of plants, metazoans and *Dictyostelium discoideum*. Moreover, a close relation between *D. discoideum* and *Brachionus calyciflorus* was observed (bootstrap supports of 88-84-69, posterior probability of 1).

## Expression profiles of the *tps* and *tre* genes

**RNAseq.** To obtain an overview of the expression of the TPS and TRE genes identified in *A. vaga*, we analyzed cDNA libraries representing hydrated *A. vaga* (control condition), early desiccated *A. vaga* (after 37h drying process) and *A. vaga* individuals rehydrated since 1.5 hours following 14 days of desiccation (Figure B in S1 File). Mapping the reads obtained from the cDNA libraries against the putative coding sequences of the four *tps* genes identified in *A. vaga* revealed few matches to *AvTpsA* and *AvTpsA'* (RPKB<1,000 for each gene in each time point) and no match to *ΨAvTpsB* and *ΨAvTpsB'*. For the TRE genes, all the six copies except the putative pseudogene *ΨAvTreD* were expressed at all time points with a high number of matching ESTs (10⁶–10⁸ RPKB, Figure D in S1 File). At the three time points, the number of reads matching the TRE genes was over thousand times greater than the number of reads matching the TPS genes, indicating a much higher expression level for trehalase compared to trehalose-6-phosphate synthase (Figure D in S1 File).

**qPCR.** Quantitative PCRs (qPCRs) were performed on the TPS and TRE genes of *A. vaga* to study in detail their expression kinetics during desiccation. We screened five time points during the drying process, desiccated state and rehydration process and compared them with hydrated (control) bdelloids (Figure B in S1 File). We found an increase of *AvTpsA* and *AvTpsA'* mRNA expression levels during drying that reached a peak of 25 fold at 37 h drying (newly desiccated state) in comparison with the control, then decreased during rehydration (Fig 5A). The difference in relative expression is significant at timepoint 37h of drying.

Similar to the expression pattern of *AvTpsA-A'*, the trehalase copies *AvTreA* and *AvTreA'* were upregulated during drying process with a peak at 37 h and a decrease down to control
<figure><img src="image_14.png" /><figcaption>**Fig 4.** Maximum likelihood phylogenetic tree based on the amino acid sequences available for trehalase (TRE) domains. The tree displayed was generated using RaxML with 1000 rapid bootstrap according to the LG+I+G evolutionary model. Numbers above the branches are bootstrap support percentages and posterior probability of 1 obtained with MrBayes (See M&M). When the node was not found in a Bayesian tree, the bootstrap value was</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Fig 4.** Maximum likelihood phylogenetic tree based on the amino acid sequences available for trehalase (TRE) domains. The tree displayed was generated using RaxML with 1000 rapid bootstrap according to the LG+I+G evolutionary model. Numbers above the branches are bootstrap support percentages and posterior probability of 1 obtained with MrBayes (See M&M). When the node was not found in a Bayesian tree, the bootstrap value was</figcaption></figure>
replaced by “-”. Bootstrap support values above 70% for MetaPIGA and NJ are indicated by lines respectively drawn above and below the ML bootstrap values and posterior probabilities. Color code: green (plant), red (metazoans), black (bacteria and archeae), cyan (protists) and brown (fungi).

doi:10.1371/journal.pone.0131313.g004

level after rehydration (Fig 5B). At 29h and 37h of drying process the relative expression is sig-
nificantly higher than the control. In contrast to *AvTreA* and *AvTreA'*, *AvTreC-C'* copies did
not show any significant difference between the tested conditions and the control (Fig 5C).
Finally, *AvTreB-B'* appeared upregulated during both drying and rehydration processes in
comparison with the control (Fig 5D).

**Discussion**

In contrast to all previous studies reporting the absence of trehalose and its biosynthetic path-
way in bdelloid rotifers [34–37], our screening of the published genome of *A. vaga* revealed the
presence of trehalose-synthesizing (trehalose-6-phosphate synthase) and hydrolyzing (treha-
lases) genes in the bdelloid rotifer *A. vaga*. Some of the copies of both the synthase and hydro-
lase gene were expressed only at specific time points during the drying-rehydration process
(Fig 5). In addition, the active site residues of the TPS domains were highly conserved in all

<figure><img src="image_15.png" /><figcaption>**Fig 5.** Relative expression of the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes in *Adineta vaga* submitted to drying or rehydration: (a.) AvTpsA&A'; (b.) AvTreA&A”; (c.) AvTreC&C' (d.) AvTreB&B'. The five time points are as follows: (1) hydrated; (2) 28h drying/humidity effect; (3) 37h drying/early desiccated state; (4) 5min rehydration after 14 days of desiccation; (5) 1h30 rehydration after 14 days of desiccation; (6) 3h rehydration after 14 days of desiccation. The expression level of each gene was normalized by reference to the housekeeping gene for the ribosomal protein L40; relative level of expression changes were calculated by reference to that that of the hydrated time point (value = 1). Error bars represent mean value 95% confidence intervals based on three replicates (one-way analysis of variance RM, \*P<0.05 \*\* P<0.01 \*\*\* P<0.001).</figcaption></figure>

doi:10.1371/journal.pone.0131313.g005
replaced by “-”. Bootstrap support values above 70% for MetaPIGA and NJ are indicated by lines respectively drawn above and below the ML bootstrap
values and posterior probabilities. Color code: green (plant), red (metazoans), black (bacteria and archeae), cyan (protists) and brown (fungi).

doi:10.1371/journal.pone.0131313.g004

level after rehydration (Fig 5B). At 29h and 37h of drying process the relative expression is sig-
nificantly higher than the control. In contrast to AvTreA and AvTreA', AvTreC-C' copies did
not show any significant difference between the tested conditions and the control (Fig 5C).
Finally, AvTreB-B' appeared upregulated during both drying and rehydration processes in
comparison with the control (Fig 5D).

Discussion

In contrast to all previous studies reporting the absence of trehalose and its biosynthetic path-
way in bdelloid rotifers [34–37], our screening of the published genome of A. vaga revealed the
presence of trehalose-synthesizing (trehalose-6-phosphate synthase) and hydrolyzing (treha-
lases) genes in the bdelloid rotifer A. vaga. Some of the copies of both the synthase and hydro-
lase gene were expressed only at specific time points during the drying-rehydration process
(Fig 5.). In addition, the active site residues of the TPS domains were highly conserved in all

<figure><img src="image_15.png" /><figcaption>Fig 5. Relative expression of the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes in *Adineta vaga* submitted to drying or rehydration: (a.) AvTpsA&A'; (b.) AvTreA&A”; (c.) AvTreC&C' (d.) AvTreB&B'. The five time points are as follows: (1) hydrated; (2) 28h drying/humidity effect; (3) 37h drying/early desiccated state; (4) 5min rehydration after 14 days of desiccation; (5) 1h30 rehydration after 14 days of desiccation; (6) 3h rehydration after 14 days of desiccation. The expression level of each gene was normalized by reference to the housekeeping gene for the ribosomal protein L40; relative level of expression changes were calculated by reference to that that of the hydrated time point (value = 1). Error bars represent mean value 95% confidence intervals based on three replicates (one-way analysis of variance RM, \*P<0.05 \*\* P<0.01 \*\*\* P<0.001).</figcaption></figure>

doi:10.1371/journal.pone.0131313.g005
the four TPS gene copies of *A. vaga* (Figure C in S1 File), especially in the *AvTpsA* pair. Moreover the modeled 3D structures of *AvTpsA* suggest that the organization of the active site is also conserved (Fig 2). The genes however appear to be acquired by HGT and as in plants, the trehalose-6-phosphate synthesized by the TPS gene might serve as a signaling molecule when bdelloids become desiccated (see below).

## Evidence for the acquisitions of the *tps* and *tre* genes of *A. vaga* by horizontal gene transfer

The four *A. vaga* TPS genes have AI>45 and are therefore of probable non-metazoan origin. This was confirmed using different phylogenetic methods: all four *A. vaga* TPS sequences fell within the plant-fungal TPS clade (including the protist sequences of *P. brassicae*) with strong support (Fig 3). *AvTpsA-A'* formed a sister clade to the clade grouping plant Class I TPS proteins and the protist sequence of *P. brassicae*. The TPS genes Ψ*AvTpsB* and Ψ*AvTpsB'* appeared to form a sister group to the plant and fungi Class I/ II TPS clade (including *P. brassicae*) but the origin of these *A. vaga* TPS domain sequences remains unresolved (Fig 3). The presence of stop codons at two different positions in Ψ*AvTpsB* and Ψ*AvTpsB'* suggests that pseudogenization happened after their acquisition. In the absence of selective pressure because of a plausible redundancy with *AvTpsA-A'*, Ψ*AvTpsB* and Ψ*AvTpsB'* may have diverged considerably from all other TPS genes studied and their position within the plant- fungi clade is therefore unresolved.

The two copies within each pair (*AvTpsA* and *AvTpsA'*, Ψ*AvTpsB* and Ψ*AvTpsB'* ) are nearly colinear but *AvTpsA-A'* does not appear to be ohnologous to *AvTpsB-B'*. Therefore, we hypothesize that bdelloid rotifers acquired two different *tps* gene copies (A and B) by two independent horizontal gene transfer events. These two independently acquired *tps* genes may have become duplicated during a gene conversion event or at the time of the whole genome duplication event that is supposed to have occurred early in bdelloid rotifer evolution [5,35], resulting in the two colinear pairs we found. However, we cannot exclude that this duplication occurred in a sexual ancestor of extant bdelloids, in which case, meiotic recombination may have been involved in the transition from one copy to two copies of each of the two genes. Determining whether this acquisition was recent or very ancient among bdelloids will require comparative analysis with other bdelloid species.

The seven *A. vaga* trehalase sequences fell within the bacterial TRE clade (Fig 4), and all the trehalase gene copies had AI>45 with one or no intron, indicating that the TRE genes of *A. vaga* may have been acquired by horizontal gene transfer from bacteria. The putative horizontal gene transfers involving TRE genes appear relatively ancient, since (1) the trehalase gene copies *AvTreA-A'* and *AvTreC-C'* on scaffolds av458 and av6 are organized in a quartet of four homologous regions (Fig 1B), and (2) the three TRE gene copies *AvTreB*, *AvTreB'* and Ψ*AvTreD* on scaffolds av367, av920 and av904 have acquired one intron each (Table 1).

Horizontal gene transfer appears to be rampant in bdelloid rotifers [4,39,46], and some genes of foreign origin were previously reported to be intact and expressed, as for the TPS and TRE genes in our study. Boschetti et al. [37], studying gene expression in response to desiccation, did not detect expression of any trehalose gene probably because the critical time point of 37h drying (at which the rotifer can be considered as having reached anhydrobiosis as its water content is reduced to about 5% [64]) was not included in their study. The authors however confirmed the existence of a high proportion (8–9%) of foreign genes in the bdelloid genome that appear to be transcribed. Those results suggest that bdelloid rotifers can capture exogenous genes and that these appear to be functional. Horizontal gene transfers in bdelloid rotifers most
the four TPS gene copies of *A. vaga* (Figure C in S1 File), especially in the *AvTpsA* pair. Moreover the modeled 3D structures of *AvTpsA* suggest that the organization of the active site is also conserved (Fig 2). The genes however appear to be acquired by HGT and as in plants, the trehalose-6-phosphate synthesized by the TPS gene might serve as a signaling molecule when bdelloids become desiccated (see below).

## Evidence for the acquisitions of the *tps* and *tre* genes of *A. vaga* by horizontal gene transfer

The four *A. vaga* TPS genes have AI>45 and are therefore of probable non-metazoan origin. This was confirmed using different phylogenetic methods: all four *A. vaga* TPS sequences fell within the plant-fungal TPS clade (including the protist sequences of *P. brassicae*) with strong support (Fig 3). *AvTpsA-A'* formed a sister clade to the clade grouping plant Class I TPS proteins and the protist sequence of *P. brassicae*. The TPS genes Ψ*AvTpsB* and Ψ*AvTpsB'* appeared to form a sister group to the plant and fungi Class I/ II TPS clade (including *P. brassicae*) but the origin of these *A. vaga* TPS domain sequences remains unresolved (Fig 3). The presence of stop codons at two different positions in Ψ*AvTpsB* and Ψ*AvTpsB'* suggests that pseudogenization happened after their acquisition. In the absence of selective pressure because of a plausible redundancy with *AvTpsA-A'*, Ψ*AvTpsB* and Ψ*AvTpsB'* may have diverged considerably from all other TPS genes studied and their position within the plant- fungi clade is therefore unresolved.

The two copies within each pair (*AvTpsA* and *AvTpsA'*, Ψ*AvTpsB* and Ψ*AvTpsB'* ) are nearly colinear but *AvTpsA-A'* does not appear to be ohnologous to *AvTpsB-B'*. Therefore, we hypothesize that bdelloid rotifers acquired two different *tps* gene copies (A and B) by two independent horizontal gene transfer events. These two independently acquired *tps* genes may have become duplicated during a gene conversion event or at the time of the whole genome duplication event that is supposed to have occurred early in bdelloid rotifer evolution [5,35], resulting in the two colinear pairs we found. However, we cannot exclude that this duplication occurred in a sexual ancestor of extant bdelloids, in which case, meiotic recombination may have been involved in the transition from one copy to two copies of each of the two genes. Determining whether this acquisition was recent or very ancient among bdelloids will require comparative analysis with other bdelloid species.

The seven *A. vaga* trehalase sequences fell within the bacterial TRE clade (Fig 4), and all the trehalase gene copies had AI>45 with one or no intron, indicating that the TRE genes of *A. vaga* may have been acquired by horizontal gene transfer from bacteria. The putative horizontal gene transfers involving TRE genes appear relatively ancient, since (1) the trehalase gene copies *AvTreA-A'* and *AvTreC-C'* on scaffolds av458 and av6 are organized in a quartet of four homologous regions (Fig 1B), and (2) the three TRE gene copies *AvTreB*, *AvTreB'* and Ψ*AvTreD* on scaffolds av367, av920 and av904 have acquired one intron each (Table 1).

Horizontal gene transfer appears to be rampant in bdelloid rotifers [4,39,46], and some genes of foreign origin were previously reported to be intact and expressed, as for the TPS and TRE genes in our study. Boschetti et al. [37], studying gene expression in response to desiccation, did not detect expression of any trehalose gene probably because the critical time point of 37h drying (at which the rotifer can be considered as having reached anhydrobiosis as its water content is reduced to about 5% [64]) was not included in their study. The authors however confirmed the existence of a high proportion (8–9%) of foreign genes in the bdelloid genome that appear to be transcribed. Those results suggest that bdelloid rotifers can capture exogenous genes and that these appear to be functional. Horizontal gene transfers in bdelloid rotifers most
likely occur during desiccation when membranes are disrupted and DNA is broken and then repaired [64,72,73].

## Why did previous studies fail to detect TPS genes in bdelloid rotifers?

Previous studies using degenerate primers were not able to amplify any TPS fragment in the three bdelloid species tested [34,35]. The degenerate primers used in those studies match poorly the four TPS gene sequences of the present study: the first eight nucleotides of the forward primer were different from the *A. vaga* TPS gene sequences, whereas numerous differences were observed between the reverse primer and the *A. vaga* TPS gene sequences (Figure E in S1 File). Indeed, Lapinski and Tunnacliffe [34] defined their primers based on an alignment of TPS sequences from bacteria, yeast and metazoans, whereas the TPS genes of *A. vaga* are more closely related to plants for *AvTpsA* and *AvTpsA'* and to fungi or plant for Ψ*AvTpsB* and Ψ*AvTpsB'*.

Previous studies also failed to detect TPS gene expression as well as trehalose accumulation. At each time point, we observed a much higher expression level for the TRE genes than for TPS (Figure D in S1 File): this strong trehalase activity in *A. vaga* may explain why previous studies could not detect trehalose [34,36]. Moreover, we found the expression of *AvTpsA-A'* and *AvTreA-A'* to be significantly upregulated during drying process, reaching a peak at 37 h (when bdelloids had just entered into desiccated state) then decreasing to the control level during rehydration (Fig 5A and 5B). *AvTreB-B'* were significantly up-regulated in both drying and rehydration compared with the control (Fig 5D), indicating that those gene copies may be constitutively expressed with an upregulation at desiccation. In any case, the low level of transcription observed for the TPS genes of *A. vaga* and their expression at very specific time points may explain why previous studies failed to detect trehalose or the expression of the genes involved it its biosynthesis. Interestingly, this biological pattern is very similar to the trehalose metabolism/catabolism pathway observed in plants, where the trehalase activity is so high that overexpression of TPS genes does not result in an increase in trehalose unless trehalase activity is disrupted [31,74].

## The absence of trehalose phosphatase in bdelloid rotifers

No fusion of TPP domains with our TPS *A. vaga* genes was observed. For *AvTpsA-A'*, we hypothesize that the absence of TPP domains could be linked to their non-metazoan origin: indeed, the plant Class I TPS genes to which they are related do not contain any TPP active site [31].

Neither trehalose-6-phosphate phosphatase nor conserved motifs of the TPP domain were detected in the genome of *A. vaga*. Plant, bacteria, fungi and yeast cells contain a wide array of nonspecific phosphatases that are able to dephosphorylate trehalose-6 phosphate (T6P), although with a lower efficiency than TPP [14,75-77]. Therefore, it is possible that in *A. vaga* the trehalose biosynthesis pathway relies on such nonspecific phosphatases instead of TPP.

In *E. coli*, TPS and TPP enzymes are encoded by distinct genes (OtsA and OtsB) [33,78] (Figure A in S1 File), while in almost all eukaryotes, the TPS proteins are fused to the TPP domains (Figure A in S1 File) [33,69]. This organization suggests that all the eukaryotic TPS and TPP fused proteins descend from a common ancestor [33]. Avonce et al. 2010 [71] recently demonstrated the existence of a natural TPS-TPP bifunctional enzyme in the bacterial species *Cytophaga hutchinsonii*. Phylogenetic analysis further suggested that this prokaryotic gene might be related to the eukaryotic TPS-TPP fused genes. Among the four predicted TPS protein sequences of *A. vaga*, only *AvTpsA-A'* was susceptible to contain putative TPP regions in their extended C-terminus but none of the conserved short motifs of the TPP-like domain
likely occur during desiccation when membranes are disrupted and DNA is broken and then repaired [64,72,73].

## Why did previous studies fail to detect TPS genes in bdelloid rotifers?

Previous studies using degenerate primers were not able to amplify any TPS fragment in the three bdelloid species tested [34,35]. The degenerate primers used in those studies match poorly the four TPS gene sequences of the present study: the first eight nucleotides of the forward primer were different from the *A. vaga* TPS gene sequences, whereas numerous differences were observed between the reverse primer and the *A. vaga* TPS gene sequences (Figure E in S1 File). Indeed, Lapinski and Tunnacliffe [34] defined their primers based on an alignment of TPS sequences from bacteria, yeast and metazoans, whereas the TPS genes of *A. vaga* are more closely related to plants for *AvTpsA* and *AvTpsA'* and to fungi or plant for Ψ*AvTpsB* and Ψ*AvTpsB'*.

Previous studies also failed to detect TPS gene expression as well as trehalose accumulation. At each time point, we observed a much higher expression level for the TRE genes than for TPS (Figure D in S1 File): this strong trehalase activity in *A. vaga* may explain why previous studies could not detect trehalose [34,36]. Moreover, we found the expression of *AvTpsA-A'* and *AvTreA-A'* to be significantly upregulated during drying process, reaching a peak at 37 h (when bdelloids had just entered into desiccated state) then decreasing to the control level during rehydration (Fig 5A and 5B). *AvTreB-B'* were significantly up-regulated in both drying and rehydration compared with the control (Fig 5D), indicating that those gene copies may be constitutively expressed with an upregulation at desiccation. In any case, the low level of transcription observed for the TPS genes of *A. vaga* and their expression at very specific time points may explain why previous studies failed to detect trehalose or the expression of the genes involved it its biosynthesis. Interestingly, this biological pattern is very similar to the trehalose metabolism/catabolism pathway observed in plants, where the trehalase activity is so high that overexpression of TPS genes does not result in an increase in trehalose unless trehalase activity is disrupted [31,74].

## The absence of trehalose phosphatase in bdelloid rotifers

No fusion of TPP domains with our TPS *A. vaga* genes was observed. For *AvTpsA-A'*, we hypothesize that the absence of TPP domains could be linked to their non-metazoan origin: indeed, the plant Class I TPS genes to which they are related do not contain any TPP active site [31].

Neither trehalose-6-phosphate phosphatase nor conserved motifs of the TPP domain were detected in the genome of *A. vaga*. Plant, bacteria, fungi and yeast cells contain a wide array of nonspecific phosphatases that are able to dephosphorylate trehalose-6 phosphate (T6P), although with a lower efficiency than TPP [14,75-77]. Therefore, it is possible that in *A. vaga* the trehalose biosynthesis pathway relies on such nonspecific phosphatases instead of TPP.

In *E. coli*, TPS and TPP enzymes are encoded by distinct genes (OtsA and OtsB) [33,78] (Figure A in S1 File), while in almost all eukaryotes, the TPS proteins are fused to the TPP domains (Figure A in S1 File) [33,69]. This organization suggests that all the eukaryotic TPS and TPP fused proteins descend from a common ancestor [33]. Avonce et al. 2010 [71] recently demonstrated the existence of a natural TPS-TPP bifunctional enzyme in the bacterial species *Cytophaga hutchinsonii*. Phylogenetic analysis further suggested that this prokaryotic gene might be related to the eukaryotic TPS-TPP fused genes. Among the four predicted TPS protein sequences of *A. vaga*, only *AvTpsA-A'* was susceptible to contain putative TPP regions in their extended C-terminus but none of the conserved short motifs of the TPP-like domain
were detected [33,68], while the other two sequences (ΨΑνΤρςΒ and ΨΑνΤρςΒ') clearly did not include TPP domains. Previous research have shown that although TPP-like domains have lost their TPP catalytic activity, these domains may still be required for correct protein conformation and stability as in the early bi-functional fusion enzyme [69].

**The plausible role of trehalose metabolism in bdelloid rotifers.** Trehalose is synthesized in many anhydrobiotic organisms and was once hypothesized to be a universal molecule that preserves cell integrity during desiccation. Nowadays, low trehalose levels found in some desiccation resistant animals like tardigrades or monogonont eggs, and in some plants, suggest that trehalose does not only act as bioprotective molecule but could be used alternatively as a signaling molecule.

In this study, we identified a series of genes in *A. vaga* that encode enzymes involved in trehalose synthesis and hydrolysis, some of which are substantially expressed. Although neither trehalose-6-phosphate phosphatase (TPP) nor conserved motifs of the TPP domain have been identified in *A. vaga*, we cannot exclude the possibility that *A. vaga* might synthesize very low amounts of endogenous trehalose using unspecific phosphatases. Interestingly, several features of the *A. vaga* trehalose metabolism, acquired through HGT, are highly similar to those observed in plants, such as the lack of trehalose accumulation, the existence of multiple copies of TPS and TRE genes, and the higher expression levels of TRE compared with TPS genes. In plants, both T6P and trehalose act as signaling molecules that induce metabolic changes, such as the accumulation of storage carbohydrates. As for many signaling molecules, the rapid degradation of trehalose may therefore be required to prevent its accumulation from interfering with the regulation of plant metabolism [15,79]. Furthermore, most Class I and Class II TPS proteins in plants have been suspected to have lost their TPS and TPP enzymatic activities but can still be detected, suggesting that other selective pressures are responsible for their persistence such as their role in other plant carbohydrate-related pathways involved in the development, sugar sensing and abiotic stress response [15,69,71]. Moreover, the synthesis of the trehalose precursor, trehalose-6-phosphate (T6P), has been shown to play a signaling role in plant growth and development and in sugar metabolism (sugar sensing under stress tolerance) [80–82]. T6P is also known as a signaling molecule in the yeast species *Saccharomyces cerevisiae* [83]. Therefore, there is growing evidence that it is not trehalose itself but more likely one of its precursors (such as T6P), or perhaps an enzyme involved in its synthesis (such as TPS), that serves as a signaling molecule controlling certain metabolic pathways [82,84]. However, it was also noticed that despite his key role in plant metabolism, T6P is present in very small quantities [85]. Our discovery of multiple copies of low-expression trehalose synthesis genes and highly expressed trehalose degradation genes in *A. vaga* suggests that trehalose and/or trehalose-6-phosphate may fulfill a signaling role in bdelloid rotifers as was proposed in plants and fungi, from which bdelloids may have acquired the TPS genes. This signaling role hypothesis is supported by the overexpression of AvTpsA-A' in early desiccated *A. vaga*, suggesting that T6P is involved as a messenger during dehydration. However, it remains necessary to confirm these hypotheses by measuring directly the appropriate metabolites and/ or enzyme activities. Finally, a recent study suggested that plant TPS may have scaffolding functions or can be involved in protein complexes [86]. Moreover, TPS from *Magnaporthe oryzae* was described as a key protein involved in enzyme activation, control of metabolic flux and transcriptional regulation of pathogenic fungus [87]. Therefore, we cannot exclude that bdelloids TPS have both catalytic and non-catalytic functions involved in the desiccation process. The discovery of TPS expression in bdelloids during their desiccation invites us to re-investigate the apparent paradox of trehalose's absence in these organisms.
were detected [33,68], while the other two sequences (ΨAvTpsB and ΨAvTpsB') clearly did not include TPP domains. Previous research have shown that although TPP-like domains have lost their TPP catalytic activity, these domains may still be required for correct protein conformation and stability as in the early bi-functional fusion enzyme [69].

**The plausible role of trehalose metabolism in bdelloid rotifers.** Trehalose is synthesized in many anhydrobiotic organisms and was once hypothesized to be a universal molecule that preserves cell integrity during desiccation. Nowadays, low trehalose levels found in some desiccation resistant animals like tardigrades or monogonont eggs, and in some plants, suggest that trehalose does not only act as bioprotective molecule but could be used alternatively as a signaling molecule.

In this study, we identified a series of genes in *A. vaga* that encode enzymes involved in trehalose synthesis and hydrolysis, some of which are substantially expressed. Although neither trehalose-6-phosphate phosphatase (TPP) nor conserved motifs of the TPP domain have been identified in *A. vaga*, we cannot exclude the possibility that *A. vaga* might synthesize very low amounts of endogenous trehalose using unspecific phosphatases. Interestingly, several features of the *A. vaga* trehalose metabolism, acquired through HGT, are highly similar to those observed in plants, such as the lack of trehalose accumulation, the existence of multiple copies of TPS and TRE genes, and the higher expression levels of TRE compared with TPS genes. In plants, both T6P and trehalose act as signaling molecules that induce metabolic changes, such as the accumulation of storage carbohydrates. As for many signaling molecules, the rapid degradation of trehalose may therefore be required to prevent its accumulation from interfering with the regulation of plant metabolism [15,79]. Furthermore, most Class I and Class II TPS proteins in plants have been suspected to have lost their TPS and TPP enzymatic activities but can still be detected, suggesting that other selective pressures are responsible for their persistence such as their role in other plant carbohydrate-related pathways involved in the development, sugar sensing and abiotic stress response [15,69,71]. Moreover, the synthesis of the trehalose precursor, trehalose-6-phosphate (T6P), has been shown to play a signaling role in plant growth and development and in sugar metabolism (sugar sensing under stress tolerance) [80–82]. T6P is also known as a signaling molecule in the yeast species *Saccharomyces cerevisiae* [83]. Therefore, there is growing evidence that it is not trehalose itself but more likely one of its precursors (such as T6P), or perhaps an enzyme involved in its synthesis (such as TPS), that serves as a signaling molecule controlling certain metabolic pathways [82,84]. However, it was also noticed that despite his key role in plant metabolism, T6P is present in very small quantities [85]. Our discovery of multiple copies of low-expression trehalose synthesis genes and highly expressed trehalose degradation genes in *A. vaga* suggests that trehalose and/or trehalose-6-phosphate may fulfill a signaling role in bdelloid rotifers as was proposed in plants and fungi, from which bdelloids may have acquired the TPS genes. This signaling role hypothesis is supported by the overexpression of AvTpsA-A' in early desiccated *A. vaga*, suggesting that T6P is involved as a messenger during dehydration. However, it remains necessary to confirm these hypotheses by measuring directly the appropriate metabolites and/ or enzyme activities. Finally, a recent study suggested that plant TPS may have scaffolding functions or can be involved in protein complexes [86]. Moreover, TPS from *Magnaporthe oryzae* was described as a key protein involved in enzyme activation, control of metabolic flux and transcriptional regulation of pathogenic fungus [87]. Therefore, we cannot exclude that bdelloids TPS have both catalytic and non-catalytic functions involved in the desiccation process. The discovery of TPS expression in bdelloids during their desiccation invites us to re-investigate the apparent paradox of trehalose's absence in these organisms.
**Supporting Information**

**S1 Dataset. qPCR data and statistical tests (assumptions of normality and homoscedasticity) performed for TPS and TRE genes.**

(XLSX)

**S1 File.** Figure A: The metabolic pathways of trehalose: 1. Five distinct trehalose synthesis pathways are present in eukaryotes, bacteria and archaea. 2. Three distinct trehalose hydrolytic pathways have been described. UDP-Glc, uridine diphosphate glucose; Glc-6P, glucose-6-phosphate; T6P, trehalose-6-phosphate; UDP, uridine diphosphate; Glc-1P, glucose-1-phosphate; Glc, glucose; ADP, adenosine triphosphate (modified from [88]). Figure B: Expression study of the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes in the bdelloid rotifer *Adineta vaga* during desiccation and rehydration: the six different time points included in the qPCR analyses are given. The cDNA libraries (in triplicates) were performed at the time points marked with "*". Figure C: Protein sequence alignment of *E. coli* OtsA with the four *Adineta vaga* TPS domains, *A. thaliana* TPS1 (GI:15218422) and *P. brassicae* TPS1 (GI:160332824). Residues important for G6P and UDP binding in *E. coli* OtsA are indicated in red or green, respectively. Mutation compared with the OtsA model are shown in dark red and dark green for residues involved in binding G6P or UDP respectively. Figure D: Quantitative expression of the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes in *A. vaga* submitted to drying and rehydration (based on triplicate RNAseq libraries). Error bars represent standard deviations. ΨAvTpsB and ΨAvTpsB' were not represented as there were no corresponding reads in the libraries. Figure E: Alignment (using DNASTAR) of degenerate (a) forward and (b) reverse primers with the four trehalose-6-phosphate genes of *A. vaga*.

(PDF)

**S2 File.** Table A: Primers used for amplifying the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes of *Adineta vaga* using PCR. Table B: Reference for sequences used in phylogeny related to trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes.

Table C: Primers used for amplifying the trehalose-6-phosphate synthase (TPS), trehalase (TRE) and L40 genes of *A. vaga* using qPCR. Table D: Nucleotidic and proteic similarities between *A. vaga* TRE sequences. The numbers in each cell represent respectively the nucleotidic and proteic similarity percentages.

(PDF)

**Acknowledgments**

The authors are grateful to E. Etoundi and J. Marescaux (UNamur) for useful discussions concerning the phylogenetic part of the project. Thanks also to E. Tinti (UNamur) for his comments and expertise in 3D structure predictions. M. Knapen is highly acknowledged for her technical help during the initiation phase of this project; special thanks to D. Mark Welch for sharing monogonont sequences.

**Author Contributions**

Conceived and designed the experiments: BH XL KVD. Performed the experiments: BH XL JM CDS. Analyzed the data: BH XL JFF LMP JM CDS KVD. Contributed reagents/materials/ analysis tools: BH XL LMP CDS JM KVD. Wrote the paper: BH XL JFF KVD.
**Supporting Information**

**S1 Dataset. qPCR data and statistical tests (assumptions of normality and homoscedasticity) performed for TPS and TRE genes.**

(XLSX)

**S1 File. Figure A: The metabolic pathways of trehalose: 1.** Five distinct trehalose synthesis pathways are present in eukaryotes, bacteria and archaea. 2. Three distinct trehalose hydrolytic pathways have been described. UDP-Glc, uridine diphosphate glucose; Glc-6P, glucose-6-phosphate; T6P, trehalose-6-phosphate; UDP, uridine diphosphate; Glc-1P, glucose-1-phosphate; Glc, glucose; ADP, adenosine triphosphate (modified from [88]). Figure B: Expression study of the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes in the bdelloid rotifer *Adineta vaga* during desiccation and rehydration: the six different time points included in the qPCR analyses are given. The cDNA libraries (in triplicates) were performed at the time points marked with "*". Figure C: Protein sequence alignment of *E. coli* OtsA with the four *Adineta vaga* TPS domains, *A. thaliana* TPS1 (GI:15218422) and *P. brassicae* TPS1 (GI:160332824). Residues important for G6P and UDP binding in *E. coli* OtsA are indicated in red or green, respectively. Mutation compared with the OtsA model are shown in dark red and dark green for residues involved in binding G6P or UDP respectively. Figure D: Quantitative expression of the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes in *A. vaga* submitted to drying and rehydration (based on triplicate RNAseq libraries). Error bars represent standard deviations. ΨAvTpsB and ΨAvTpsB' were not represented as there were no corresponding reads in the libraries. Figure E: Alignment (using DNASTAR) of degenerate (a) forward and (b) reverse primers with the four trehalose-6-phosphate genes of *A. vaga*.

(PDF)

**S2 File. Table A: Primers used for amplifying the trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes of *Adineta vaga* using PCR.** Table B: Reference for sequences used in phylogeny related to trehalose-6-phosphate synthase (TPS) and trehalase (TRE) genes.

Table C: Primers used for amplifying the trehalose-6-phosphate synthase (TPS), trehalase (TRE) and L40 genes of *A. vaga* using qPCR. Table D: Nucleotidic and proteic similarities between *A. vaga* TRE sequences. The numbers in each cell represent respectively the nucleotidic and proteic similarity percentages.

(PDF)

**Acknowledgments**

The authors are grateful to E. Etoundi and J. Marescaux (UNamur) for useful discussions concerning the phylogenetic part of the project. Thanks also to E. Tinti (UNamur) for his comments and expertise in 3D structure predictions. M. Knapen is highly acknowledged for her technical help during the initiation phase of this project; special thanks to D. Mark Welch for sharing monogonont sequences.

**Author Contributions**

Conceived and designed the experiments: BH XL KVD. Performed the experiments: BH XL JM CDS. Analyzed the data: BH XL JFF LMP JM CDS KVD. Contributed reagents/materials/ analysis tools: BH XL LMP CDS JM KVD. Wrote the paper: BH XL JFF KVD.
References

1. Hsu WS (1956) Oogenesis in Habrotrocha tridens (Milne). Biol Bull 111: 364.

2. Hsu WS (1956) Oogenesis in the Bdelloidea rotifer *Philodina roseola* Ehrenberg. Cellule 57: 283-296.

3. Boschetti C, Carr A, Crisp A, Eyres I, Wang-Koh Y, Lubzens E, et al. (2012) Biochemical diversification through foreign gene expression in bdelloid rotifers. PLoS Genet 8: e1003035. doi: 10.1371/journal.pgen.1003035 PMID: 23166508

4. Hespeels B, Flot J-F, Derzelle A, Van Doninck K (2014) Evidence for ancient horizontal gene acquisi-
tions in bdelloid rotifers of the genus Adineta. In: Pontarotti P, editor. Evolutionary Biology: Genome
Evolution, Speciation, Coevolution and Origin of Life. Springer International Publishing. pp. 207-225.
doi: 10.1007/978-3-319-07623-2

5. Segers H (2007) Annotated checklist of the rotifers (Phylum Rotifera), with notes on nomenclature, tax-
onomy and distribution. Zootaxa 1564: 104.

6. Ricci C (1987) Ecology of bdelloids: how to be successful. Hydrobiologia 147: 117-127. doi: 10.1007/BF00025734

7. Van Leeuwenhoek A (1702) Green Weeds growing in Water, and some Animalcula found about them.
Philos Trans 23: 1304-1311.

8. Denekamp N, Thorne M, Clark M, Kube M, Reinhardt R, Lubzens E (2009) Discovering genes associ-
ated with dormancy in the monogonont rotifer *Brachionus plicatilis*. BMC Genomics 10: 108. doi: 10.
1186/1471-2164-10-108 PMID: 19284654

9. Hengherr S, Heyer AG, Brümmer F, Schill RO (2011) Trehalose and vitreous states: desiccation tolerance of dormant stages of the crustaceans Triops and Daphnia. Physiol Biochem Zool 84: 147-153. doi: 10.1086/658499 PMID: 21460525

10. Hengherr S, Schill RO, Clegg JS (2011) Mechanisms associated with cellular desiccation tolerance in the animal extremophile Artemia. Physiol Biochem Zool 84: 249-257. doi: 10.1086/659314 PMID: 21527815

11. Gusev O, Nakahara Y, Vanyagina V, Malutina L, Cornette R, Sakashita T, et al. (2010) Anhydrobiosis-associated nuclear DNA damage and repair in the sleeping chironomid: Linkage with radioresistance. PLoS One 5: e14008. doi: 10.1371/journal.pone.0014008 PMID: 21103355

12. Alpert P (2005) The limits and frontiers of desiccation-tolerant life. Integr Comp Biol 45: 685-695. doi: 10.1093/icb/45.5.685 PMID: 21676818

13. Potts M, Slaughter SM, Hunneke F-U, Garst JF, Helm RF (2005) Desiccation tolerance of prokaryotes: application of principles to human cells. Integr Comp Biol 45: 800-809. doi: 10.1093/icb/45.5.800 PMID: 21676831

14. Elbein AD, Pan YT, Pastuszak I, Carroll D (2003) New insights on trehalose: a multifunctional molecule. Glycobiology 13: 17R-27R. doi: 10.1093/glycob/cwg047 PMID: 12626396

15. Fernandez O, Béthencourt L, Quero A, Sangwan RS, Clément C (2010) Trehalose and plant stress responses: friend or foe? Trends Plant Sci 15: 409-417. doi: 10.1016/j.tplants.2010.04.004 PMID: 20494608

16. Ohtake S, Wang Y (2011) Trehalose: current use and future applications. J Pharm Sci 100: 2020-2053. doi: 10.1002/jps PMID: 21337544

17. Potts M (2001) Desiccation tolerance: a simple process? Trends Microbiol 9: 553-559. PMID: 11825716

18. Jain NK, Roy I (2009) Effect of trehalose on protein structure. Protein Sci 18: 24-36. doi: 10.1002/pro.3 PMID: 19177348

19. Zaparty M, Hagemann A, Bräsen C, Hensel R, Lupas AN, Brinkmann H, et al. (2013) The first prokaryotic trehalose synthase complex identified in the hyperthermophilic crenarchaeon Thermoproteus tenax. PLoS One 8: e61354. doi: 10.1371/journal.pone.0061354 PMID: 23626675

20. Pampurova S, Verschooten K, Avonce N, Van Dijck P (2014) Functional screening of a cDNA library from the desiccation-tolerant plant Selaginella lepidophylla in yeast mutants identifies trehalose biosynthesis genes of plant and microbial origin. J Plant Res 127: 803-813. doi: 10.1007/s10265-014-0663-x PMID: 25246071

21. Sakurai M, Furuki T, Akao K-I, Tanaka D, Nakahara Y, Kikawada T, et al. (2008) Vitrification is essential for anhydrobiosis in an African chironomid, Polypedilum vanderplanki. Proc Natl Acad Sci U S A 105: 5093-5098. doi: 10.1073/pnas.0706197105 PMID: 18362351

22. Madin KAC, Crowe JH (1975) Anhydrobiosis in nematodes: Carbohydrate and lipid metabolism during dehydration. J Exp Zool 193: 335-342. doi: 10.1002/jez.1401930309
References

1. Hsu WS (1956) Oogenesis in Habrotrocha tridens (Milne). Biol Bull 111: 364.

2. Hsu WS (1956) Oogenesis in the Bdelloidea rotifer *Philodina roseola* Ehrenberg. Cellule 57: 283-296.

3. Boschetti C, Carr A, Crisp A, Eyres I, Wang-Koh Y, Lubzens E, et al. (2012) Biochemical diversification through foreign gene expression in bdelloid rotifers. PLoS Genet 8: e1003035. doi: 10.1371/journal.pgen.1003035 PMID: 23166508

4. Hespeels B, Flot J-F, Derzelle A, Van Doninck K (2014) Evidence for ancient horizontal gene acquisi-
tions in bdelloid rotifers of the genus Adineta. In: Pontarotti P, editor. Evolutionary Biology: Genome
Evolution, Speciation, Coevolution and Origin of Life. Springer International Publishing. pp. 207-225.
doi: 10.1007/978-3-319-07623-2

5. Segers H (2007) Annotated checklist of the rotifers (Phylum Rotifera), with notes on nomenclature, tax-
onomy and distribution. Zootaxa 1564: 104.

6. Ricci C (1987) Ecology of bdelloids: how to be successful. Hydrobiologia 147: 117-127. doi: 10.1007/BF00025734

7. Van Leeuwenhoek A (1702) Green Weeds growing in Water, and some Animalcula found about them.
Philos Trans 23: 1304-1311.

8. Denekamp N, Thorne M, Clark M, Kube M, Reinhardt R, Lubzens E (2009) Discovering genes associ-
ated with dormancy in the monogonont rotifer *Brachionus plicatilis*. BMC Genomics 10: 108. doi: 10.
1186/1471-2164-10-108 PMID: 19284654

9. Hengherr S, Heyer AG, Brümmer F, Schill RO (2011) Trehalose and vitreous states: desiccation tolerance of dormant stages of the crustaceans Triops and Daphnia. Physiol Biochem Zool 84: 147-153. doi: 10.1086/658499 PMID: 21460525

10. Hengherr S, Schill RO, Clegg JS (2011) Mechanisms associated with cellular desiccation tolerance in the animal extremophile Artemia. Physiol Biochem Zool 84: 249-257. doi: 10.1086/659314 PMID: 21527815

11. Gusev O, Nakahara Y, Vanyagina V, Malutina L, Cornette R, Sakashita T, et al. (2010) Anhydrobiosis-associated nuclear DNA damage and repair in the sleeping chironomid: Linkage with radioresistance. PLoS One 5: e14008. doi: 10.1371/journal.pone.0014008 PMID: 21103355

12. Alpert P (2005) The limits and frontiers of desiccation-tolerant life. Integr Comp Biol 45: 685-695. doi: 10.1093/icb/45.5.685 PMID: 21676818

13. Potts M, Slaughter SM, Hunneke F-U, Garst JF, Helm RF (2005) Desiccation tolerance of prokaryotes: application of principles to human cells. Integr Comp Biol 45: 800-809. doi: 10.1093/icb/45.5.800 PMID: 21676831

14. Elbein AD, Pan YT, Pastuszak I, Carroll D (2003) New insights on trehalose: a multifunctional molecule. Glycobiology 13: 17R-27R. doi: 10.1093/glycob/cwg047 PMID: 12626396

15. Fernandez O, Béthencourt L, Quero A, Sangwan RS, Clément C (2010) Trehalose and plant stress responses: friend or foe? Trends Plant Sci 15: 409-417. doi: 10.1016/j.tplants.2010.04.004 PMID: 20494608

16. Ohtake S, Wang Y (2011) Trehalose: current use and future applications. J Pharm Sci 100: 2020-2053. doi: 10.1002/jps PMID: 21337544

17. Potts M (2001) Desiccation tolerance: a simple process? Trends Microbiol 9: 553-559. PMID: 11825716

18. Jain NK, Roy I (2009) Effect of trehalose on protein structure. Protein Sci 18: 24-36. doi: 10.1002/pro.3 PMID: 19177348

19. Zaparty M, Hagemann A, Bräsen C, Hensel R, Lupas AN, Brinkmann H, et al. (2013) The first prokaryotic trehalose synthase complex identified in the hyperthermophilic crenarchaeon Thermoproteus tenax. PLoS One 8: e61354. doi: 10.1371/journal.pone.0061354 PMID: 23626675

20. Pampurova S, Verschooten K, Avonce N, Van Dijck P (2014) Functional screening of a cDNA library from the desiccation-tolerant plant Selaginella lepidophylla in yeast mutants identifies trehalose biosynthesis genes of plant and microbial origin. J Plant Res 127: 803-813. doi: 10.1007/s10265-014-0663-x PMID: 25246071

21. Sakurai M, Furuki T, Akao K-I, Tanaka D, Nakahara Y, Kikawada T, et al. (2008) Vitrification is essential for anhydrobiosis in an African chironomid, Polypedilum vanderplanki. Proc Natl Acad Sci U S A 105: 5093-5098. doi: 10.1073/pnas.0706197105 PMID: 18362351

22. Madin KAC, Crowe JH (1975) Anhydrobiosis in nematodes: Carbohydrate and lipid metabolism during dehydration. J Exp Zool 193: 335-342. doi: 10.1002/jez.1401930309
23. Adhikari BN, Wall DH, Adams BJ (2009) Desiccation survival in an Antarctic nematode: molecular analysis using expressed sequenced tags. BMC Genomics 10: 69. doi: 10.1186/1471-2164-10-69 PMID: 19203352

24. Hengherr S, Heyer AG, Köhler H-R, Schill RO (2008) Trehalose and anhydrobiosis in tardigrades—evidence for divergence in responses to dehydration. FEBS J 275: 281–2888. doi: 10.1111/j.1742-4658. 2007.06198.x PMID: 18070104

25. Jönsson K, Persson O (2010) Trehalose in three species of desiccation tolerant tardigrades. Open Zool J: 1–5.

26. Goddijn O, Dun K Van (1999) Trehalose metabolism in plants. Trends Plant Sci 1385: 315–319.

27. Van Houtte H, Vandesteene L, López-Galvis L, Lemmens L, Kissel E, Carpentier S, et al. (2013) Over-expression of the trehalase gene AtTRE1 leads to increased drought stress tolerance in Arabidopsis and is involved in abscisic acid-induced stomatal closure. Plant Physiol 161: 1158–1171. doi: 10.1104/ pp.112.211391 PMID: 23341362

28. Delorge I, Janiak M, Carpentier S, Van Dijck P (2014) Fine tuning of trehalose biosynthesis and hydro-lysis as novel tools for the generation of abiotic stress tolerant plants. Front Plant Sci 5: 147. doi: 10. 3389/fpls.2014.00147 PMID: 24782885

29. Li H-W, Zang B-S, Deng X-W, Wang X-P (2011) Overexpression of the trehalose-6-phosphate synthase gene OsTPS1 enhances abiotic stress tolerance in rice. Planta 234: 1007–1018. doi: 10. 1007/s00425-011-1458-0 PMID: 21698458

30. Ponnu J, Wahl V, Schmid M (2011) Trehalose-6-phosphate: connecting plant metabolism and development. Front Plant Sci 2. doi: 10.3389/fpls.2011.00070

31. Lunn JE, Delorge I, Figueroa CM, Van Dijck P, Stitt M (2014) Trehalose metabolism in plants. Plant J: 544–567. doi: 10.1111/tpj.12509 PMID: 24645920

32. Tsai AY- L, Gazzarrini S (2014) Trehalose-6-phosphate and SnRK1 kinases in plant development and signaling: the emerging picture. Front Plant Sci 5: 119. doi: 10.3389/fpls.2014.00119 PMID: 24744765

33. Avonce N, Mendoza-Vargas A, Morett E, Iturriaga G (2006) Insights on the evolution of trehalose bio-synthesis. BMC Evol Biol 6: 109. doi: 10.1186/1471-2148-6-109 PMID: 17178000

34. Lapinski J, Tunnacliffe A (2003) Anhydrobiosis without trehalose in bdelloid rotifers. FEBS Lett 553: 387–390. PMID: 14572656

35. Tunnacliffe A, Lapinski J (2003) Resurrecting Van Leeuwenhoek's rotifers: a reappraisal of the role of disaccharides in anhydrobiosis. Philos Trans R Soc London Ser B Biol Sci 358: 1755–1771.

36. Caprioli M, Katholm AK, Melone G, Ramlov H, Ricci C, Santo N (2004) Trehalose in desiccated rotifers: a comparison between a bdelloid and a monogonont species. Comp Biochem Physiol A Mol Integr Physiol 139: 527–532. PMID: 15596399

37. Boschetti C (2011) Foreign genes and novel hydrophilic protein genes participate in the desiccation response of the bdelloid rotifer *Adineta ricciae*. J Exp Biol 214: 59–68. doi: 10.1242/jeb.050328 PMID: 21147969

38. Suga K, Mark Welch D, Tanaka Y, Sakakura Y, Hagiwara A (2007) Analysis of expressed sequence tags of the cyclically parthenogenetic rotifer Brachionus plicatilis. PLoS One 2: e671. doi: 10.1371/ journal.pone.0000671 PMID: 17668053

39. Flot J-F, Hespeels B, Li X, Noel B, Arkhipova I, Danchin EGJ, et al. (2013) Genomic evidence for ameiotic evolution in the bdelloid rotifer *Adineta vaga*. Nature 500: 453–457. doi: 10.1038/nature12326 PMID: 23873043

40. Gertz EM, Yu Y-K, Agarwala R, Schäffer A a, Altschul SF (2006) Composition-based statistics and translated nucleotide searches: improving the TBLASTN module of BLAST. BMC Biol 4: 41. doi: 10. 1186/1741-7007-4-41 PMID: 17156431

41. Hall T (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41: 95–98.

42. Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, et al. (2012) MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res 40: e49. doi: 10.1093/nar/ gkr1293 PMID: 22217600

43. Birney E, Clamp M, Durbin R (2004) GeneWise and Genomewise. Genome Res 14: 988–995. doi: 10. 1101/gr.1865504 PMID: 15123596

44. Friesen VL (2000) Introns. In: Baker Allan, editor. Molecular Methods in Ecology. Wiley-Blackwell. pp. 274–290.

45. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402. doi: 10.1093/nar/25.17.3389 PMID: 9254694
23. Adhikari BN, Wall DH, Adams BJ (2009) Desiccation survival in an Antarctic nematode: molecular analysis using expressed sequenced tags. BMC Genomics 10: 69. doi: 10.1186/1471-2164-10-69 PMID: 19203352

24. Hengherr S, Heyer AG, Köhler H-R, Schill RO (2008) Trehalose and anhydrobiosis in tardigrades—evidence for divergence in responses to dehydration. FEBS J 275: 281–2888. doi: 10.1111/j.1742-4658. 2007.06198.x PMID: 18070104

25. Jönsson K, Persson O (2010) Trehalose in three species of desiccation tolerant tardigrades. Open Zool J: 1–5.

26. Goddijn O, Dun K Van (1999) Trehalose metabolism in plants. Trends Plant Sci 1385: 315–319.

27. Van Houtte H, Vandesteene L, López-Galvis L, Lemmens L, Kissel E, Carpentier S, et al. (2013) Over-expression of the trehalase gene AtTRE1 leads to increased drought stress tolerance in Arabidopsis and is involved in abscisic acid-induced stomatal closure. Plant Physiol 161: 1158–1171. doi: 10.1104/ pp.112.211391 PMID: 23341362

28. Delorge I, Janiak M, Carpentier S, Van Dijck P (2014) Fine tuning of trehalose biosynthesis and hydro-lysis as novel tools for the generation of abiotic stress tolerant plants. Front Plant Sci 5: 147. doi: 10. 3389/fpls.2014.00147 PMID: 24782885

29. Li H-W, Zang B-S, Deng X-W, Wang X-P (2011) Overexpression of the trehalose-6-phosphate synthase gene OsTPS1 enhances abiotic stress tolerance in rice. Planta 234: 1007–1018. doi: 10. 1007/s00425-011-1458-0 PMID: 21698458

30. Ponnu J, Wahl V, Schmid M (2011) Trehalose-6-phosphate: connecting plant metabolism and development. Front Plant Sci 2. doi: 10.3389/fpls.2011.00070

31. Lunn JE, Delorge I, Figueroa CM, Van Dijck P, Stitt M (2014) Trehalose metabolism in plants. Plant J: 544–567. doi: 10.1111/tpj.12509 PMID: 24645920

32. Tsai AY- L, Gazzarrini S (2014) Trehalose-6-phosphate and SnRK1 kinases in plant development and signaling: the emerging picture. Front Plant Sci 5: 119. doi: 10.3389/fpls.2014.00119 PMID: 24744765

33. Avonce N, Mendoza-Vargas A, Morett E, Iturriaga G (2006) Insights on the evolution of trehalose bio-synthesis. BMC Evol Biol 6: 109. doi: 10.1186/1471-2148-6-109 PMID: 17178000

34. Lapinski J, Tunnacliffe A (2003) Anhydrobiosis without trehalose in bdelloid rotifers. FEBS Lett 553: 387–390. PMID: 14572656

35. Tunnacliffe A, Lapinski J (2003) Resurrecting Van Leeuwenhoek's rotifers: a reappraisal of the role of disaccharides in anhydrobiosis. Philos Trans R Soc London Ser B Biol Sci 358: 1755–1771.

36. Caprioli M, Katholm AK, Melone G, Ramlov H, Ricci C, Santo N (2004) Trehalose in desiccated rotifers: a comparison between a bdelloid and a monogonont species. Comp Biochem Physiol A Mol Integr Physiol 139: 527–532. PMID: 15596399

37. Boschetti C (2011) Foreign genes and novel hydrophilic protein genes participate in the desiccation response of the bdelloid rotifer *Adineta ricciae*. J Exp Biol 214: 59–68. doi: 10.1242/jeb.050328 PMID: 21147969

38. Suga K, Mark Welch D, Tanaka Y, Sakakura Y, Hagiwara A (2007) Analysis of expressed sequence tags of the cyclically parthenogenetic rotifer Brachionus plicatilis. PLoS One 2: e671. doi: 10.1371/ journal.pone.0000671 PMID: 17668053

39. Flot J-F, Hespeels B, Li X, Noel B, Arkhipova I, Danchin EGJ, et al. (2013) Genomic evidence for ameiotic evolution in the bdelloid rotifer *Adineta vaga*. Nature 500: 453–457. doi: 10.1038/nature12326 PMID: 23873043

40. Gertz EM, Yu Y-K, Agarwala R, Schäffer A a, Altschul SF (2006) Composition-based statistics and translated nucleotide searches: improving the TBLASTN module of BLAST. BMC Biol 4: 41. doi: 10. 1186/1741-7007-4-41 PMID: 17156431

41. Hall T (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41: 95–98.

42. Wang Y, Tang H, Debarry JD, Tan X, Li J, Wang X, et al. (2012) MCScanX: A toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res 40: e49. doi: 10.1093/nar/ gkr1293 PMID: 22217600

43. Birney E, Clamp M, Durbin R (2004) GeneWise and Genomewise. Genome Res 14: 988–995. doi: 10. 1101/gr.1865504 PMID: 15123596

44. Friesen VL (2000) Introns. In: Baker Allan, editor. Molecular Methods in Ecology. Wiley-Blackwell. pp. 274–290.

45. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389–3402. doi: 10.1093/nar/25.17.3389 PMID: 9254694
46. Gladyshev E a, Meselson M, Arkhipova IR (2008) Massive horizontal gene transfer in bdelloid rotifers. Science 320: 1210–1213. doi: 10.1126/science.1156407 PMID: 18511688

47. Hur JH, Van Doninck K, Mandigo ML, Meselson M (2009) Degenerate tetraploidy was established before bdelloid rotifer families diverged. Mol Biol Evol 26: 375–383. doi: 10.1093/molbev/msn260 PMID: 18996928

48. Gibson R, Turkenburg J, Charnock S (2002) Insights into trehalose synthesis provided by the structure of the retaining glucosyltransferase OtsA. Chem Biol 9: 1337–1346. PMID: 12498887

49. Gibson RP, Tarling C a, Roberts S, Withers SG, Davies GJ (2004) The donor subsite of trehalose-6-phosphate synthase: binary complexes with UDP-glucose and UDP-2-deoxy-2-fluoro-glucose at 2 Å resolution. J Biol Chem 279: 1950–1955. doi: 10.1074/jbc.M307643200 PMID: 14570926

50. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797. doi: 10.1093/nar/gkh340

51. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 30: 2725–2729. doi: 10.1093/molbev/mst197 PMID: 24132122

52. Schwede T (2003) SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res 31: 3381–3385. doi: 10.1093/nar/gkg520 PMID: 12824332

53. Arnold K, Bordoli L, Kopp J, Schwede T (2006) The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics 22: 195–201. doi: 10.1093/bioinformatics/bti770 PMID: 16301204

54. Biasini M, Bienert S, Waterhouse A, Arnold K, Studer G, Schmidt T, et al. (2014) SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res 42: W252–W258. doi: 10.1093/nar/gku340 PMID: 24782522

55. Kelley LA, Sternberg MJ (2009) Protein structure prediction on the Web: a case study using the Phyre server. Nat Protoc 4: 363–371. doi: 10.1038/nprot.2009.2 PMID: 19247286

56. Darriba D, Taboada GL, Doallo R, Posada D (2011) ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics 27: 1164–1165. doi: 10.1093/bioinformatics/btr088 PMID: 21335321

57. Schwarz G (1978) Estimating the dimension of a model. Ann Stat 6: 461–464.

58. Akaike H (1974) A new look at the statistical model identification. Autom Control IEEE Trans 19: 716–723.

59. Stamatakis A (2006) RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690. doi: 10.1093/bioinformatics/btl446 PMID: 16928733

60. Silvestro D, Michalak I (2011) raxmlGUI: a graphical front-end for RAxML. Org Divers Evol 12: 335–337. doi: 10.1007/s13127-011-0056-0

61. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574. doi: 10.1093/bioinformatics/btg180 PMID: 12912839

62. Milne I, Lindner D, Bayer M, Husmeier D, McGuire G, Marshall DF, et al. (2009) TOPALi v2: a rich graphical interface for evolutionary analyses of multiple alignments on HPC clusters and multi-core desktops. Bioinformatics 25: 126–127. doi: 10.1093/bioinformatics/btn575 PMID: 18984599

63. Helaers R, Milinkovitch MC (2010) MetaPIGA v2.0: maximum likelihood large phylogeny estimation using the metapopulation genetic algorithm and other stochastic heuristics. BMC Bioinformatics 11: 379. doi: 10.1186/1471-2105-11-379 PMID: 20633263

64. Hespeels B, Knapen M, Hanot-Mambres D, Heuskin A- C, Pineux F, Lucas S, et al. (2014) Gateway to genetic exchange? DNA double-strand breaks in the bdelloid rotifer *Adineta vaga* submitted to desiccation. J Evol Biol 27: 1334–1345. doi: 10.1111/jeb.12326 PMID: 25105197

65. Howe KL, Chothia T, Durbin R (2002) GAZE: a generic framework for the integration of gene-prediction data by dynamic programming. Genome Res 12: 1418–1427. doi: 10.1101/gr.149502 PMID: 12213779

66. Bustin S a, Benes V, Garson J a, Hellemans J, Huggett J, Kubista M, et al. (2009) The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem 55: 611–622. doi: 10.1373/clinchem.2008.112797 PMID: 19246619

67. Van Doninck K, Mandigo ML, Hur JH, Wang P, Guglielmini J, Milinkovitch MC, et al. (2009) Phylogenomics of unusual histone H2A variants in bdelloid rotifers. PLoS Genet 5: e1000401. doi: 10.1371/journal.pgen.1000401 PMID: 19266019

68. Rao KN, Kumaran D, Seetharaman J, Bonanno JB, Burley SK, Swaminathan S (2006) Crystal structure of trehalose-6-phosphate phosphatase-related protein: Biochemical and biological implications: 1735–1744. doi: 10.1110/ps.062096606.ily
46. Gladyshev E a, Meselson M, Arkhipova IR (2008) Massive horizontal gene transfer in bdelloid rotifers. Science 320: 1210–1213. doi: 10.1126/science.1156407 PMID: 18511688

47. Hur JH, Van Doninck K, Mandigo ML, Meselson M (2009) Degenerate tetraploidy was established before bdelloid rotifer families diverged. Mol Biol Evol 26: 375–383. doi: 10.1093/molbev/msn260 PMID: 18996928

48. Gibson R, Turkenburg J, Charnock S (2002) Insights into trehalose synthesis provided by the structure of the retaining glucosyltransferase OtsA. Chem Biol 9: 1337–1346. PMID: 12498887

49. Gibson RP, Tarling C a, Roberts S, Withers SG, Davies GJ (2004) The donor subsite of trehalose-6-phosphate synthase: binary complexes with UDP-glucose and UDP-2-deoxy-2-fluoro-glucose at 2 Å resolution. J Biol Chem 279: 1950–1955. doi: 10.1074/jbc.M307643200 PMID: 14570926

50. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797. doi: 10.1093/nar/gkh340

51. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S (2013) MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol 30: 2725–2729. doi: 10.1093/molbev/mst197 PMID: 24132122

52. Schwede T (2003) SWISS-MODEL: an automated protein homology-modeling server. Nucleic Acids Res 31: 3381–3385. doi: 10.1093/nar/gkg520 PMID: 12824332

53. Arnold K, Bordoli L, Kopp J, Schwede T (2006) The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics 22: 195–201. doi: 10.1093/bioinformatics/bti770 PMID: 16301204

54. Biasini M, Bienert S, Waterhouse A, Arnold K, Studer G, Schmidt T, et al. (2014) SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res 42: W252–W258. doi: 10.1093/nar/gku340 PMID: 24782522

55. Kelley LA, Sternberg MJ (2009) Protein structure prediction on the Web: a case study using the Phyre server. Nat Protoc 4: 363–371. doi: 10.1038/nprot.2009.2 PMID: 19247286

56. Darriba D, Taboada GL, Doallo R, Posada D (2011) ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics 27: 1164–1165. doi: 10.1093/bioinformatics/btr088 PMID: 21335321

57. Schwarz G (1978) Estimating the dimension of a model. Ann Stat 6: 461–464.

58. Akaike H (1974) A new look at the statistical model identification. Autom Control IEEE Trans 19: 716–723.

59. Stamatakis A (2006) RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690. doi: 10.1093/bioinformatics/btl446 PMID: 16928733

60. Silvestro D, Michalak I (2011) raxmlGUI: a graphical front-end for RAxML. Org Divers Evol 12: 335–337. doi: 10.1007/s13127-011-0056-0

61. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574. doi: 10.1093/bioinformatics/btg180 PMID: 12912839

62. Milne I, Lindner D, Bayer M, Husmeier D, McGuire G, Marshall DF, et al. (2009) TOPALi v2: a rich graphical interface for evolutionary analyses of multiple alignments on HPC clusters and multi-core desktops. Bioinformatics 25: 126–127. doi: 10.1093/bioinformatics/btn575 PMID: 18984599

63. Helaers R, Milinkovitch MC (2010) MetaPIGA v2.0: maximum likelihood large phylogeny estimation using the metapopulation genetic algorithm and other stochastic heuristics. BMC Bioinformatics 11: 379. doi: 10.1186/1471-2105-11-379 PMID: 20633263

64. Hespeels B, Knapen M, Hanot-Mambres D, Heuskin A- C, Pineux F, Lucas S, et al. (2014) Gateway to genetic exchange? DNA double-strand breaks in the bdelloid rotifer *Adineta vaga* submitted to desiccation. J Evol Biol 27: 1334–1345. doi: 10.1111/jeb.12326 PMID: 25105197

65. Howe KL, Chothia T, Durbin R (2002) GAZE: a generic framework for the integration of gene-prediction data by dynamic programming. Genome Res 12: 1418–1427. doi: 10.1101/gr.149502 PMID: 12213779

66. Bustin S a, Benes V, Garson J a, Hellemans J, Huggett J, Kubista M, et al. (2009) The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem 55: 611–622. doi: 10.1373/clinchem.2008.112797 PMID: 19246619

67. Van Doninck K, Mandigo ML, Hur JH, Wang P, Guglielmini J, Milinkovitch MC, et al. (2009) Phylogenomics of unusual histone H2A variants in bdelloid rotifers. PLoS Genet 5: e1000401. doi: 10.1371/journal.pgen.1000401 PMID: 19266019

68. Rao KN, Kumaran D, Seetharaman J, Bonanno JB, Burley SK, Swaminathan S (2006) Crystal structure of trehalose-6-phosphate phosphatase-related protein: Biochemical and biological implications: 1735–1744. doi: 10.1110/ps.062096606.ily
69. Vandesteene L, Ramon M, Le Roy K, Van Dijck P, Rolland F (2010) A single active trehalose-6-P synthase (TPS) and a family of putative regulatory TPS-like proteins in Arabidopsis. Mol Plant 3: 406–419. doi: 10.1093/mp/ssp114 PMID: 20100798

70. Vandesteene L, López-Galvis L, Vanneste K, Feil R, Maere S, Lammens W, et al. (2012) Expansive evolution of the trehalose-6-phosphate phosphatase gene family in Arabidopsis. Plant Physiol 160: 884–896. doi: 10.1104/pp.112.201400 PMID: 22855938

71. Avonce N, Wuyts J, Verschooten K, Vandesteene L, Van Dijck P (2010) The Cytophaga hutchinsonii ChTPSP: First characterized bifunctional TPS-TPP protein as putative ancestor of all eukaryotic trehalose biosynthesis proteins. Mol Biol Evol 27: 359–369. doi: 10.1093/molbev/msp241 PMID: 19812028

72. Gladyshev E, Meselson M (2008) Extreme resistance of bdelloid rotifers to ionizing radiation. Proc Natl Acad Sci U S A 105: 5139–5144. doi: 10.1073/pnas.0800966105 PMID: 18362355

73. Gladyshev E, Arkhipova I (2010) Genome structure of bdelloid rotifers: shaped by asexuality or desiccation? J Hered 101: 85–93. doi: 10.1093/jhered/esq008

74. Müller J, Aeschbacher R a, Wingler A, Boller T, Wiemken A (2001) Trehalose and trehalase in Arabidopsis. Plant Physiol 125: 1086–1093. PMID: 11161063

75. Goddijn OJ, Verwoerd TC, Voogd E, Krutwagen RW, De Graff PTHM, Poels J, ... & Pen J (1997) Inhibition of trehalase activity enhances trehalose accumulation in transgenic plants. Plant Physiol 113: 181–190. PMID: 9008394

76. Dijck P Van, Rop L De, Szlufcik K, Van Ael E, Thevelein JM (2002) Disruption of the Candida albicans TPS2 gene encoding trehalose-6-phosphate phosphatase decreases infectivity without affecting hypha formation. Infect Immun 70: 1772–1782. doi: 10.1128/IAI.70.4.1772 PMID: 11895938

77. Puttikamonkul S, Willger SD, Grahl N, Perfect JR, Movahed N, Bothner B, et al. (2010) Trehalose 6-phosphate phosphatase is required for cell wall integrity and fungal virulence but not trehalose biosynthesis in the human fungal pathogen Aspergillus fumigatus. Mol Microbiol 77: 891–911. doi: 10.1111/j.1365-2958.2010.07254.x PMID: 20545865

78. Kaasen I, Falkenberg P, Styrvold OB, Strøm A R (1992) Molecular cloning and physical mapping of the otsBA genes, which encode the osmoregulatory trehalose pathway of Escherichia coli: evidence that transcription is activated by katF (AppR). J Bacteriol 174: 889–898. PMID: 1310094

79. Wingler A (2002) The function of trehalose biosynthesis in plants. Phytochemistry 60: 437–440. PMID: 12052507

80. Schluepmann H, Pellny T, van Dijken A, Smeekens S, Paul M (2003) Trehalose 6-phosphate is indispensable for carbohydrate utilization and growth in Arabidopsis thaliana. Proc Natl Acad Sci U S A 100: 6849–6854. doi: 10.1073/pnas.1132018100 PMID: 12748379

81. Satoh-Nagasawa N, Nagasawa N, Malcomber S, Sakai H, Jackson D (2006) Trehalose metabolism and signaling. Nature 441: 227–230. doi: 10.1038/nature04725 PMID: 16688177

82. Yadav UP, Ivakov A, Feil R, Duan GY, Walther D, Giavalisco P, et al. (2014) The sucrose-trehalose 6-phosphate (Tre6P) nexus: specificity and mechanisms of sucrose signalling by Tre6P. J Exp Bot 65: 1051–1068. doi: 10.1093/jxb/ert457 PMID: 24420566

83. Blázquez MA, Santos E, Flores CL, Martínez-Zapater JM, Salinas J, Gancedo C (1998) Isolation and molecular characterization of the Arabidopsis TPS1 gene, encoding trehalose-6-phosphate synthase. Plant J 13: 685–689. PMID: 9681010

84. Avonce N, Leyman B (2004) The Arabidopsis trehalose-6-P synthase AtTPS1 gene is a regulator of glucose, abscisic acid, and stress signaling. Plant Physiol 136: 3649–3659. doi: 10.1104/pp.104.052084.1 PMID: 15516499

85. Lunn JE, Feil R, Hendriks JHM, Gibon Y, Morcuende R, Osuna D, et al. (2006) Sugar-induced increases in trehalose 6-phosphate are correlated with redox activation of ADP-glucose pyrophosphorylase and higher rates of starch synthesis in Arabidopsis thaliana. Biochem J 397: 139–148. doi: 10.1042/BJ20060083 PMID: 16551270

86. Delorge I, Figueroa CM, Feil R, Lunn JE, Van Dijck P (2014) TREHALOSE-6-PHOSPHATE SYNTHASE 1 is not the only active TPS in Arabidopsis thaliana. Biochem J. doi: 10.1042/BJ20141322

87. Fernandez J, Wilson R (2011) The sugar sensor, trehalose-6-phosphate synthase (Tps1), regulates primary and secondary metabolism during infection by the rice blast fungus: Will Magnaporthe oryzae’s "sweet tooth" become its "Achilles’ heel"? Mycology 2: 46–53. doi: 10.1080/21501203.2011.563431

88. De Maesschalck C (2009) Unraveling the function of trehalose biosynthesis genes in lateral root development University of Ghent, Belgium.
69. Vandesteene L, Ramon M, Le Roy K, Van Dijck P, Rolland F (2010) A single active trehalose-6-P synthase (TPS) and a family of putative regulatory TPS-like proteins in Arabidopsis. Mol Plant 3: 406–419. doi: 10.1093/mp/ssp114 PMID: 20100798

70. Vandesteene L, López-Galvis L, Vanneste K, Feil R, Maere S, Lammens W, et al. (2012) Expansive evolution of the trehalose-6-phosphate phosphatase gene family in Arabidopsis. Plant Physiol 160: 884–896. doi: 10.1104/pp.112.201400 PMID: 22855938

71. Avonce N, Wuyts J, Verschooten K, Vandesteene L, Van Dijck P (2010) The Cytophaga hutchinsonii ChTPSP: First characterized bifunctional TPS-TPP protein as putative ancestor of all eukaryotic trehalose biosynthesis proteins. Mol Biol Evol 27: 359–369. doi: 10.1093/molbev/msp241 PMID: 19812028

72. Gladyshev E, Meselson M (2008) Extreme resistance of bdelloid rotifers to ionizing radiation. Proc Natl Acad Sci U S A 105: 5139–5144. doi: 10.1073/pnas.0800966105 PMID: 18362355

73. Gladyshev E, Arkhipova I (2010) Genome structure of bdelloid rotifers: shaped by asexuality or desiccation? J Hered 101: 85–93. doi: 10.1093/jhered/esq008

74. Müller J, Aeschbacher R a, Wingler A, Boller T, Wiemken A (2001) Trehalose and trehalase in Arabidopsis. Plant Physiol 125: 1086–1093. PMID: 11161063

75. Goddijn OJ, Verwoerd TC, Voogd E, Krutwagen RW, De Graff PTHM, Poels J, ... & Pen J (1997) Inhibition of trehalase activity enhances trehalose accumulation in transgenic plants. Plant Physiol 113: 181–190. PMID: 9008394

76. Dijck P Van, Rop L De, Szlufcik K, Van Ael E, Thevelein JM (2002) Disruption of the Candida albicans TPS2 gene encoding trehalose-6-phosphate phosphatase decreases infectivity without affecting hypha formation. Infect Immun 70: 1772–1782. doi: 10.1128/IAI.70.4.1772 PMID: 11895938

77. Puttikamonkul S, Willger SD, Grah N, Perfect JR, Movahed N, Bothner B, et al. (2010) Trehalose 6-phosphate phosphatase is required for cell wall integrity and fungal virulence but not trehalose biosynthesis in the human fungal pathogen Aspergillus fumigatus. Mol Microbiol 77: 891–911. doi: 10.1111/j.1365-2958.2010.07254.x PMID: 20545865

78. Kaasen I, Falkenberg P, Styrvold OB, Strøm A R (1992) Molecular cloning and physical mapping of the otsBA genes, which encode the osmoregulatory trehalose pathway of Escherichia coli: evidence that transcription is activated by katF (AppR). J Bacteriol 174: 889–898. PMID: 1310094

79. Wingler A (2002) The function of trehalose biosynthesis in plants. Phytochemistry 60: 437–440. PMID: 12052507

80. Schluepmann H, Pellny T, van Dijken A, Smeekens S, Paul M (2003) Trehalose 6-phosphate is indispensable for carbohydrate utilization and growth in Arabidopsis thaliana. Proc Natl Acad Sci U S A 100: 6849–6854. doi: 10.1073/pnas.1132018100 PMID: 12748379

81. Satoh-Nagasawa N, Nagasawa N, Malcomber S, Sakai H, Jackson D (2006) Trehalose metabolism and signaling. Nature 441: 227–230. doi: 10.1038/nature04725 PMID: 16688177

82. Yadav UP, Ivakov A, Feil R, Duan GY, Walther D, Giavalisco P, et al. (2014) The sucrose-trehalose 6-phosphate (Tre6P) nexus: specificity and mechanisms of sucrose signalling by Tre6P. J Exp Bot 65: 1051–1068. doi: 10.1093/jxb/ert457 PMID: 24420566

83. Blázquez MA, Santos E, Flores CL, Martínez-Zapater JM, Salinas J, Gancedo C (1998) Isolation and molecular characterization of the Arabidopsis TPS1 gene, encoding trehalose-6-phosphate synthase. Plant J 13: 685–689. PMID: 9681010

84. Avonce N, Leyman B (2004) The Arabidopsis trehalose-6-P synthase AtTPS1 gene is a regulator of glucose, abscisic acid, and stress signaling. Plant Physiol 136: 3649–3659. doi: 10.1104/pp.104.052084.1 PMID: 15516499

85. Lunn JE, Feil R, Hendriks JHM, Gibon Y, Morcuende R, Osuna D, et al. (2006) Sugar-induced increases in trehalose 6-phosphate are correlated with redox activation of ADPglucose pyrophosphorylase and higher rates of starch synthesis in Arabidopsis thaliana. Biochem J 397: 139–148. doi: 10.1042/BJ20060083 PMID: 16551270

86. Delorge I, Figueroa CM, Feil R, Lunn JE, Van Dijck P (2014) TREHALOSE-6-PHOSPHATE SYNTHASE 1 is not the only active TPS in Arabidopsis thaliana. Biochem J. doi: 10.1042/BJ20141322

87. Fernandez J, Wilson R (2011) The sugar sensor, trehalose-6-phosphate synthase (Tps1), regulates primary and secondary metabolism during infection by the rice blast fungus: Will Magnaporthe oryzae's "sweet tooth" become its "Achilles' heel"? Mycology 2: 46–53. doi: 10.1080/21501203.2011.563431

88. De Maesschalck C (2009) Unraveling the function of trehalose biosynthesis genes in lateral root development University of Ghent, Belgium.