A gene with a thousand alleles: The hyper-variable effectors of plant-parasitic nematodes.

Pathogens are engaged in a fierce evolutionary arms race with their host. The genes at the forefront of the engagement between kingdoms are often part of diverse and highly mutable gene families. Even in this context, we discovered unprecedented variation in the hyper-variable (HYP) effectors of plant-parasitic nematodes. HYP effectors are single-gene loci that potentially harbor thousands of alleles. Alleles vary in the organization, as well as the number, of motifs within a central hyper-variable domain (HVD). We dramatically expand the HYP repertoire of two plant-parasitic nematodes and define distinct species-specific "rules" underlying the apparently flawless genetic rearrangements. Finally, by analyzing the HYPs in 68 individual nematodes, we unexpectedly found that despite the huge number of alleles, most individuals are germline homozygous. These data support a mechanism of programmed genetic variation, termed HVD editing, where alterations are locus specific, strictly governed by rules, and theoretically produce thousands of variants without errors.

Oncogenic c-Myc induces replication stress by increasing cohesins chromatin occupancy in a CTCF-dependent manner.

Oncogene-induced replication stress is a crucial driver of genomic instability and one of the key events contributing to the onset and evolution of cancer. Despite its critical role in cancer, the mechanisms that generate oncogene-induced replication stress remain not fully understood. Here, we report that an oncogenic c-Myc-dependent increase in cohesins on DNA contributes to the induction of replication stress. Accumulation of cohesins on chromatin is not sufficient to cause replication stress, but also requires cohesins to accumulate at specific sites in a CTCF-dependent manner. We propose that the increased accumulation of cohesins at CTCF site interferes with the progression of replication forks, contributing to oncogene-induced replication stress. This is different from, and independent of, previously suggested mechanisms of oncogene-induced replication stress. This, together with the reported protective role of cohesins in preventing replication stress-induced DNA damage, supports a double-edge involvement of cohesins in causing and tolerating oncogene-induced replication stress.

Genomic Signatures Associated with Transitions to Viviparity in Cyprinodontiformes.

The transition from oviparity to viviparity has occurred independently over 150 times across vertebrates, presenting one of the most compelling cases of phenotypic convergence. However, whether the repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common signature in genomic divergence remains largely unknown. Although recent investigations into the evolution of viviparity have demonstrated striking similarity among the genes and molecular pathways involved across disparate vertebrate groups, quantitative tests for genome-wide convergent have provided ambivalent answers. Here, we investigate the potential role of molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo genomes and utilized publicly available genomes of viviparous and oviparous species to test for molecular convergence across both coding and noncoding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and in rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, both statistical power and biological confounds could constrain our ability to detect significant correlated evolution. We therefore identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. Motif enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development, and immunity occurred alongside the evolution of viviparity. Overall, however, our findings indicate that independent transitions to viviparity in these fish are not strongly associated with an excess of molecular convergence, but a few genes show convincing evidence of convergent evolution.

THAPBI PICT-a fast, cautious, and accurate metabarcoding analysis pipeline.

THAPBI PICT is an open source software pipeline for metabarcoding analysis of Illumina paired-end reads, including cases of multiplexing where more than one amplicon is amplified per DNA sample. Initially a Phytophthora ITS1 Classification Tool (PICT), we demonstrate using worked examples with our own and public data sets how, with appropriate primer settings and a custom database, it can be applied to other amplicons and organisms, and used for reanalysis of existing datasets. The core dataflow of the implementation is (i) data reduction to unique marker sequences, often called amplicon sequence variants (ASVs), (ii) dynamic thresholds for discarding low abundance sequences to remove noise and artifacts (rather than error correction by default), before (iii) classification using a curated reference database. The default classifier assigns a label to each query sequence based on a database match that is either perfect, or a single base pair edit away (substitution, deletion or insertion). Abundance thresholds for inclusion can be set by the user or automatically using per-batch negative or synthetic control samples. Output is designed for practical interpretation by non-specialists and includes a read report (ASVs with classification and counts per sample), sample report (samples with counts per species classification), and a topological graph of ASVs as nodes with short edit distances as edges. Source code available from https://github.com/peterjc/thapbi-pict/ with documentation including installation instructions.

Misregulation of cell cycle-dependent methylation of budding yeast CENP-A contributes to chromosomal instability.

Centromere (CEN) identity is specified epigenetically by specialized nucleosomes containing evolutionarily conserved CEN-specific histone H3 variant CENP-A (Cse4 in Saccharomyces cerevisiae, CENP-A in humans), which is essential for faithful chromosome segregation. However, the epigenetic mechanisms that regulate Cse4 function have not been fully defined. In this study, we show that cell cycle-dependent methylation of Cse4-R37 regulates kinetochore function and high-fidelity chromosome segregation. We generated a custom antibody that specifically recognizes methylated Cse4-R37 and showed that methylation of Cse4 is cell cycle regulated with maximum levels of methylated Cse4-R37 and its enrichment at the CEN chromatin occur in the mitotic cells. Methyl-mimic cse4-R37F mutant exhibits synthetic lethality with kinetochore mutants, reduced levels of CEN-associated kinetochore proteins and chromosome instability (CIN), suggesting that mimicking the methylation of Cse4-R37 throughout the cell cycle is detrimental to faithful chromosome segregation. Our results showed that SPOUT methyltransferase Upa1 contributes to methylation of Cse4-R37 and overexpression of UPA1 leads to CIN phenotype. In summary, our studies have defined a role for cell cycle-regulated methylation of Cse4 in high-fidelity chromosome segregation and highlight an important role of epigenetic modifications such as methylation of kinetochore proteins in preventing CIN, an important hallmark of human cancers.

Comparative genomics among cyst nematodes reveals distinct evolutionary histories among effector families and an irregular distribution of effector-associated promoter motifs.

Potato cyst nematodes (PCNs), an umbrella term used for two species, Globodera pallida and G. rostochiensis, belong worldwide to the most harmful pathogens of potato. Pathotype-specific host plant resistances are essential for PCN control. However, the poor delineation of G. pallida pathotypes has hampered the efficient use of available host plant resistances. Long-read sequencing technology allowed us to generate a new reference genome of G. pallida population D383 and, as compared to the current reference, the new genome assembly is 42 times less fragmented. For comparison of diversification patterns of six effector families between G. pallida and G. rostochiensis, an additional reference genome was generated for an outgroup, the beet cyst nematode Heterodera schachtii (IRS population). Large evolutionary contrasts in effector family topologies were observed. While VAPs (venom allergen-like proteins) diversified before the split between the three cyst nematode species, the families GLAND5 and GLAND13 only expanded in PCNs after their separation from the genus Heterodera. Although DNA motifs in the promoter regions thought to be involved in the orchestration of effector expression ("DOG boxes") were present in all three cyst nematode species, their presence is not a necessity for dorsal gland-produced effectors. Notably, DOG box dosage was only loosely correlated with the expression level of individual effector variants. Comparison of the G. pallida genome with those of two other cyst nematodes underlined the fundamental differences in evolutionary history between effector families. Resequencing of PCN populations with different virulence characteristics will allow for the linking of these characteristics to the composition of the effector repertoire as well as for the mapping of PCN diversification patterns resulting from extreme anthropogenic range expansion.

Comprehensive characterization of the Hsp70 interactome reveals novel client proteins and interactions mediated by posttranslational modifications.

Hsp70 interactions are critical for cellular viability and the response to stress. Previous attempts to characterize Hsp70 interactions have been limited by their transient nature and the inability of current technologies to distinguish direct versus bridged interactions. We report the novel use of cross-linking mass spectrometry (XL-MS) to comprehensively characterize the Saccharomyces cerevisiae (budding yeast) Hsp70 protein interactome. Using this approach, we have gained fundamental new insights into Hsp70 function, including definitive evidence of Hsp70 self-association as well as multipoint interaction with its client proteins. In addition to identifying a novel set of direct Hsp70 interactors that can be used to probe chaperone function in cells, we have also identified a suite of posttranslational modification (PTM)-associated Hsp70 interactions. The majority of these PTMs have not been previously reported and appear to be critical in the regulation of client protein function. These data indicate that one of the mechanisms by which PTMs contribute to protein function is by facilitating interaction with chaperones. Taken together, we propose that XL-MS analysis of chaperone complexes may be used as a unique way to identify biologically important PTMs on client proteins.

The genome and lifestage-specific transcriptomes of a plant-parasitic nematode and its host reveal susceptibility genes involved in trans-kingdom synthesis of vitamin B5.

Plant-parasitic nematodes are a major threat to crop production in all agricultural systems. The scarcity of classical resistance genes highlights a pressing need to find new ways to develop nematode-resistant germplasm. Here, we sequence and assemble a high-quality phased genome of the model cyst nematode Heterodera schachtii to provide a platform for the first system-wide dual analysis of host and parasite gene expression over time, covering all major parasitism stages. Analysis of the hologenome of the plant-nematode infection site identified metabolic pathways that were incomplete in the parasite but complemented by the host. Using a combination of bioinformatic, genetic, and biochemical approaches, we show that a highly atypical completion of vitamin B5 biosynthesis by the parasitic animal, putatively enabled by a horizontal gene transfer from a bacterium, is required for full pathogenicity. Knockout of either plant-encoded or now nematode-encoded steps in the pathway significantly reduces parasitic success. Our experiments establish a reference for cyst nematodes, further our understanding of the evolution of plant-parasitism by nematodes, and show that congruent differential expression of metabolic pathways in the infection hologenome represents a new way to find nematode susceptibility genes. The approach identifies genome-editing-amenable targets for future development of nematode-resistant crops.

Light-induced asymmetries in embryonic retinal gene expression are mediated by the vascular system and extracellular matrix.

Left-right asymmetries in the nervous system (lateralisation) influence a broad range of behaviours, from social responses to navigation and language. The role and pathways of endogenous and environmental mechanisms in the ontogeny of lateralisation remains to be established. The domestic chick is a model of both endogenous and experience-induced lateralisation driven by light exposure. Following the endogenous rightward rotation of the embryo, the asymmetrical position in the egg results in a greater exposure of the right eye to environmental light. To identify the genetic pathways activated by asymmetric light stimulation, and their time course, we exposed embryos to different light regimes: darkness, 6 h of light and 24 h of light. We used RNA-seq to compare gene expression in the right and left retinas and telencephalon. We detected differential gene expression in right vs left retina after 6 h of light exposure. This difference was absent in the darkness condition and had already disappeared by 24 h of light exposure, suggesting that light-induced activation is a self-terminating phenomenon. This transient effect of light exposure was associated with a downregulation of the sensitive-period mediator gene DIO2 (iodothyronine deiodinase 2) in the right retina. No differences between genes expressed in the right vs. left telencephalon were detected. Gene networks associated with lateralisation were connected to vascularisation, cell motility, and the extracellular matrix. Interestingly, we know that the extracellular matrix-including the differentially expressed PDGFRB gene-is involved in morphogenesis, sensitive periods, and in the endogenous chiral mechanism of primary cilia, that drives lateralisation. Our data show a similarity between endogenous and experience-driven lateralisation, identifying functional gene networks that affect lateralisation in a specific time window.

De Novo Assembly of Plasmodium knowlesi Genomes From Clinical Samples Explains the Counterintuitive Intrachromosomal Organization of Variant SICAvar and kir Multiple Gene Family Members.

Plasmodium knowlesi, a malaria parasite of Old World macaque monkeys, is used extensively to model Plasmodium biology. Recently, P. knowlesi was found in the human population of Southeast Asia, particularly Malaysia. P. knowlesi causes uncomplicated to severe and fatal malaria in the human host with features in common with the more prevalent and virulent malaria caused by Plasmodium falciparum. As such, P. knowlesi presents a unique opportunity to develop experimental translational model systems for malaria pathophysiology informed by clinical data from same-species human infections. Experimental lines of P. knowlesi represent well-characterized genetically stable parasites, and to maximize their utility as a backdrop for understanding malaria pathophysiology, genetically diverse contemporary clinical isolates, essentially wild-type, require comparable characterization. The Oxford Nanopore PCR-free long-read sequencing platform was used to sequence and de novo assemble P. knowlesi genomes from frozen clinical samples. The sequencing platform and assembly pipelines were designed to facilitate capturing data and describing, for the first time, P. knowlesi schizont-infected cell agglutination (SICA) var and Knowlesi-Interspersed Repeats (kir) multiple gene families in parasites acquired from nature. The SICAvar gene family members code for antigenically variant proteins analogous to the virulence-associated P. falciparum erythrocyte membrane protein (PfEMP1) multiple var gene family. Evidence presented here suggests that the SICAvar family members have arisen through a process of gene duplication, selection pressure, and variation. Highly evolving genes including PfEMP1family members tend to be restricted to relatively unstable sub-telomeric regions that drive change with core genes protected in genetically stable intrachromosomal locations. The comparable SICAvar and kir gene family members are counter-intuitively located across chromosomes. Here, we demonstrate that, in contrast to conserved core genes, SICAvar and kir genes occupy otherwise gene-sparse chromosomal locations that accommodate rapid evolution and change. The novel methods presented here offer the malaria research community not only new tools to generate comprehensive genome sequence data from small clinical samples but also new insight into the complexity of clinically important real-world parasites.

Tumour-associated antigenic peptides are present in the HLA class I ligandome of cancer cell line derived extracellular vesicles.

The recent success of monoclonal antibody checkpoint inhibitor therapies that enhance the ability of CD8+ T cells to detect cancer-related antigenic peptides has refocused the need to fully understand the repertoire of peptides being presented to the immune system. Whilst the peptide ligandome presented by cell surface human leucocyte antigen class I (HLA-I) molecules on cancer cells has been studied extensively, the ligandome of extracellular vesicles (EVs) remains poorly defined. Here, we report the HLA-I ligandome of both the cell surface and EVs from eight breast cancer cell lines (MCF7, MDA-MB-231, MDA-MB-361, MDA-MB-415, MDA-MB-453, HCC 1806, HCC 1395, and HCC 1954), and additionally the melanoma cell line ESTDAB-056 and the multiple myeloma line RPMI 8226. Utilizing HLA-I immunoisolation and mass spectrometry, we detected a total of 6574 peptides from the cell surface and 2461 peptides from the EVs of the cell lines studied. Within the EV HLA-I ligandome, we identified 150 peptides derived from tumour associated antigenic proteins, of which 19 peptides have been shown to elicit T-cell responses in previous studies. Our data thus show the prevalence of clinically relevant tumour-associated antigenic peptides in the HLA-I ligandome presented on EV.

ScreenGarden: a shinyR application for fast and easy analysis of plate-based high-throughput screens.

BACKGROUND: Colony growth on solid media is a simple and effective measure for high-throughput genomic experiments such as yeast two-hybrid, synthetic dosage lethality and Synthetic Physical Interaction screens. The development of robotic pinning tools has facilitated the experimental design of these assays, and different imaging software can be used to automatically measure colony sizes on plates. However, comparison to control plates and statistical data analysis is often laborious and pinning issues or plate specific growth effects can lead to the detection of false-positive growth defects. RESULTS: We have developed ScreenGarden, a shinyR application, to enable easy, quick and robust data analysis of plate-based high throughput assays. The code allows comparisons of different formats of data and different sized arrays of colonies. A comparison of ScreenGarden with previous analysis tools shows that it performs, at least, equivalently. The software can be run either via a website or offline via the RStudio program; the code is available and can be modified by expert uses to customise the analysis. CONCLUSIONS: ScreenGarden provides a simple, fast and effective tool to analyse colony growth data from genomic experiments.

Genome sequence of the English grain aphid, Sitobion avenae and its endosymbiont Buchnera aphidicola.

The English grain aphid, Sitobion avenae, is a major agricultural pest of wheat, barley and oats, and one of the principal vectors of barley yellow dwarf virus leading to significant reductions in grain yield, annually. Emerging resistance to and increasing regulation of insecticides has resulted in limited options for their control. Using PacBio HiFi data, we have produced a high-quality draft assembly of the S. avenae genome; generating a primary assembly with a total assembly size of 475.7 Mb, and an alternate assembly with a total assembly size of 430.8 Mb. Our primary assembly was highly contiguous with only 326 contigs and a contig N50 of 15.95 Mb. Assembly completeness was estimated at 97.7% using BUSCO analysis and 31,007 and 29,037 protein-coding genes were predicted from the primary and alternate assemblies, respectively. This assembly, which is to our knowledge the first for an insecticide resistant clonal lineage of English grain aphid, will provide novel insight into the molecular and mechanistic determinants of resistance and will facilitate future research into mechanisms of viral transmission and aphid behavior.

Forced association of SARS-CoV-2 proteins with the yeast proteome perturb vesicle trafficking.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the highly infectious coronavirus disease COVID-19. Extensive research has been performed in recent months to better understand how SARS-CoV-2 infects and manipulates its host to identify potential drug targets and support patient recovery from COVID-19. However, the function of many SARS-CoV-2 proteins remains uncharacterised. Here we used the Synthetic Physical Interactions (SPI) method to recruit SARS-CoV-2 proteins to most of the budding yeast proteome to identify conserved pathways which are affected by SARS-CoV-2 proteins. The set of yeast proteins that result in growth defects when associated with the viral proteins have homologous functions that overlap those identified in studies performed in mammalian cells. Specifically, we were able to show that recruiting the SARS-CoV-2 NSP1 protein to HOPS, a vesicle-docking complex, is sufficient to perturb membrane trafficking in yeast consistent with the hijacking of the endoplasmic-reticulum-Golgi intermediate compartment trafficking pathway during viral infection of mammalian cells. These data demonstrate that the yeast SPI method is a rapid way to identify potential functions of ectopic viral proteins.

Draft genome assemblies for tree pathogens Phytophthora pseudosyringae and Phytophthora boehmeriae.

Species of Phytophthora, plant pathogenic eukaryotic microbes, can cause disease on many tree species. Genome sequencing of species from this genus has helped to determine components of their pathogenicity arsenal. Here, we sequenced genomes for two widely distributed species, Phytophthora pseudosyringae and Phytophthora boehmeriae, yielding genome assemblies of 49 and 40 Mb, respectively. We identified more than 270 candidate disease promoting RXLR effector coding genes for each species, and hundreds of genes encoding candidate plant cell wall degrading carbohydrate active enzymes (CAZymes). These data boost genome sequence representation across the Phytophthora genus, and form resources for further study of Phytophthora pathogenesis.

Comparative genomics of two inbred lines of the potato cyst nematode Globodera rostochiensis reveals disparate effector family-specific diversification patterns.

BACKGROUND: Potato cyst nematodes belong to the most harmful pathogens in potato, and durable management of these parasites largely depends on host-plant resistances. These resistances are pathotype specific. The current Globodera rostochiensis pathotype scheme that defines five pathotypes (Ro1 - Ro5) is both fundamentally and practically of limited value. Hence, resistant potato varieties are used worldwide in a poorly informed manner. RESULTS: We generated two novel reference genomes of G. rostochiensis inbred lines derived from a Ro1 and a Ro5 population. These genome sequences comprise 173 and 189 scaffolds respectively, marking a ≈ 24-fold reduction in fragmentation as compared to the current reference genome. We provide copy number variations for 19 effector families. Four dorsal gland effector families were investigated in more detail. SPRYSECs, known to be implicated in plant defence suppression, constitute by far the most diversified family studied herein with 60 and 99 variants in Ro1 and Ro5 distributed over 18 and 26 scaffolds. In contrast, CLEs, effectors involved in feeding site induction, show strong physical clustering. The 10 and 16 variants cluster on respectively 2 and 1 scaffolds. Given that pathotypes are defined by their effectoromes, we pinpoint the disparate nature of the contributing effector families in terms of sequence diversification and loss and gain of variants. CONCLUSIONS: Two novel reference genomes allow for nearly complete inventories of effector diversification and physical organisation within and between pathotypes. Combined with insights we provide on effector family-specific diversification patterns, this constitutes a basis for an effectorome-based virulence scheme for this notorious pathogen.

The GpIA7 effector from the potato cyst nematode Globodera pallida targets potato EBP1 and interferes with the plant cell cycle.

The potato cyst nematode Globodera pallida acquires all of its nutrients from an elaborate feeding site that it establishes in a host plant root. Normal development of the root cells is re-programmed in a process coordinated by secreted nematode effector proteins. The biological function of the G. pallida GpIA7 effector was investigated in this study. GpIA7 is specifically expressed in the subventral pharyngeal glands of pre-parasitic stage nematodes. Ectopic expression of GpIA7 in potato plants affected plant growth and development, suggesting a potential role for this effector in feeding site establishment. Potato plants overexpressing GpIA7 were shorter, with reduced tuber weight and delayed flowering. We provide evidence that GpIA7 associates with the plant growth regulator StEBP1 (ErbB-3 epidermal growth factor receptor-binding protein 1). GpIA7 modulates the regulatory function of StEBP1, altering the expression level of downstream target genes, including ribonucleotide reductase 2, cyclin D3;1, and retinoblastoma related 1, which are down-regulated in plants overexpressing GpIA7. We provide an insight into the molecular mechanism used by the nematode to manipulate the host cell cycle and demonstrate that this may rely, at least in part, on hindering the function of host EBP1.

Cdc7-mediated phosphorylation of Cse4 regulates high-fidelity chromosome segregation in budding yeast.

Faithful chromosome segregation maintains chromosomal stability as errors in this process contribute to chromosomal instability (CIN), which has been observed in many diseases including cancer. Epigenetic regulation of kinetochore proteins such as Cse4 (CENP-A in humans) plays a critical role in high-fidelity chromosome segregation. Here we show that Cse4 is a substrate of evolutionarily conserved Cdc7 kinase, and that Cdc7-mediated phosphorylation of Cse4 prevents CIN. We determined that Cdc7 phosphorylates Cse4 in vitro and interacts with Cse4 in vivo in a cell cycle-dependent manner. Cdc7 is required for kinetochore integrity as reduced levels of CEN-associated Cse4, a faster exchange of Cse4 at the metaphase kinetochores, and defects in chromosome segregation, are observed in a cdc7-7 strain. Phosphorylation of Cse4 by Cdc7 is important for cell survival as constitutive association of a kinase-dead variant of Cdc7 (cdc7-kd) with Cse4 at the kinetochore leads to growth defects. Moreover, phospho-deficient mutations of Cse4 for consensus Cdc7 target sites contribute to CIN phenotype. In summary, our results have defined a role for Cdc7-mediated phosphorylation of Cse4 in faithful chromosome segregation.

The GpIA7 effector from the potato cyst nematode Globodera pallida targets potato EBP1 and interferes with the plant cell cycle programme.

The potato cyst nematode Globodera pallida acquires all of its nutrients from an elaborate feeding site that it establishes in a host plant root. Normal development of the root cells is re-programmed in a process coordinated by secreted nematode effector proteins. The biological function of the G. pallida GpIA7 effector was investigated in this study. GpIA7 is specifically expressed in the subventral pharyngeal glands of pre-parasitic stage nematodes. Ectopic expression of GpIA7 in potato plants affected plant growth and development, suggesting a potential role for this effector in feeding site establishment. Potato plants overexpressing GpIA7 were shorter, with reduced tuber weight and delayed flowering. We provide evidence that GpIA7 associates with the plant growth regulator StEBP1 (ErbB-3 epidermal growth factor receptor-binding protein 1). GpIA7 modulates the regulatory function of StEBP1, altering the expression level of downstream target genes, including ribonucleotide reductase 2, cyclin D3;1 and retinoblastoma related 1, which are downregulated in plants overexpressing GpIA7. We provide an insight into the molecular mechanism used by the nematode to manipulate the host cell cycle and provide evidence that this may rely, at least in part, on hindering the function of host EBP1.

Genome Expression Dynamics Reveal the Parasitism Regulatory Landscape of the Root-Knot Nematode Meloidogyne incognita and a Promoter Motif Associated with Effector Genes.

Root-knot nematodes (genus Meloidogyne) are the major contributor to crop losses caused by nematodes. These nematodes secrete effector proteins into the plant, derived from two sets of pharyngeal gland cells, to manipulate host physiology and immunity. Successful completion of the life cycle, involving successive molts from egg to adult, covers morphologically and functionally distinct stages and will require precise control of gene expression, including effector genes. The details of how root-knot nematodes regulate transcription remain sparse. Here, we report a life stage-specific transcriptome of Meloidogyne incognita. Combined with an available annotated genome, we explore the spatio-temporal regulation of gene expression. We reveal gene expression clusters and predicted functions that accompany the major developmental transitions. Focusing on effectors, we identify a putative cis-regulatory motif associated with expression in the dorsal glands, providing an insight into effector regulation. We combine the presence of this motif with several other criteria to predict a novel set of putative dorsal gland effectors. Finally, we show this motif, and thereby its utility, is broadly conserved across the Meloidogyne genus, and we name it Mel-DOG. Taken together, we provide the first genome-wide analysis of spatio-temporal gene expression in a root-knot nematode and identify a new set of candidate effector genes that will guide future functional analyses.