Genus Caenorhabditis Research Articles

Life history in Caenorhabditis elegans: from molecular genetics to evolutionary ecology.

Life history is defined by traits that reflect key components of fitness, especially those relating to reproduction and survival. Research in life history seeks to unravel the relationships among these traits and understand how life history strategies evolve to maximize fitness. As such, life history research integrates the study of the genetic and developmental mechanisms underlying trait determination with the evolutionary and ecological context of Darwinian fitness. As a leading model organism for molecular and developmental genetics, Caenorhabditis elegans is unmatched in the characterization of life history-related processes, including developmental timing and plasticity, reproductive behaviors, sex determination, stress tolerance, and aging. Building on recent studies of natural populations and ecology, the combination of C. elegans' historical research strengths with new insights into trait variation now positions it as a uniquely valuable model for life history research. In this review, we summarize the contributions of C. elegans and related species to life history and its evolution. We begin by reviewing the key characteristics of C. elegans life history, with an emphasis on its distinctive reproductive strategies and notable life cycle plasticity. Next, we explore intraspecific variation in life history traits and its underlying genetic architecture. Finally, we provide an overview of how C. elegans has guided research on major life history transitions both within the genus Caenorhabditis and across the broader phylum Nematoda. While C. elegans is relatively new to life history research, significant progress has been made by leveraging its distinctive biological traits, establishing it as a highly cross-disciplinary system for life history studies.

Widespread changes in gene expression accompany body size evolution in nematodes.

Body size is a fundamental trait that drives multiple evolutionary and ecological patterns. Caenorhabditis inopinata is a fig-associated nematode that is exceptionally large relative to other members of the genus, including Caenorhabditis elegans. We previously showed that C. inopinata is large primarily due to postembryonic cell size expansion that occurs during the larval-to-adult transition. Here, we describe gene expression patterns in C. elegans and C. inopinata throughout this developmental period to understand the transcriptional basis of body size change. We performed RNA-seq in both species across the L3, L4, and adult stages. Most genes are differentially expressed across all developmental stages, consistent with C. inopinata's divergent ecology and morphology. We also used a model comparison approach to identify orthologues with divergent dynamics across this developmental period between the 2 species. This included genes connected to neurons, behavior, stress response, developmental timing, and small RNA/chromatin regulation. Multiple hypodermal collagens were also observed to harbor divergent developmental dynamics across this period, and genes important for molting and body morphology were also detected. Genes associated with transforming growth factor β signaling revealed idiosyncratic and unexpected transcriptional patterns given their role in body size regulation in C. elegans. This widespread transcriptional divergence between these species is unexpected and maybe a signature of the ecological and morphological divergence of C. inopinata. Alternatively, transcriptional turnover may be the rule in the Caenorhabditis genus, indicative of widespread developmental system drift among species. This work lays the foundation for future functional genetic studies interrogating the bases of body size evolution in this group.

Evolution of neuropeptide Y/RFamide-like receptors in nematodes

The Neuropeptide Y/RFamide-like receptors belong to the Rhodopsin-like G protein-coupled receptors G protein-coupled receptors (GPCRs) and are involved in functions such as locomotion, feeding and reproduction. With 41 described receptors they form the best-studied group of neuropeptide GPCRs in Caenorhabditis elegans. In order to understand the expansion of the Neuropeptide Y/RFamide-like receptor family in nematodes, we started from the sequences of selected receptor paralogs in C. elegans as query and surveyed the corresponding orthologous sequences in another 159 representative nematode target genomes. To this end we employed a automated pipeline based on ExonMatchSolver, a tool that solves the paralog-to-contig assignment problem. Utilizing subclass-specific HMMs we were able to detect a total of 1557 Neuropeptide Y/RFamide-like receptor sequences (1100 NPRs, 375 FRPRs and 82 C09F12.3) in the 159 target nematode genomes investigated here. These sequences demonstrate a good conservation of the Neuropeptide Y/RFamide-like receptors across the Nematoda and highlight the diversification of the family in nematode evolution. No other genus shares all Neuropeptide Y/RFamide-like receptors with the genus Caenorhabditis. At the same time, we observe large numbers of clade specific duplications and losses of family members across the phylum Nematoda.

Pervasive conservation of intron number and other genetic elements revealed by a chromosome-level genomic assembly of the hyper-polymorphic nematode Caenorhabditis brenneri.

With within-species genetic diversity estimates that span the gambit of that seen across the entirety of animals, the Caenorhabditis genus of nematodes holds unique potential to provide insights into how population size and reproductive strategies influence gene and genome organization and evolution. Our study focuses on Caenorhabditis brenneri, currently known as one of the most genetically diverse nematodes within its genus and metazoan phyla. Here, we present a high-quality gapless genome assembly and annotation for C. brenneri, revealing a common nematode chromosome arrangement characterized by gene-dense central regions and repeat rich peripheral parts. Comparison of C. brenneri with other nematodes from the 'Elegans' group revealed conserved macrosynteny but a lack of microsynteny, characterized by frequent rearrangements and low correlation iof orthogroup sizes, indicative of high rates of gene turnover. We also assessed genome organization within corresponding syntenic blocks in selfing and outcrossing species, affirming that selfing species predominantly experience loss of both genes and intergenic DNA. Comparison of gene structures revealed strikingly small number of shared introns across species, yet consistent distributions of intron number and length, regardless of population size or reproductive mode, suggesting that their evolutionary dynamics are primarily reflective of functional constraints. Our study provides valuable insights into genome evolution and expands the nematode genome resources with the highly genetically diverse C. brenneri, facilitating research into various aspects of nematode biology and evolutionary processes.

Gene duplication and evolutionary plasticity of lin-12/Notch gene function in Caenorhabditis.

Gene duplication is an important substrate for the evolution of new gene functions, but the impacts of gene duplicates on their own activities and on the developmental networks in which they act are poorly understood. Here, we use a natural experiment of lin-12/Notch gene duplication within the nematode genus Caenorhabditis, combined with characterization of loss- and gain-of-function mutations, to uncover functional distinctions between the duplicate genes in 1 species (Caenorhabditis briggsae) and their single-copy ortholog in Caenorhabditis elegans. First, using improved genomic sequence and gene model characterization, we confirm that the C. briggsae genome includes 2 complete lin-12 genes, whereas most other genes encoding proteins that participate in the LIN-12 signaling pathway retain a one-to-one orthology with C. elegans. We use CRISPR-mediated genome editing to introduce alleles predicted to cause gain-of-function (gf) or loss-of-function (lf) into each C. briggsae gene and find that the gf mutations uncover functional distinctions not apparent from the lf alleles. Specifically, Cbr-lin-12.1(gf), but not Cbr-lin-12.2(gf), causes developmental defects similar to those observed in Cel-lin-12(gf). In contrast to Cel-lin-12(gf), however, the Cbr-lin-12.1(gf) alleles do not cause dominant phenotypes as compared to the wild type, and the mutant phenotype is observed only when 2 gf alleles are present. Our results demonstrate that gene duplicates can exhibit differential capacities to compensate for each other and to interfere with normal development, and uncover coincident gene duplication and evolution of developmental sensitivity to LIN-12/Notch activity.

Wild-type Caenorhabditis sinica, a model nematode for speciation and evolution, massively found in Vietnam

Caenorhabditis sinica is a male/female species, in which the genetic diversity is possibly high. Thus, the species is advantageous for the understanding of the mechanisms of diversity, evolution, and adaptation within the nematode genus Caenorhabditis. Previous studies reported the geographic distribution of C. sinica only in China; nonetheless, this should be more convincing for its ecology. We allegedly surveyed the Caenorhabditis species and their diversity in the forests of Vietnam. We found 59 C. sinica strains, and this therefore indicates the vast variation that is 85.88% and 100% identical within isolated wild-type strains. They unequally present in habitats and are enormous in the northern forest, and rarer in the southern forest. The comparison of 18S rDNA barcode sequences from 59 C. sinica isolates by nucleotide sequence alignment showed the consistent diversity among the strains in or off the same ecologies, and all are comparable with the ever-first isolated strain JU727.

Comparative Genomics of Sex, Chromosomes, and Sex Chromosomes in Caenorhabditis elegans and Other Nematodes.

The nematode phylum has evolved a remarkable diversity of reproductive modes, including the repeated emergence of asexuality and hermaphroditism across divergent clades. The species-richness and small genome size of nematodes make them ideal systems for investigating the genome-wide causes and consequences of such major transitions. The availability of functional annotations for most Caenorhabditis elegans genes further allows the linking of patterns of gene content evolution with biological processes. Such gene-centric studies were recently complemented by investigations of chromosome evolution that made use of the first chromosome-scale genome assemblies outside the Caenorhabditis genus. This review highlights recent comparative genomic studies of reproductive mode evolution addressing the hybrid origin of asexuality and the parallel gene loss following the emergence of hermaphroditism. It further summarizes ongoing efforts to characterize ancient linkage blocks called Nigon elements, which form central units of chromosome evolution. Fusions between Nigon elements have been demonstrated to impact recombination and speciation. Finally, multiple recent fusions between autosomal and thesex-linked Nigon element reveal insights into the dynamic evolution of sex chromosomes across various timescales.

A germline-specific role for unconventional components of the γ-tubulin complex in Caenorhabditis elegans.

The γ-tubulin complex (γTuC) is a widely conserved microtubule nucleator, but some of its components, namely GCP4, GCP5 and GCP6 (also known as TUBGCP4, TUBGCP5 and TUBGCP6, respectively), have not been detected in Caenorhabditis elegans. Here, we identified two γTuC-associated proteins in C. elegans, GTAP-1 and GTAP-2, for which apparent orthologs were detected only in the genus Caenorhabditis. GTAP-1 and GTAP-2 were found to localize at centrosomes and the plasma membrane of the germline, and their centrosomal localization was interdependent. In early C. elegans embryos, whereas the conserved γTuC component MZT-1 (also known as MOZART1 and MZT1) was essential for the localization of centrosomal γ-tubulin, depletion of GTAP-1 and/or GTAP-2 caused up to 50% reduction of centrosomal γ-tubulin and precocious disassembly of spindle poles during mitotic telophase. In the adult germline, GTAP-1 and GTAP-2 contributed to efficient recruitment of the γTuC to the plasma membrane. Depletion of GTAP-1, but not of GTAP-2, severely disrupted both the microtubule array and the honeycomb-like structure of the adult germline. We propose that GTAP-1 and GTAP-2 are unconventional components of the γTuC that contribute to the organization of both centrosomal and non-centrosomal microtubules by targeting the γTuC to specific subcellular sites in a tissue-specific manner.

Daf-42 is an evolutionarily young gene essential for dauer development in Caenorhabditis elegans.

Under adverse environmental conditions, nematodes arrest into dauer, an alternative developmental stage for diapause. Dauer endures unfavorable environments and interacts with host animals to access favorable environments, thus playing a critical role in survival. Here, we report that in Caenorhabditis elegans, daf-42 is essential for development into the dauer stage, as the null mutant of daf-42 exhibited a "no viable dauer" phenotype in which no viable dauers were obtained in any dauer-inducing conditions. Long-term time lapse microscopy of synchronized larvae revealed that daf-42 is involved in developmental changes from the pre-dauer L2d stage to the dauer stage. daf-42 encodes large, disordered proteins of various sizes that are expressed in and secreted from the seam cells within a narrow time window shortly before the molt into dauer stage. Transcriptome analysis showed that the transcription of genes involved in larval physiology and dauer metabolism is highly affected by the daf-42 mutation. Contrary to the notion that essential genes that control the life and death of an organism may be well conserved across diverse species, daf-42 is an evolutionarily young gene conserved only in the Caenorhabditis genus. Our study shows that dauer formation is a vital process that is controlled not only by conserved genes but also by newly emerged genes, providing important insights into evolutionary mechanisms.

Widespread sex ratio polymorphism inCaenorhabditisnematodes

Although equal sex ratio is ubiquitous and represents an equilibrium in evolutionary theory, biased sex ratios are predicted for certain local conditions. Cases of sex ratio bias have been mostly reported for single species, but little is known about its evolution above the species level. Here, we surveyed progeny sex ratios in 23 species of the nematode genus Caenorhabditis, including 19 for which we tested multiple strains. For the species with multiple strains, five species had female-biased and two had non-biased sex ratios in all strains, respectively. The other 12 species showed polymorphic sex ratios across strains. Female-biased sex ratios could be due to sperm competition whereby X-bearing sperm outcompete nullo-X sperm during fertilization. In this model, when sperm are limited allowing all sperm to be used, sex ratios are expected to be equal. However, in assays limiting mating to a few hours, most strains showed similarly biased sex ratios compared with unlimited mating experiments, except that one C. becei strain showed significantly reduced female bias compared with unlimited mating. Our study shows frequent polymorphism in sex ratios within Caenorhabditis species and that sperm competition alone cannot explain the sex ratio bias.

Widespread transposon co-option in the Caenorhabditis germline regulatory network.

The movement of selfish DNA elements can lead to widespread genomic alterations with potential to create novel functions. We show that transposon expansions in Caenorhabditis nematodes led to extensive rewiring of germline transcriptional regulation. We find that about one-third of Caenorhabditis elegans germline-specific promoters have been co-opted from two related miniature inverted repeat transposable elements (TEs), CERP2 and CELE2. These promoters are regulated by HIM-17, a THAP domain-containing transcription factor related to a transposase. Expansion of CERP2 occurred before radiation of the Caenorhabditis genus, as did fixation of mutations in HIM-17 through positive selection, whereas CELE2 expanded only in C. elegans. Through comparative analyses in Caenorhabditis briggsae, we find not only evolutionary conservation of most CERP2 co-opted promoters but also a substantial fraction that are species-specific. Our work reveals the emergence and evolutionary conservation of a novel transcriptional network driven by TE co-option with a major impact on regulatory evolution.

The compact genome of Caenorhabditis niphades n. sp., isolated from a wood-boring weevil, Niphades variegatus

BackgroundThe first metazoan genome sequenced, that of Caenorhabditis elegans, has motivated animal genome evolution studies. To date > 50 species from the genus Caenorhabditis have been sequenced, allowing research on genome variation.ResultsIn the present study, we describe a new gonochoristic species, Caenorhabditis niphades n. sp., previously referred as C. sp. 36, isolated from adult weevils (Niphades variegatus), with whom they appear to be tightly associated during its life cycle. Along with a species description, we sequenced the genome of C. niphades n. sp. and produced a chromosome-level assembly. A genome comparison highlighted that C. niphades n. sp. has the smallest genome (59 Mbp) so far sequenced in the Elegans supergroup, despite being closely related to a species with an exceptionally large genome, C. japonica.ConclusionsThe compact genome of C. niphades n. sp. can serve as a key resource for comparative evolutionary studies of genome and gene number expansions in Caenorhabditis species.

Complete Mitogenome of Cruznema Tripartitum Confirms Highly Conserved Gene Arrangement within Family Rhabditidae.

Mitochondrial genomes have widely been used as molecular markers in understanding the patterns and processes of nematode evolution. The species in genus Cruznema are free-living bacterivores as well as parasites of crickets and mollusks. The complete mitochondrial genome of C. tripartitum was determined through high-throughput sequencing as the first sequenced representative of the genus Cruznema. The genome is comprised of 14,067 bp nucleotides, and includes 12 protein-coding, two rRNA, and 22 tRNA genes. Phylogenetic analyses based on amino acid data support C. tripartitum as a sister to the clade containing Caenorhabditis elegans and Oscheius chongmingensis. The analysis of gene arrangement suggested that C. tripartitum shares the same gene order with O. chongmingensis, Litoditis marina, Diplocapter coronatus, genus Caenorhabditis, and Pristionchus pacificus. Thus, the mitochondrial gene arrangement is highly conserved in the family Rhabditidae as well as some species in Diplogasteridae.

Chromosome-Level Reference Genomes for Two Strains of Caenorhabditis briggsae: An Improved Platform for Comparative Genomics.

The publication of the Caenorhabditis briggsae reference genome in 2003 enabled the first comparative genomics studies between C. elegans and C. briggsae, shedding light on the evolution of genome content and structure in the Caenorhabditis genus. However, despite being widely used, the currently available C. briggsae reference genome is substantially less complete and structurally accurate than the C. elegans reference genome. Here, we used high-coverage Oxford Nanopore long-read and chromosome-conformation capture data to generate chromosome-level reference genomes for two C. briggsae strains: QX1410, a new reference strain closely related to the laboratory AF16 strain, and VX34, a highly divergent strain isolated in China. We also sequenced 99 recombinant inbred lines generated from reciprocal crosses between QX1410 and VX34 to create a recombination map and identify chromosomal domains. Additionally, we used both short- and long-read RNA sequencing data to generate high-quality gene annotations. By comparing these new reference genomes to the current reference, we reveal that hyper-divergent haplotypes cover large portions of the C. briggsae genome, similar to recent reports in C. elegans and C. tropicalis. We also show that the genomes of selfing Caenorhabditis species have undergone more rearrangement than their outcrossing relatives, which has biased previous estimates of rearrangement rate in Caenorhabditis. These new genomes provide a substantially improved platform for comparative genomics in Caenorhabditis and narrow the gap between the quality of genomic resources available for C. elegans and C. briggsae.

Genome-wide characterization, evolution, structure, and expression analysis of the F-box genes in Caenorhabditis

BackgroundF-box proteins represent a diverse class of adaptor proteins of the ubiquitin-proteasome system (UPS) that play critical roles in the cell cycle, signal transduction, and immune response by removing or modifying cellular regulators. Among closely related organisms of the Caenorhabditis genus, remarkable divergence in F-box gene copy numbers was caused by sizeable species-specific expansion and contraction. Although F-box gene number expansion plays a vital role in shaping genomic diversity, little is known about molecular evolutionary mechanisms responsible for substantial differences in gene number of F-box genes and their functional diversification in Caenorhabditis. Here, we performed a comprehensive evolution and underlying mechanism analysis of F-box genes in five species of Caenorhabditis genus, including C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei.ResultsHerein, we identified and characterized 594, 192, 377, 39, 1426 F-box homologs encoding putative F-box proteins in the genome of C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei, respectively. Our work suggested that extensive species-specific tandem duplication followed by a small amount of gene loss was the primary mechanism responsible for F-box gene number divergence in Caenorhabditis genus. After F-box gene duplication events occurred, multiple mechanisms have contributed to gene structure divergence, including exon/intron gain/loss, exonization/pseudoexonization, exon/intron boundaries alteration, exon splits, and intron elongation by tandem repeats. Based on high-throughput RNA sequencing data analysis, we proposed that F-box gene functions have diversified by sub-functionalization through highly divergent stage-specific expression patterns in Caenorhabditis species.ConclusionsMassive species-specific tandem duplications and occasional gene loss drove the rapid evolution of the F-box gene family in Caenorhabditis, leading to complex gene structural variation and diversified functions affecting growth and development within and among Caenorhabditis species. In summary, our findings outline the evolution of F-box genes in the Caenorhabditis genome and lay the foundation for future functional studies.

Balancing selection maintains hyper-divergent haplotypes in Caenorhabditis elegans.

Across diverse taxa, selfing species have evolved independently from outcrossing species thousands of times. The transition from outcrossing to selfing significantly decreases the effective population size, effective recombination rate, and heterozygosity within a species. These changes lead to a reduction in genetic diversity, and therefore adaptive potential, by intensifying the effects of random genetic drift and linked selection. Within the nematode genus Caenorhabditis, selfing has evolved at least three times and all three species, including in the model organism Caenorhabditis elegans, show substantially reduced genetic diversity relative to outcrossing species. Selfing and outcrossing Caenorhabditis species are often found in the same niches, but we still do not know how selfing species with limited genetic diversity can adapt to these environments. Here, we examine the whole-genome sequences from 609 wild C. elegans strains isolated worldwide and show that genetic variation is concentrated in punctuated hyper-divergent regions that cover 20% of the C. elegans reference genome. These regions are enriched in environmental response genes that mediate sensory perception, pathogen response, and xenobiotic stress response. Population genomic evidence suggests that genetic diversity in these regions has been maintained by long-term balancing selection. Using long- read genome assemblies for 15 wild strains, we show that hyper-divergent haplotypes contain unique sets of genes and show levels of divergence comparable to levels found between Caenorhabditis species that diverged millions of years ago. These results provide an example for how species can avoid the evolutionary “dead end” associated with selfing.

Additional description and genome analyses of Caenorhabditis auriculariae representing the basal lineage of genus Caenorhabditis

Caenorhabditis auriculariae, which was morphologically described in 1999, was re-isolated from a Platydema mushroom-associated beetle. Based on the re-isolated materials, some morphological characteristics were re-examined and ascribed to the species. In addition, to clarify phylogenetic relationships with other Caenorhabditis species and biological features of the nematode, the whole genome was sequenced and assembled into 109.5 Mb with 16,279 predicted protein-coding genes. Molecular phylogenetic analyses based on ribosomal RNA and 269 single-copy genes revealed the species is closely related to C. sonorae and C. monodelphis placing them at the most basal clade of the genus. C. auriculariae has morphological characteristics clearly differed from those two species and harbours a number of species-specific gene families, indicating its usefulness as a new outgroup species for Caenorhabditis evolutionary studies. A comparison of carbohydrate-active enzyme (CAZy) repertoires in genomes, which we found useful to speculate about the lifestyle of Caenorhabditis nematodes, suggested that C. auriculariae likely has a life-cycle with tight-association with insects.

A cohort of Caenorhabditis species lacking the highly conserved let-7 microRNA.

The let-7 gene encodes a highly conserved microRNA with critical functions integral to cell fate specification and developmental progression in diverse animals. In Caenorhabditis elegans, let-7 is a component of the heterochronic (developmental timing) gene regulatory network, and loss-of-function mutations of let-7 result in lethality during the larval to adult transition due to misregulation of the conserved let-7 target, lin-41. To date, no bilaterian animal lacking let-7 has been characterized. In this study, we identify a cohort of nematode species within the genus Caenorhabditis, closely related to C. elegans, that lack the let-7 microRNA, owing to absence of the let-7 gene. Using Caenorhabditis sulstoni as a representative let-7-lacking species to characterize normal larval development in the absence of let-7, we demonstrate that, except for the lack of let-7, the heterochronic gene network is otherwise functionally conserved. We also report that species lacking let-7 contain a group of divergent let-7 paralogs—also known as the let-7-family of microRNAs—that have apparently assumed the role of targeting the LIN-41 mRNA.

Cultivation of Caenorhabditis elegans on new cheap monoxenic media without peptone.

The study of species biodiversity within the Caenorhabditis genus of nematodes would be facilitated by the isolation of as many species as possible. So far, over 50 species have been found, usually associated with decaying vegetation or soil samples, with many from Africa, South America and Southeast Asia. Scientists based in these regions can contribute to Caenorhabditis sampling and their proximity would allow intensive sampling, which would be useful for understanding the natural history of these species. However, severely limited research budgets are often a constraint for these local scientists. In this study, we aimed to find a more economical, alternative growth media to rear Caenorhabditis and related species. We tested 25 media permutations using cheaper substitutes for the reagents found in the standard nematode growth media (NGM) and found three media combinations that performed comparably to NGM with respect to the reproduction and longevity of C. elegans. These new media should facilitate the isolation and characterization of Caenorhabditis and other free-living nematodes for the researchers in the poorer regions such as Africa, South America, and Southeast Asia where nematode diversity appears high.

Gene expression analysis of oxidative stress tolerance in the entomopathogenic nematode Heterorhabditis bacteriophora

Summary The entomopathogenic nematode (EPN) Heterorhabditis bacteriophora is used as a biological control agent against diverse insect pests. Nematode-based products contain third-stage Dauer juveniles (DJ), which are adapted to long-term survival. Their longevity can be limited by stress conditions during industrial production up to field application. Oxidative stress has been recently reported to have a direct influence on the H. bacteriophora DJ longevity. Thus, understanding mechanisms by which DJ respond to oxidative stress can provide insights to improve DJ longevity. In this study, we carried out a comparative transcriptomic analysis on the early stage of oxidative stress induction (4 h) in two H. bacteriophora inbred lines with contrasting oxidative stress tolerance, HU2-IL1 (stress-tolerant) and PT1-IL1 (stress-sensitive). For assessing the transcriptome, MACE, a versatile RNA-seq variant was applied. Our de novo transcriptome assembly generated more than 20 000 transcripts, from which 10 290 were linked to 9776 different Uniprot accessions. The majority of the annotated transcripts presented high homology to parasitic nematodes within the genus Ancylostoma, whereas homology to the genus Caenorhabditis was negligible. A total of 630 and 461 transcripts were up-regulated (log2 fold-change (FC) ⩾ 2.0) in the stress-sensitive and the stress-tolerant line, respectively. The proportion of down-regulated transcripts was higher for both lines. However, down-regulation in the stress-sensitive line (5207 transcripts) exhibited a larger proportion than in the tolerant line (1844 transcripts), which indicates that targeted suppression of biological processes is also a crucial factor for the survival of H. bacteriophora under oxidative stress. Our global view of the transcriptome remodelling under oxidative stress suggests that the stress-sensitive line fails to maintain vital biological processes in contrast to the tolerant line. Interestingly, both lines activated similar biological processes directly involved in detoxification of reactive oxygen species, indicating that expression changes on transcript isoforms are of high relevance in this context. This study will open ways for the selection of DJ longevity predictor genes and allow the design of molecular markers for the breeding of improved lines.