Caenorhabditis monodelphis sp. n.: defining the stem morphology and genomics of the genus Caenorhabditis
© The Author(s) 2017
Received: 20 October 2016
Accepted: 18 January 2017
Published: 23 February 2017
The genus Caenorhabditis has been central to our understanding of metazoan biology. The best-known species, Caenorhabditis elegans, is but one member of a genus with around 50 known species, and knowledge of these species will place the singular example of C. elegans in a rich phylogenetic context. How did the model come to be as it is today, and what are the dynamics of change in the genus?
As part of this effort to “put C. elegans in its place”, we here describe the morphology and genome of Caenorhabditis monodelphis sp. n., previously known as Caenorhabditis sp. 1. Like many other Caenorhabditis, C. monodelphis sp. n. has a phoretic association with a transport host, in this case with the fungivorous beetle Cis castaneus. Using genomic data, we place C. monodelphis sp. n. as sister to all other Caenorhabditis for which genome data are available. Using this genome phylogeny, we reconstruct the stemspecies morphological pattern of Caenorhabditis.
With the morphological and genomic description of C. monodelphis sp. n., another key species for evolutionary and developmental studies within Caenorhabditis becomes available. The most important characters are its early diverging position, unique morphology for the genus and its similarities with the hypothetical ancestor of Caenorhabditis.
KeywordsTaxonomy Systematics Evolution Genome Phylogeny Description
The nematode genus Caenorhabditis includes the well-known model organism C. elegans, which has provided key insights into molecular and developmental biology . Over the past ten years, numerous new Caenorhabditis species have been discovered and described [2, 3]. These putative new taxa are generally indistinguishable morphologically, and thus the most recent descriptions of new species within Caenorhabditis have been based on DNA sequences and mating tests only . This streamlined species-description methodology has been driven by the need to have names to attach to real biological entities, and the fact that traditional taxonomy has been unable to keep up with species discovery. The method is relatively simple to implement, and delivers taxa that have a biological reality . However, as the number of species discovered in Caenorhabditis grows, traditional, morphological descriptions are still valuable for the understanding of patterns of trait evolution and inference of ecological functions [4, 5]. Although morphology cannot be used to definitively delineate species, it should not be abandoned all together.
M-A Félix, C Braendle and AD Cutter  provided new species name designations for 15 biological species, considerably increasing the number of named Caenorhabditis species in laboratory culture. However, several key Caenorhabditis species remain undescribed. A well-known but undescribed species of Caenorhabditis, informally referred to as Caenorhabditis sp. 1, has been analysed in several evolutionary and developmental studies [3, 6–8]. C. sp. 1 was previously found only once inside a fruiting body of the fungus Ganoderma applanatum (Pers.) Pat. (Polyporaceae), growing on the stump of tree in Berlin, Germany. Galleries inside the fungus were frequently visited by beetles of the species Cis castaneus (Ciidae), a beetle with a host preference for Ganoderma . Associations between nematodes and insects, where the nematode uses the insect as a transport carrier (phoretism), have already been described for several Caenorhabditis species, including Caenorhabditis angaria, C. remanei, and C. bovis, and similar phoretic associations could be expected for many or possibly all other Caenorhabditis species .
Here we use both morphological and molecular analyses to characterise and describe C. sp. 1 as a new species, Caenorhabditis monodelphis sp. n., and explore its relationship with the beetle Cis castaneus. Molecular phylogenetic analysis based on whole genome sequencing of an inbred derivative of the type strain affirms the placement of C. monodelphis sp. n. as sister to other analysed Caenorhabditis, and we analyse the evolution of phenotypic traits to infer those present in the hypothetical ancestor of Caenorhabditis.
Isolation and culture
Caenorhabditis monodelphis sp. n. (strain SB341) was originally isolated from fruiting bodies of Ganoderma applanatum (Pers.) Pat. 1887 collected in Berlin-Grunewald, Germany (April, 2001) and later from four locations in Belgium (strain DSC001 collected from 51°06'24"N, 3°18'13"E, March 2014, strain DSC002 collected from 50°52'7"N, 4°06'54", February 2014, and an uncultured population 51°02'41"N, 3°27'17" June 2014) and from one location in the Botanical Garden in Oslo, Norway (strain JU2884; 59°55'04"N 10°46'01"E, 22 July 2015). These collections were from the same mushroom species. We also found C. monodelphis sp. n. in the fruiting body of Fomes fomentarius (L.) Fr. 1849 (50°43'02"N, 4°05'06"E, February 2015). Strain SB341 was chosen as type.
Nematodes were extracted from the fruiting bodies of G. applanatum using the modified Baermann method . Dauer larvae were isolated from the beetle Cis castaneus (Herbst, 1793) that had been extracted from the same mushroom from multiple locations (except the type population and from 51°02'41"N, 3°27'17"). Adults and dauer larvae were picked out and cultured on nutrient agar plates seeded with E. coli OP50 at 15 °C.
Cultures of nematodes from two populations (strain SB341 and DSC001) were used for the description. Measurements and drawings were made with an Olympus BX51 equipped with differential interference contrast (DIC). Light microscopic images were taken with a Nikon DS-FI2 camera. For Scanning Electron Microscopy (SEM), two fixation methods were used. For the first fixation method, live animals were fixed in a microwave in Trump’s fixative (2% paraformaldehyde + 2.5% glutaraldehyde in a 0.1 M Sorenson buffer) for a few seconds. Specimens were subsequently washed three times in double-distilled water. For the second method, specimens were put in a refrigerator at 4 °C for 1 h, then Trump’s fixative was added and specimens were left overnight at 4 °C. The specimens were then washed with a 0.2 M phosphate buffer followed by 1 h post-fixation in a 1% OsO4 solution at room temperature and subsequently washed 4 times in double-distilled water. For both methods, the specimens were dehydrated by passing them through a graded ethanol concentration series of 30, 50, 75, 95% (20 min each) and 3x 100% (10 min each). The specimens were critical point-dried with liquid CO2, mounted on stubs with carbon discs and coated with gold (25 nm) before observation with a JSM-840 EM (JEOL, Tokyo, Japan) at 15 kV. Sperm cells were observed in the female post-uterine sac with Transmission Electron Microscopy (TEM), processing samples as described , except for ultramicrotomy with a Leica EM UC7 and 1 h 1% osmium postfixation (Slos et al. unpublished).
For DNA barcoding analyses, temporary slides of individual nematodes were made in tap water and digital light microscope pictures were taken as a morphological voucher. The nematode was then transferred to a PCR tube with a solution containing 10 μl 0.05 M NaOH and 1 μl Tween20, heated for 15 min at 95 °C, and 40 μl of double-distilled water was added. PCR was carried out targeting either the 28S (large subunit) ribosomal RNA gene (nLSU) or the ribosomal internal transcribed spacer 2 (ITS2) locus, and PCR products were cleaned and sequenced directly. Forward and reverse primers for the nLSU were D2Ab (ACAAGTACCGTGAGGGAAAGTTG) and D3b (TCGGAAGGAACCAGCTACTA). For ITS2 we used VRAIN2F (CTTTGTACACACCGCCCGTCGCT) and VRAIN2R (TTTCACTCGCCGTTACTAAGG GAATC). The sequences obtained were 100% identical to published sequences for Caenorhabditis sp. 1 .
Genomic DNA was extracted from an inbred strain, JU1667, of C. monodelphis sp. n. (derived from strain SB341), maintained on E. coli OP50, using the proteinase K-spin column protocol (detailed in Additional file 1). Total RNA from the same culture was also extracted (methods detailed in Additional file 1). Two paired-end genomic libraries (insert sizes of 300 bp and 600 bp, respectively) and a single paired-end RNA-seq library (insert size 180 bp) were constructed using TruSeq reagents and sequenced on the Illumina HiSeq 2000 by Edinburgh Genomics. We obtained 124.3 million genomic read pairs (100 base, paired end) and 46.2 million pairs of RNA-Seq reads (also 100 base, paired end).
De novo genome assembly and gene prediction
Details of software versions and parameters are available (see Additional file 2). We performed initial quality control of our genomic sequence data using FastQC  and used Skewer  to remove low quality (Phred score < 30) and adapter sequence. Using blobtools , we generated taxon-annotated GC-coverage (TAGC) plots to identify and remove bacterial contamination. Sequence data were assembled with CLC assembler (CLCBio, Copenhagen, Denmark) and reads mapped back to this assembly using CLC mapper. Each assembly contig was compared to the NCBI Nucleotide (nt) database using megablast from the NCBI BLAST+ suite . Genomic read pairs were aligned to genome references from five E. coli (strains: BL21 (DE3), ETEC H10407, K12 substr. DH10B, K-12 substr. MC4100 and B str. REL606) using Bowtie , and aligned pairs discarded. We identified laboratory-induced contamination with Caenorhabditis elegans in the 600 bp insert library data. To remove this, we aligned read pairs of the uncontaminated 300 bp-insert library to the C. elegans N2 reference genome. Regions of similarity between the genomes of C. monodelphis sp. n. and C. elegans (i.e. those regions of C. elegans with aligned C. monodelphis sp. n. reads) were masked with Ns using BEDtools . Read pairs of the 600 bp-insert library were subsequently aligned to this masked C. elegans reference and any aligned read pairs discarded.
Cleaned sequence data were assembled with ABySS  (k = 83) and contigs were scaffolded with transcript evidence using SCUBAT . RepeatModeler  was used to identify repetitive regions which were then masked using RepeatMasker . RNA-Seq read pairs were aligned to the assembly using STAR , and the resulting BAM file was used to guide the prediction of protein-coding genes by BRAKER .
Gene structure comparisons
Genome sequences and annotation GFFs were downloaded from WormBase  and imported into a custom Ensembl database (version 84) . Using the Ensembl Perl API, the canonical transcript from each protein-coding gene was identified and exon and intron statistics were calculated. To compare the gene structures of C. monodelphis sp. n. with that of C. elegans, we identified all orthologous clusters (details below) in which C. monodelphis sp. n. and C. elegans proteins were present as single-copy. Exon and intron statistics were calculated for each gene pair, as described previously. Plots were generated using the ggplot2 package  and GenePainter .
Pairwise comparisons of protein sequences derived from genomic data for 23 species of Caenorhabditis and two outgroup species, Oscheius tipulae and Heterorhabditis bacteriophora, (see Additional file 3 for details) were performed using NCBI BLAST+  and clustered into orthologous groups using OrthoFinder . The sequences of 303 one-to-one orthologues (allowing for up to two species to have missing data) were extracted and aligned using ClustalOmega . Poorly aligned regions were removed from the alignments using trimAL  and trimmed alignments concatenated using FASconCAT  to yield a supermatrix. We performed maximum-likelihood (ML) analysis using RAxML  (PROTGTR + Γ substitution model) with 1,000 bootstrap replicates. Bayesian analysis was performed using PhyloBayes  (CAT-GTR), with two independent Markov chains, and convergence was assessed using Tracer .
This published work and the nomenclatural acts it contains have been registered in Zoobank: http://zoobank.org/urn:lsid:zoobank.org:pub:0E6F137B-9975-4A8E-91F2-D588A572076E. The LSID for this publication is: urn:lsid:zoobank.org:pub:0E6F137B-9975-4A8E-91 F2-D588A572076E.
Here we provide a formal description of SB341 as the type strain of C. monodelphis sp. n.
Measurements (in μm) of heat relaxed specimens of Caenorhabditis monodelphis sp. n.
870 ± 105
694 ± 36
456 ± 24
17.1 ± 0.8
22 ± 1.6
23 ± 1.2
4.9 ± 0.5
4.1 ± 0.3
3.6 ± 0.1
20.5 ± 2.6
22 ± 2.3
9.8 ± 0.7
1.99 ± 0.17
1.8 ± 0.2
3.9 ± 0.30
65 ± 1.8
51 ± 6.9
32 ± 3
20 ± 0.6
27 ± 2.3
27 ± 2
21 ± 1.1
1.9 ± 0.6
1.2 ± 0.2
0.6 ± 0.1
2.5 ± 0.2
2.4 ± 0.2
10 ± 0.7
9.7 ± 0.9
15 ± 1.5
15 ± 1.5
12.42 ± 1.6
13 ± 1.2
150 ± 6.7
141 ± 9.2
107 ± 3.3
55 ± 3.1
52 ± 3.6
26 ± 2.1
22.8 ± 1.1
39 ± 3.2
40 ± 4.9
Nerve ring to terminal bulb
11 ± 4.9
19 ± 3.3
Terminal bulb length
30 ± 1.8
27 ± 1.7
Diameter of median bulb
22 ± 2.5
17 ± 1.3
9 ± 0.5
Diameter of terminal bulb
25 ± 2
19 ± 1
11 ± 0.4
Anterior end to deirid
150 ± 8
150 ± 8.3
Postdeirid to anus
170 ± 29.8
141 ± 14
651 ± 100
494 ± 32
25 ± 2.6
24 ± 1.9
Anal body width
22 ± 2.1
17 ± 1.1
12 ± 0.6
43 ± 4.3
32 ± 3.2
46 ± 2.4
Anus to phasmid distance
26 ± 2.2
303 ± 68
342 ± 44
Gonad flexure length
226 ± 67
46 ± 6.8
45 ± 6.8
9.8 ± 1.3
53 ± 3.1
29 ± 2.9
25 ± 1
15 ± 0.9
Small species (female 0.72 - 1.04 mm, male 0.65 – 0.77 mm); cuticle thin, ca. 1 μm wide and finely annulated, 0.8 μm wide at midbody. Lateral field inconspicuous, about 9% of body width, consisting one ridge that can be traced anteriorly to the level of the median bulb and posteriorly at level of rectum in females and about 1½ spicules length anterior of the cloacal aperture in males. Six lips slightly protruding, each with one apical papilliform labial sensillum and a second circle of four sublateral cephalic sensilla in both sexes; amphids opening on the lateral lips, hardly discernible. Buccal tube long and slender, more than twice the width in lip region, pharyngeal sleeve envelopes nearly half of the stoma, the anterior as well as the posterior end of the tube appear slightly thickened, cheilostom inconspicuous, arcade cells forming the gymnostom sometimes visible; glottoid apparatus completely absent. Pharynx with a prominent median bulb, diameter more than 90% of diameter of terminal bulb; terminal bulb pyriform, with double chambered haustrulum, the anterior chamber smallish; cardia conspicuous, opens funnel-like in intestine. Nerve ring encircles isthmus in its anterior part in living specimens, more to the middle of the isthmus in heat relaxed or preserved specimens; deirids usually conspicuous in the lateral field at level of beginning of terminal bulb, sometimes not visible in heat relaxed animals; pore of excretory-secretory system hard to discern posterior of deirid level. Two gland cells ventral and slightly posterior of terminal bulb conspicuous in live specimens. Lateral canals visible in live specimens extending anteriorly to two stoma length from the anterior end and ending at rectum level in the female. Postdeirids usually very conspicuous dorsally of the lateral field at about 75% of body length in both sexes and about half the length between vulva and beginning of rectum (or at level of posterior end of uterus remnant) in females, sometimes not visible in heat relaxed specimens.
Maximum body diameter clearly anterior of the vulva, vulva position 65% body length, a transverse slit, bordered in both ends by cuticular longitudinal flaps, vulva lips moderately protruding, four diagonal vulval muscles conspicuous; one pseudocoelomocyte exists anterior of gonad flexure ventrally. Genital tracts asymmetrical; posterior branch rudimentary, sac like, on the left hand side of intestine, without flexure, almost as long as body diameter at the level of the vulva, containing spermatozoa (Fig. 2); anterior branch right of intestine, reflexed dorsally close to the pharynx, flexure more than half the length of the gonad (measured from vulva to flexure); at the flexure oocytes in several rows, downstream in one row, oocytes predominantly growing in the last position, where granules are stored inside; sphincter between oviduct and uterus, only a few sperm cells in oviduct, most of them in uterus and blind sac; oviparous, one egg at a time in uterus (rarely two), segmentation starts in the uterus. Rectum a little S-shaped, rectal gland cells very small, posterior anal lip slightly protuberant. Tail short, panagrolaimid, dorsally convex, with offset tip tapering, smooth to somewhat telescope-like by cuticle forming a sleeve-like structure; tail tip with tiny hooks, mostly one dorsal, but also subventral (compare with Poikilolaimus); opening of phasmids located at 60–65% of tail length, shortly anterior of tip, phasmid glands not reaching anus level.
Testis right of intestine, ventrally reflexed in a certain distance posterior of pharynx; flexure relatively short. One pseudocoelomocyte between pharynx and flexure ventrally. Bursa well developed, peloderan, anteriorly open, with smooth margin and sometimes terminally indented, posterior part of velum transversely striated. Nine pairs of genital papillae (GP) present, two of them anterior of the cloaca, genital papilla 1 (GP1) and GP2 spaced, GP3 to GP6 and GP7 to GP9 clustered, GP5 and GP7 point to the dorsal side of the velum, GP6 slightly bottle shaped, GP8 and GP9 fused at base, GP2 and GP8 not reaching the margin of velum. Phasmids forming small tubercles to the ventral side posterior of the last GP; formula of GPs: v1,v2/(v3,v4,ad,v5) (pd,v6,v7)ph. Precloacal sensillum small, precloacal lip simple (according to type A of W Sudhaus and K Kiontke ), postcloacal sensilla long filamentous. Spicules short and stout, tawny, separate, slightly curved, with prominent head; shaft with a transverse seam, with a prominent longitudinal ridge, a dorsal lamella, and an oval “window”, the tip notched. Gubernaculum dorsally projecting, flexible, in the distal part following the contour of the spicules, spoon shaped in ventral view.
Unsheathed, mouth closed; stoma long, slender. Pharyngeal sleeve covering about half of the stoma; pharynx with well-developed median and terminal bulbs; corpus length ca. 52% of pharynx length. Nerve ring somewhat in the middle between the middle and terminal bulb. Genital primordium at about 60% of body length, elongated oval in shape. Tail conical. Amphids, lateral lines, position excretory pore, deirids and phasmids not observed.
In one female a second set of “sensilla” were observed a short distance posterior to postdeirids, possibly a duplication of the postdeirids.
Type carrier and locality
Holotype and paratypes of Caenorhabditis monodelphis sp. n. were isolated from the tunnels of Cis castaneus (Herbst, 1793) (Ciidae, Coleoptera) in the bracket fungus Ganoderma applanatum (Polyporales) on a stump of the common beech (Fagus sylvatica) a few centimetres above the ground in Berlin-Grunewald in April 2001. The same sample included individuals of Diploscapter sp., Plectus sp., Oscheius dolichura and one individual dorylaimid and mononchid.
Holotype male (collection number WT 3684) and five female and four male paratypes (WT 3685, WT 3686) are deposited in the National Plant Protection Organization Wageningen, The Netherlands. In addition, four female and four male paratypes, are deposited in the collection of Museum Voor Dierkunde at Ghent University, Ghent, Belgium, five female and three male paratypes in Museum für Naturkunde an der Humboldt-Universität zu Berlin, Berlin, Germany. Additional paratypes are available in the UGent Nematode Collection (slides UGnem158, 159 & 160) of the Nematology Research Unit, Department of Biology, Ghent University, Ghent, Belgium.
Diagnosis and relationship
Caenorhabditis monodelphis sp. n. can be recognised as a Caenorhabditis based on the thickened GP6 and the clearly visible postdeirids. Caenorhabditis monodelphis sp. n. is distinguished from all other described Caenorhabditis species by the presence of a monodelphic genital tract in the female with a blind sac posterior the vulva, a panagrolaimid female tail shape, adults with only one ridge on the lateral field, a very long and slender stoma without visible glottoid apparatus and male with short, stout spicule with bifurcate tip.
Ecology and biology
Caenorhabditis monodelphis sp. n. is a gonochoristic species with both males and females. Females are oviparous and carry only one egg (rarely two eggs). Development from egg to adult took about 5–6 days in juice prepared from brown algae at room temperature. Development from dauer larva to adults was completed in less than 3 days at 20 °C on NA seeded with OP50. The lifespan of adults is at minimum 14 days for males and 17 days for females. One pair of adults produced 167 offspring in 8 days and the daily production of fertile eggs was 6–31 (mean 18; n = 14). After the reproductive phase, females lived 9–14 days (n = 3) with males present.
Caenorhabditis monodelphis sp. n. has until now only been found in Ganoderma and Fomes in Germany and Belgium in relation with the ciid beetle Cis castaneus. The Ganoderma carrying C. monodelphis sp. n. from Oslo was not investigated for the presence of C. castaneus. In fungal fruiting bodies lacking the beetle C. monodelphis sp. n. was not found. Dauers of C. monodelphis sp. n. were found under the elytra of the beetle, but were not found internally when the beetle was further dissected. These findings indicate a phoretic association with the beetle. As only dauer larvae were isolated from beetles, while adults and larvae were present in the fruiting bodies, we infer that C. monodelphis sp. n. exit from dauer within the mushroom, develop to adulthood and start to reproduce. The food source of the species in natural conditions is not known, but they survive and reproduce easily on E. coli OP50 in culture.
Genome sequence of an inbred strain of Caenorhabditis monodelphis sp. n.
Genome assembly statistics for C. monodelphis sp. n. and other Caenorhabditis species
CEGMA complete/partial (%)
Gene structure comparison of orthologous gene pairs from C. monodelphis sp. N. and C. elegans
C. monodelphis sp. n.
Gene length (bp)
Coding exon length (bp)
Coding exon count (n)
CDS span (bp)a
Intron length (bp)
Total intron span per gene (bp)
C. monodelphis sp. n. is sister to other known Caenorhabditis
Stemspecies pattern reconstruction
Our phylogenetic analyses were based on species with whole genome data available, and thus did not include the full known diversity of the genus. The stemspecies pattern was reconstructed based on ingroup and outgroup comparison. Previous molecular phylogenetic analyses of Caenorhabditis species using a small number of marker genes  placed C. monodelphis sp. n. and C. sonorae  as sister species, again arising at the base of the genus.
The following morphological synapomorphies can be hypothesised to support a C. monodelphis sp. n. – C. sonorae clade: mouth opening triangular (Fig. 4b), spicule having a complicated tip (notched or dentated) and a longish thin walled “window” in the blade (Figs. 1i, 4l), postcloacal sensilla being filiform (Fig. 4k), and the female tail shortened to less than three times anal body width. Other similarities between both these species are plesiomorphic.
By ingroup comparison we reconstruct the following characters of the stemspecies of Anarhabditis without differentiating them into apo- or plesiomorphies (on apomorphies see the legend of Fig. 7):
– adults of small size (less than 1 mm); − lips not offset from anterior end; − four cephalic sensilla present in male and female; − stoma with pharyngeal sleeve (stegostom length nearly that of gymnostom); − median bulb of pharynx strongly developed, corpus intima with transverse ridging, terminal bulb with double haustrulum; − gonochoristic; − female tail elongate conoid; − gonads amphidelphic, the anterior branch right and the posterior left of intestine; − vulva at midbody, a transverse slit; − oviparous, usually only one egg at a time in the uteri; − male gonad on the right side, reflexed to the ventral; − bursa peloderan and anteriorly open, oval-shaped in ventral view, with smooth margin, terminally not notched; − 9 pairs of even genital papillae, two precloacal largely spaced, GP3–6 evenly spaced, the last three GPs forming a tight cluster; GP1, GP5 and GP7 terminate on the dorsal surface of the bursa velum; − phasmids open behind GP9, inconspicuous; − bursa formula thus v1,v2/v3,v4,ad,v5 (pd,v6,v7)ph; − male tail tip present; − 1 + 2 circumcloacal sensilla inconspicuous, precloacal lip simple; − spicules separate, stout, head not rounded, behind the shaft a slight ventral projection, dorsal part of blade weakly cuticularised (velum), its tip possibly not even (argued below); − gubernaculum simple spatulate; − dauerlarvae with double cuticle (ensheathed), not waving.
Taxonomy of Caenorhabditis monodelphis sp. n.
Caenorhabditis monodelphis sp. n. is a new species of Caenorhabditis supported by its phylogenetic position as inferred from 303 molecular markers, morphology, habitat and specific association with Cis castaneus (Coleoptera). Morphologically, it could be confused with “Protorhabditis” species because of the absence of a clear glottoid apparatus. A glottoid apparatus has been lost 5–6 times independently within “Rhabditidae”  and, as illustrated here, also in C. monodelphis sp. n. This species resembles species from “Protorhabditis” with a very long stoma without glottoid apparatus, but differs from the Oxyuroides-group within “Protorhabditis” in having an open bursa and GP1 not anterior of the bursa. It is differentiated from the Xylocola-group within “Protorhabditis” in having nine genital papillae.
Previously, Caenorhabditis has been characterised as having the following apomorphic characteristics: the presence of a dorsal velum on the spicule, a lateral field with three ridges, an unsheathed dauer juvenile and a slightly thickened GP6 . With the discovery and description of C. monodelphis sp. n. the number of lateral ridges is no longer an apomorphic character of Caenorhabditis, since C. monodelphis sp. n. only has one lateral ridge.
Association with fungivorous beetles
Species of Caenorhabditis are known to occur in soil, compost, cadavers of insects, some plant material and the intestine of birds , and can most easily be isolated from rotting fruits, flowers and stems . Caenorhabditis elegans has also been found infesting cultures of the mushroom Agaricus bisporus . Wild mushrooms are an under-explored habitat for this genus, but our limited geographical sampling indicates that they could be an important habitat. Caenorhabditis monodelphis sp. n. was present in galleries made by Cis castaneus inside Ganoderma applanatum in Belgium, Norway, Germany and in an old fruiting body of Fomes fomentarius in Belgium. Although the true distribution of C. monodelphis sp. n. is not yet known, it is expected that this species will be found throughout Europe where Ganoderma (or in lesser extent Fomes) co-occurs with the mycophagous beetle Cis castaneus.
That Caenorhabditis species have phoretic relationships with insects and other invertebrates is well known . For C. monodelphis sp. n., all records are from mushroom fruiting bodies that are also inhabited by different insect groups, and dauer larvae were found under the elytra of Cis castaneus. Based on this evidence, we propose that C. monodelphis sp. n. propagates in galleries generated by Cisidae and the dauer larvae are transported by these beetles to uninfested mushrooms. Records of C. monodelphis sp. n. in both Ganoderma and Fomes, respectively the preferred  and the known  host indicate a beetle-specific rather than a mushroom-specific relationship. The only other known Caenorhabditis species which appears to be phoretically associated with fungivorous organisms, most likely insects, is C. auriculariae Tsuda & Futai,  of the Elegans super-group. This species was found only once in the fruit bodies of Auricularia polytricha (Agaricomycetes) in Japan, but the vector needed to infest the mushroom is unknown . C. elegans was also found to infest cultures of the champignon mushroom Agaricus bisporus , but most likely originated from mushroom compost where it can be frequently found. Several samples of different mushrooms on wood in Europe, USA and Japan did not yield other Caenorhabditis spp. However, given that many more insect species are known to feed and reproduce on mushrooms  and Rhabditida are known to use insects as a phoretic transport carrier , it is possible that mushroom species are habitats for many other rhabditid species, including new species of Caenorhabditis.
Genome sequence and gene structures of C. monodelphis sp. n.
Using next generation sequencing technologies and advanced bioinformatics toolkits, we have generated a good first-draft genome sequence for an inbred line derived from the type strain of C. monodelphis sp. n.. Although assembly metrics and CEGMA scores indicate the assembly is relatively contiguous and complete, it is likely that a proportion of C. monodelphis sp. n. genes are assembled only partially. This may have affected gene prediction, with the number of predicted gene models (17,180) being lower than estimates from most other Caenorhabditis species with available sequence data . Comparisons of orthologous gene pairs revealed a significant divergence in gene structure between C. monodelphis sp. n. and C. elegans. C. monodelphis sp. n. genes are typically longer, contain more coding exons and a longer span of introns than C. elegans genes (Table 3). This increase in gene length may, in part, account for the difference in genome span between C. monodelphis sp. n. and C. elegans. The clear trend towards more coding-exons in C. monodelphis sp. n. relative to C. elegans (Fig. 5) could be explained by extensive intron loss or gain in either species. Previous studies using a small number of genes have shown that intron losses have been far more common in Caenorhabditis evolution than intron gains [7, 52, 53]. Thus, it is possible that the gene structures seen in C. monodelphis sp. n. reflect an intron-rich ancestral state, and intron loss has predominated during the evolution of C. elegans. In Pristionchus pacificus, which is distantly related to Caenorhabditis, genes typically have roughly twice as many introns as their orthologues in C. elegans . Further analysis using genomes from more closely related outgroup species and other Caenorhabditis species is necessary before we can infer the dynamics of intron evolution in the genus.
Phylogenetic analysis of 303 clusters of putatively orthologous protein sequences derived from whole genome sequence data of 23 species of Caenorhabditis and two outgroup species resulted in a well resolved phylogenetic diagram and confirmation of C. monodelphis sp. n. as basal to all other analysed Caenorhabditis species (Fig. 6). The topology is largely congruent with previously published analyses performed using a smaller number of molecular loci . However, in contrast to the analyses of Kiontke et. al.  and Felix et. al.  which show C. brenneri and C. doughertyi as sister species, our phylogenetic hypothesis places C. doughertyi as more closely related to C. wallacei and C. tropicalis. This node, however, has low bootstrap support. Genome sequencing projects for several Caenorhabditis species, including those from the currently under-sampled Drosophilae super-group, are currently underway. These data will be essential to resolving the phylogenetic relationships of this important genus where morphology can be misinformative and/or misleading.
Reconstruction of stemspecies pattern
Details of our inference of ASR depends on the placement of P. elaphri in “Protorhabditis” versus as sister taxon of Anarhabditis (because of its distinct pharynx morphology) (Fig. 7). Molecular data resolving this issue are urgently required. Caenorhabditis monodelphis sp. n. and P. elaphri share a conspicuously long and narrow stoma due to an extended stegostom (long pharyngeal sleeve) without a glottoid apparatus (bulging of the three sectors of metastegostom) and one ridge in the lateral field. Based on current evidence, we interpret these peculiar similarities as homologous and thus as further characters of the Anarhabditis stemspecies as well as the Caenorhabditis stemspecies. The narrowing of the buccal cavity could have restricted the formation of sectoral swellings of the metastegostom, so that the typical glottoid apparatus disappeared. This happened in parallel in the rhabditid Matthesonema eremitum . The hypothesis of a reduction of the glottoid apparatus and its denticles in the stemspecies of Anarhabditis is in conflict with the structure of the metastegostom in most species of Caenorhabditis, where it looks like a transformation of a glottoid apparatus , and in C. sonorae is credibly described as a glottoid apparatus . To solve this conflict we must assume a partial reversion both in C. sonorae and in the sister-lineage of C. sonorae/C. monodelphis sp. n. However, instead of proposing two independent reversions, the possibility of an independent reduction of the glottoid apparatus in Protoscapter and C. monodelphis sp. n. remains an equally parsimonious alternative. A reinvestigation of P. elaphri could resolve this question.
In the stemspecies pattern of Anarhabditis the morphology of the tip of the spicules remains unclear. In the description of P. elaphri some drawings show the tip to be nearly pointed , but in other drawings (Figure twelve m of ) it is terminally notched. In P. tristis, I Andrássy  depicted a small terminal hook. These characters were not mentioned in the text in either species’ description. Nevertheless, the dentation of the spicule tips in the first branching Caenorhabditis sonorae/C. monodelphis sp. n. is different and distinct enough to judge this character as synapomorphic for these sister species (Fig. 7). Starting from the characters of the last common stemspecies of both these species, in C. sonorae the lateral ridge must have been reduced, so that its lateral field is smooth, and the male tail tip was retracted, so that the tail ends obtusely between the last GPs. In C. monodelphis sp. n. both these characters remain plesiomorphic, but the female posterior gonad branch is in the process of reduction. The ecological requirements of C. sonorae (inhabitant of cactus rot) and C. monodelphis sp. n. (living in the tunnels of Ciidae beetles in bracket fungi) are so different, that no statement on the ecology of their last common stemspecies is possible. However, as P. elaphri and C. monodelphis sp. n. exhibit a phoretic relationship with beetles and their dauer larvae seek a place under the elytra, we cautiously suggest that this behaviour could be found in the stemspecies of Anarhabditis and that of Caenorhabditis, respectively.
Transformations from the stemspecies pattern of Anarhabditis to Caenorhabditis can be traced in the cladogram (Fig. 7). With respect to the hypothesis of the stemspecies pattern of Caenorhabditis formulated by W Sudhaus and K Kiontke  only the character of the lateral field must be revised: a single ridge in the lateral field of adults must be assumed in the stemspecies pattern of Anarhabditis and of Caenorhabditis, respectively. Therefore, the evolution of three lateral cuticular ridges must have occurred first within Caenorhabditis (Fig. 7).
Degenerative evolution towards monodelphy
Uniquely for Caenorhabditis species, in C. monodelphis sp. n. the posterior female gonad branch has been reduced to a blind sac without gamete forming function. This vestigial branch serves mainly in storing sperm. In contrast to most mono-prodelphic rhabditids, the vulva is not shifted posteriorly in C. monodelphis sp. n.. A relict posterior gonad together with a nearly median vulva also occurs in Oscheius guentheri (Sudhaus & Hooper, 1994)  and an undescribed Diplogastrellus species from India (Sudhaus, unpublished data). Remarkable, in all these cases the anterior branch does not extend into the body posterior to the vulva, in contrast to monodelphic cephalobids, panagrolaimids and the rhabditid Rhabpanus. In Rhabpanus ossiculus Massey,  and R. uniquus Tahseen, Sultana, Khan & Hussain,  the prodelphic reflexed gonad reaches almost to the rectum while the vulva is located at 65–69% of body length and a short post-uterine sac filled with sperm is present [58, 59]. In contrast to species of Acrobeloides, Cephalobus, Mesorhabditis and Panagrolaimus, the posterior branch of the gonad of O. guentheri is not reduced by apoptosis of the distal tip cell , and the vestigial branch is very variable within this species . These patterns argue for a relatively recent reduction in O. guentheri. Based on the similarities (in the female gonad and nearly median vulva) between C. monodelphis sp. n. and O. guentheri, the gonadal system of female C. monodelphis sp. n. may also represent a relatively recent evolutionary shift.
The basal position and the unique characters of C. monodelphis sp. n. in the genus Caenorhabditis and its similarity with the hypothetical ancestor of Caenorhabditis makes C. monodelphis sp. n. a key species for future evolutionary and developmental studies within Caenorhabditis. Importantly we present here, alongside traditional morpological diagnosis of this new species a complete genome draft, which we believe is the first time this has been done for a metazoan species description. Release of the draft genome sequence of C. monodelphis sp. n., along with its formal description will, we hope, promote forward- and reverse-genetic analyses of its biology. In particular, CRISPR-Cas9 gene editing technologies, which require sequence knowledge for design of targeting oligonucleotides, are immediately facilitated.
While publication of marker sequence alongside species description is becoming commonplace , formal publication of whole genome data alongside species descriptions has historically been limited to prokaryotic taxa. In Eukaryota, this practice is just gaining traction, with the recent publication of the description of a fungal taxon with genome data (Epichloë inebrians, an ergot fungus ). Additionally, novel arthropod taxa used in phylogenomic analyses have had species descriptions published independently, but near-concurrently, with their genome data (Mengenilla moldrzyki, a strepsipteran insect [63, 64] or transcriptome data (the centipede Eupolybothrus cavernicolus ). For Caenorhabditis species, where morphology can be misinformative and/or misleading, phylogenomic analyses – and thus determination of genome sequence – will be essential for resolution of relationships. We suggest that genome scale data allied to species description should become commonplace.
Named after the monodelphic reproductive system in the female.
We are grateful to Dr. Karin Kiontke for providing a culture of Caenorhabditis monodelphis sp. n. We thank Dr. Marie-Anne Félix for sharing new data about Caenorhbaditis monodelphis sp. n. found in Norway and for the isolation of DNA and RNA samples from her inbred strain JU1667. We are thankful for Marjolein Couvreur for providing SEM images and Myriam Claeys for providing the TEM image. Sequencing was performed by Edinburgh Genomics, The University of Edinburgh. We also thank Roderic Page, Chrstopher Schardl, Christoph Bleidorn and Jason Stajich for responses on twitter concerning genome data allied to Eukaryotic species descriptions.
This work was supported by by a special research fund UGent 01 N02312 and the Foundation for Scientific Research, Flanders grant FWOKAN2013001201. Edinburgh Genomics is partly supported through core grants from NERC (R8/H10/56), MRC (MR/K001744/1) and BBSRC (BB/J004243/1). L.S. is funded by a Baillie Gifford Studentship, University of Edinburgh.
Availability of data and materials
Genome and transcriptome sequence read data of C. monodelphis sp. n. JU1677 are available from the European Nucleotide Archive and NCBI Short Read Archive under the accession PRJEB7905. The genome assembly and annotations are available to browse and download at ensembl.caenorhabditis.org and download.caenorhabditis.org, respectively. Accessions and links to genome-derived protein sequence data used in phylogenetic analyses are available in Additional file 3. Data files associated with this study have been deposited in Zenodo under the accession 10.5281/zenodo.160693.
DS, WS, LS and MB conceived the study. DS, WS, LS collected and analysed data and wrote the manuscript. All authors provided comments on early drafts of the manuscript. WB, LS and MB funded this study. All authors read, revised, and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Girard LR, Fiedler TJ, Harris TW, Carvalho F, Antoshechkin I, Han M, et al. WormBook: the online review of Caenorhabditis elegans biology. Nucleic Acids Res. 2007;35:D472–5.View ArticlePubMedGoogle Scholar
- Félix M-A, Braendle C, Cutter AD. A streamlined system for species diagnosis in Caenorhabditis (Nematoda: Rhabditidae) with name designations for 15 distinct biological species. PLoS ONE. 2014;9:e94723.
- Kiontke K, Félix M-A, Ailion M, Rockman M, Braendle C, Penigault J-B, et al. A phylogeny and molecular barcodes for Caenorhabditis, with numerous new species from rotting fruits. BMC Evol Biol. 2011;11:339.View ArticlePubMedPubMed CentralGoogle Scholar
- Abebe E, Mekete T, Thomas WK. A critique of current methods in nematode taxonomy. Afr J Biotechnol. 2013;10:312–23.Google Scholar
- Huys R, Llewellyn-Hughes J, Olson PD, Nagasawa K. Small subunit rDNA and Bayesian inference reveal Pectenophilus ornatus (Copepoda incertae sedis) as highly transformed Mytilicolidae, and support assignment of Chondracanthidae and Xarifiidae to Lichomolgoidea (Cyclopoida). Biol J Linn Soc. 2006;87:403–25.View ArticleGoogle Scholar
- Kiontke K, Barrière A, Kolotuev I, Podbilewicz B, Sommer R, Fitch DHA, et al. Trends, stasis, and drift in the evolution of nematode vulva development. Curr Biol. 2007;17:1925–37.View ArticlePubMedGoogle Scholar
- Kiontke K, Gavin NP, Raynes Y, Roehrig C, Piano F, Fitch DHA. Caenorhabditis phylogeny predicts convergence of hermaphroditism and extensive intron loss. Proc Natl Acad Sci U S A. 2004;101:9003–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Nuez I, Félix M-A. Evolution of susceptibility to ingested double-stranded RNAs in Caenorhabditis nematodes. PLoS ONE. 2012;7:e29811.View ArticlePubMedPubMed CentralGoogle Scholar
- Guevara R, Rayner ADM, Reynolds SE. Orientation of specialist and generalist fungivorous ciid beetles to host and non-host odours. Physiol Entomol. 2000;25:288–95.View ArticleGoogle Scholar
- Kiontke K, Sudhaus W. Ecology of Caenorhabditis species. In: Community TCeR: WormBook, editor. WormBook. 2006.Google Scholar
- Hooper DJ. Extraction of free-living stages from soil. In: Laboratory Methods for Work with Plant and Soil Nematodes. Edited by Southey JF. London: Her Majesty’s Stationery Office; 1986: 5–30.
- Yushin VV, Claeys M, Bert W. Ultrastructural immunogold localization of major sperm protein (MSP) in spermatogenic cells of the nematode Acrobeles complexus (Nematoda, Rhabditida). Micron. 2016;89:43–55.
- Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.Google Scholar
- Jiang HS, Lei R, Ding SW, Zhu SF. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinformatics. 2014;15:182.View ArticlePubMedPubMed CentralGoogle Scholar
- Laetsch D. Blobtools. https://github.com/DRL/blobtools.
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST plus: architecture and applications. BMC Bioinformatics. 2009;10:421.
- Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Meth. 2012;9:357–U354.View ArticleGoogle Scholar
- Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Simpson JT, Wong K, Jackman SD, Schein JE, Jones SJM, Birol I. ABySS: a parallel assembler for short read sequence data. Genome Res. 2009;19:1117–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Koutsovoulos G. SCUBAT2. https://github.com/GDKO/SCUBAT2.
- Smit AF, Hubley R. RepeatModeler Open-1.0. 2008–2015. http://www.repeatmasker.org.
- Smit AF, Hubley R, Green P. RepeatMasker Open-4.0. 1996–2010. http://www.repeatmasker.org.
- Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.View ArticlePubMedGoogle Scholar
- Hoff KJ, Lange S, Lomsadze A, Borodovsky M, Stanke M. BRAKER1: unsupervised RNA-Seq-based genome annotation with GeneMark-ET and AUGUSTUS. Bioinformatics. 2016;32:767–9.View ArticlePubMedGoogle Scholar
- Howe KL, Bolt BJ, Cain S, Chan J, Chen WJ, Davis P, et al. WormBase 2016: expanding to enable helminth genomic research. Nucleic Acids Res. 2016;44:D774–80.View ArticlePubMedGoogle Scholar
- Aken BL, Ayling S, Barrell D, Clarke L, Curwen V, Fairley S, et al. The Ensembl gene annotation system. Database. 2016;2016:baw093.
- Wickham H. ggplot2: Elegant Graphics for Data Analysis: Springer-Verlag New York: Springer; 2009.
- Mühlhausen S, Hellkamp M, Kollmar M. GenePainter v. 2.0 resolves the taxonomic distribution of intron positions. Bioinformatics. 2015;31:1302–4.View ArticlePubMedGoogle Scholar
- Emms DM, Kelly S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 2015;16:157.View ArticlePubMedPubMed CentralGoogle Scholar
- Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li WZ, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using clustal omega. Mol Syst Biol. 2011;7:539.View ArticlePubMedPubMed CentralGoogle Scholar
- Capella-Gutierrez S, Silla-Martinez JM, Gabaldon T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics. 2009;25:1972–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Kueck P, Longo GC. FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studies. Front Zool. 2014;11:81.View ArticleGoogle Scholar
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Lartillot N, Lepage T, Blanquart S. PhyloBayes 3: a Bayesian software package for phylogenetic reconstruction and molecular dating. Bioinformatics. 2009;25:2286–8.View ArticlePubMedGoogle Scholar
- Rambaut A, Suchard M, Xie D, Drummond A. Tracer v1. 6. 2014. http://tree.bio.ed.ac.uk/software/tracer/.
- Kiontke K, Fitch DHA. The phylogenetic relationships of Caenorhabditis and other rhabditids. In: Community TCeR: WormBook, editor. WormBook. 2005.Google Scholar
- Félix M-A. Cryptic quantitative evolution of the vulva intercellular signaling network in Caenorhabditis. Curr Biol. 2007;17:103–14.View ArticlePubMedGoogle Scholar
- Akimkina T, Yook K, Curnock S, Hodgkin J. Genome characterization, analysis of virulence and transformation of Microbacterium nematophilum, a coryneform pathogen of the nematode Caenorhabditis elegans. FEMS Microbiol Lett. 2006;264:145–51.View ArticlePubMedGoogle Scholar
- Sudhaus W, Kiontke K. Phylogeny of Rhabditis subgenus Caenorhabditis (Rhabditidae, Nematoda). J Zool Syst Evol Res. 1996;34:217–33.
- Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23:1061–7.View ArticlePubMedGoogle Scholar
- Kiontke K. Description of Rhabditis (Caenorhabditis) drosophilae n. sp. and R. (C.) sonorae n. sp. (Nematoda: Rhabditida) from saguaro cactus rot in Arizona. Fundam Appl Nematol. 1997;20:305–15.
- Sudhaus W. Phylogenetic systematisation and catalogue of paraphyletic “Rhabditidae” (Secernentea, Nematoda). J Nematode Morphol Syst. 2011;14:113–78.
- Sudhaus W. Vergleichende Untersuchungen zur Phylogenie, Systematik, Ökologie, Biologie und Ethologie der Rhabditidae (Nematoda). Zoologica. 1976;43:1–229.Google Scholar
- Sudhaus W, Fitch D. Comparative studies on the phylogeny and systematics of the Rhabditidae (Nematoda). J Nematol. 2001;33:1–70.
- Hirschmann H. Die Nematoden der Wassergrenze mittelfränkischer Gewässer. Zoologische Jahrbücher (Systematik). 1952;81:313–407.Google Scholar
- Sudhaus W. Order Rhabditina: “Rhabditidae”. In: Schmidt-Rhaesa A, editor. Handbook of zoology gastrotricha, cycloneuralia and gnathifera, vol. 2. Nematodath ed. Berlin, Boston: Walter De Gruyter; 2014. p. 537–55.
- Grewal PS, Richardson PN. Effects of Caenorhabditis elegans (Nematoda: Rhabditidae) on yield and quality of the cultivated mushroom Agaricus bisporus. Ann Appl Biol. 1991;118:381–94.
- Økland B. Insect fauna compared between six polypore species in a southern Norwegian spruce forest. Fauna Norv Ser B. 1995;42:21–6.Google Scholar
- Tsuda K, Futai K. Description of Caenorhabditis auriculariae n. sp. (Nematoda: Rhabditida) from fruiting bodies of Auricularia polytricha. Jpn J Nematol. 1999;29:18–23.
- Hammond PM, Lawrence JF. Appendix - Mycophagy in Insects: a Summary. In: Insect-fungus Interactions. Wilding N, Collins NM, Hammond PM, Webber JF, editors, vol. 14. London: Academic Press; 1989. p. 275–324.
- Timper P, Davies KG. Biotic interactions. In: Nematode Behaviour. Gaugler R, Bilgrami AL, editors. Wallingford, UK: CABI; 2004. p. 277–308.
- Cho S, Jin S-W, Cohen A, Ellis RE. A phylogeny of Caenorhabditis reveals frequent loss of introns during nematode evolution. Genome Res. 2004;14:1207–20.
- Hoogewijs D, De Henau S, Dewilde S, Moens L, Couvreur M, Borgonie G, et al. The Caenorhabditis globin gene family reveals extensive nematode-specific radiation and diversification. BMC Evol Biol. 2008;8:1.
- Dieterich C, Clifton SW, Schuster LN, Chinwalla A, Delehaunty K, Dinkelacker I, et al. The Pistionchus pacificus genome provides a unique perspective on nematode lifestyle and parasitism. Nat Genet. 2008;40:1193–8.
- Sudhaus W. Matthesonema eremitum n. sp. (Nematoda, Rhabditida) associated with hermit crabs (Coenobita) from the Philippines and its phylogenetic implications. Nematologica. 1986;32:247–55.
- Andrássy I. Erd- und Süßwassernematoden aus Bulgarien. Acta Zool Acad Sci Hung. 1958;4:1–88.Google Scholar
- Sudhaus W, Hooper DJ. Rhabditis (Oscheius) guentheri sp. n., an unusual species with reduced posterior ovary, with observations on the Dolichura and Insectivora groups (Nematoda: Rhabditidae). Nematologica. 1994;40:508–33.View ArticleGoogle Scholar
- Massey CL. Two new genera of nematodes parasitic in the eastern subterranean termite, Reticulitermes flavipes. J Invertebr Pathol. 1971;17:238–42.View ArticlePubMedGoogle Scholar
- Tahseen Q, Sultana R, Khan R, Hussain A. Description of two new and one known species of the closely related genera Artigas, 1927 and Massey, 1971 (Nematoda: Rhabditidae) with a discussion on their relationships. Nematology. 2012;14:555–70.
- Félix MA, Sternberg PW. Symmetry breakage in the development of one-armed gonads in nematodes. Development. 1996;122:2129–42.PubMedGoogle Scholar
- Sommer R, Carta LK, Kim S-y, Sternberg PW. Morphological, genetic and molecular description of Pristionchus pacificus sp. n. (Nematoda: Neodiplogasteridae). Fundam Appl Nematol. 1996;19:511–22.
- Chen L, Li XZ, Li CJ, Swoboda GA, Young CA, Sugawara K, et al. Two distinct Epichloë species symbiotic with Achnatherum inebrians, drunken horse grass. Mycologia. 2015;107:863–73.View ArticlePubMedGoogle Scholar
- Niehuis O, Hartig G, Grath S, Pohl H, Lehmann J, Tafer H, et al. Genomic and morphological evidence converge to resolve the enigma of Strepsiptera. Curr Biol. 2012;22:1309–13.
- Pohl H, Niehuis O, Gloyna K, Misof B, Beutel R. A new species of Mengenilla (Insecta, Strepsiptera) from Tunisia. ZooKeys. 2012;198:79–102.View ArticleGoogle Scholar
- Edmunds SC, Hunter CI, Smith V, Stoev P, Penev L. Biodiversity research in the “big data” era: GigaScience and pensoft work together to publish the most data-rich species description. GigaScience. 2013;2:14.View ArticlePubMedPubMed CentralGoogle Scholar