Nematodes have diverse reproductive strategies, which make them ideal subjects for comparative studies to address how mating systems evolve. Here we present the sex ratios and mating dynamics of the free-living nematode Rhabditis sp. SB347, in which males, females and hermaphrodites co-exist. The three sexes are produced by both selfing and outcrossing, and females tend to appear early in a mother’s progeny. Males prefer mating with females over hermaphrodites, which our results suggest is related to the female-specific production of the sex pheromones ascr#1 and ascr#9. We discuss the parallels between this system and that of parasitic nematodes that exhibit alternation between uniparental and biparental reproduction.

SummaryThree key steps in meiosis allow diploid organisms to produce haploid gametes: (1) homologous chromosomes (homologs) pair and undergo crossovers; (2) homologs segregate to opposite poles; and (3) sister chromatids segregate to opposite poles. The XX/XO sex determination system found in many nematodes [1] facilitates the study of meiosis because variation is easily recognized [2, 3, 4]. Here we show that meiotic segregation of X chromosomes in the trioecious nematode Auanema rhodensis [5] varies according to sex (hermaphrodite, female, or male) and type of gametogenesis (oogenesis or spermatogenesis). In this species, XO males exclusively produce X-bearing sperm [6, 7]. The unpaired X precociously separates into sister chromatids, which co-segregate with the autosome set to generate a functional haplo-X sperm. The other set of autosomes is discarded into a residual body. Here we explore the X chromosome behavior in female and hermaphrodite meioses. Whereas X chromosomes segregate following the canonical pattern during XX female oogenesis to yield haplo-X oocytes, during XX hermaphrodite oogenesis they segregate to the first polar body to yield nullo-X oocytes. Thus, crosses between XX hermaphrodites and males yield exclusively male progeny. During hermaphrodite spermatogenesis, the sister chromatids of the X chromosomes separate during meiosis I, and homologous X chromatids segregate to the functional sperm to create diplo-X sperm. Given these intra-species, intra-individual, and intra-gametogenesis variations in the meiotic program, A. rhodensis is an ideal model for studying the plasticity of meiosis and how it can be modulated.

Trioecy, a mating system in which males, females and hermaphrodites co-exist, is a useful system to investigate the origin and maintenance of alternative mating strategies. In the trioecious nematode Auanema rhodensis , males have one X chromosome (XO), whereas females and hermaphrodites have two (XX). The female vs. hermaphrodite sex determination mechanisms have remained elusive. In this study, RNA-seq analyses show a 20% difference between the L2 hermaphrodite and female gene expression profiles. RNAi experiments targeting the DM ( doublesex / mab-3 ) domain transcription factor dmd-10/11 suggest that the hermaphrodite sexual fate requires the upregulation of this gene. The genetic linkage map (GLM) shows that there is chromosome-wide heterozygosity for the X chromosome in F2 hermaphrodite-derived lines originated from crosses between two parental inbred strains. These results confirm the lack of recombination of the X chromosome in hermaphrodites, as previously reported. We also describe conserved chromosome elements (Nigon elements), which have been mostly maintained throughout the evolution of Rhabditina nematodes. The seven-chromosome karyotype of A. rhodensis , instead of the typical six found in other rhabditine species, derives from fusion/rearrangements events involving three Nigon elements. The A. rhodensis X chromosome is the smallest and most polymorphic with the least proportion of conserved genes. This may reflect its atypical mode of father-to-son transmission and its lack of recombination in hermaphrodites and males. In conclusion, this study provides a framework for studying the evolution of chromosomes in rhabditine nematodes, as well as possible mechanisms for the sex determination in a three-sexed species.

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

One limitation of the widely used RNA-seq method is that long transcripts are represented by more reads than shorter transcripts, resulting in a biased estimation of expression levels. The 3′ RNA-seq method, which yields only one sequence per transcript, bypasses this limitation. Here, RNA was extracted from two samples, in which we expected to find differentially expressed genes. Each was processed by both traditional and 3′ RNA-seq protocols. Both methods yielded similar differentially expressed genes and estimated expression levels in a comparable way, confirming they both represent valid tools for RNA-seq analysis. Notably, however, we identified more differentially expressed transcripts with the 3′ RNA-seq method, suggesting a greater power to detect expression variation using this method. Hence, when little genomic information is available for the species studied, the standard RNA-seq presents a better cost-benefit compromise, whereas for model species, the 3′ RNA-seq method might more accurately detect differential expression.

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, the number of chromosomes is related to the number of centromeres. However, in some groups, such as rhabditid nematodes, centromeric function is distributed across multiple sites on each chromosome. These holocentric chromosomes might, a priori, be expected to be permissive of large-scale chromosomal rearrangement, as chromosomal fragments could still partition correctly and fusions would not generate lethal conflict between multiple centromeres. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60 Mb O. tipulae genome is resolved into six chromosomal molecules. We find evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes we identify seven ancestral chromosomal elements (Nigon elements), and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex-chromosome associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

Background The emergence of insecticide resistance is a fast-paced example of the evolutionary process of natural selection. In this study, we investigated the molecular basis of resistance in the myiasis-causing fly Cochliomyia hominivorax (Diptera: Calliphoridae) to dimethyl-organophosphate (OP) insecticides. Methods By sequencing the RNA from surviving larvae treated with dimethyl-OP (resistant condition) and non-treated larvae (control condition), we identified genes displaying condition-specific polymorphisms, as well as those differentially expressed. Results Both analyses revealed that resistant individuals have altered expression and allele-specific expression of genes involved in proteolysis (specifically serine-endopeptidase), olfactory perception and cuticle metabolism, among others. We also confirmed that resistant individuals carry almost invariably the Trp251Ser mutation in the esterase E3, known to confer OP and Pyrethroid resistance. Interestingly, genes involved in metabolic and detoxifying processes (notably cytochrome P450s) were found under-expressed in resistant individuals. An exception to this were esterases, which were found up-regulated. Conclusions These observations suggest that reduced penetration and aversion to dimethyl-OP contaminated food may be important complementary strategies of resistant individuals. The specific genes and processes found are an important starting point for future functional studies. Their role in insecticide resistance merits consideration to better the current pest management strategies.

The New World Screwworm, Cochliomyia hominivorax (Calliphoridae), is the most important myiasis-causing species in America. Screwworm myiasis is a zoonosis that can cause severe lesions in livestock, domesticated and wild animals, and occasionally in people. Beyond the sanitary problems associated with this species, these infestations negatively impact economic sectors, such as the cattle industry. Here, we present a chromosome-scale assembly of C. hominivorax's genome, organized in 6 chromosome-length and 515 unplaced scaffolds and spanning 534 Mb. There was a clear correspondence between the D. melanogaster linkage groups A-E and the chromosomal-scale scaffolds. Chromosome Quotient (CQ) analysis identified a single scaffold from the X chromosome that contains most of the orthologs of genes that are on the D. melanogaster fourth chromosome (linkage group F or dot chromosome). CQ analysis also identified potential X and Y unplaced scaffolds and genes. Y-linkage for selected regions was confirmed by PCR with male and female DNA. Some of the long chromosome-scale scaffolds include Y-linked sequences, suggesting misassembly of these regions. These resources will provide a basis for future studies aiming at understanding the biology and evolution of this devastating obligate parasite.