Transcriptome-wide selection of a reliable set of reference genes for gene expression studies in potato cyst nematodes (Globodera spp.)

Selection and Validation of Reference Genes for Gene Expression Studies in Euonymus japonicus Based on RNA Sequencing

Euonymus japonicus is one of the most low-temperature-tolerant evergreen broad-leaved tree species in the world and is widely used in urban greening. However, there are very few molecular biology studies on its low-temperature tolerance mechanism. So far, no researcher has selected and reported on its reference genes. In this study, 21 candidate reference genes (12 traditional housekeeping genes and 9 other genes) were initially selected based on gene expression and coefficient of variation (CV) through RNA-Seq (unpublished data), and qRT-PCR was used to detect the expression levels of candidate reference genes in three different groups of samples (leaves under different temperature stresses, leaves of plants at different growth stages, and different organs). After further evaluating the expression stability of these genes using geNorm, NormFinder, Bestkeeper, and RefFind, the results show that the traditional housekeeping gene eIF5A and the new reference gene RTNLB1 have good stability in the three different groups of samples, so they are reference genes with universality. In addition, we used eIF5A and RTNLB1 as reference genes to calibrate the expression pattern of the target gene EjMAH1, which confirmed this view. This article is the first to select and report on the reference gene of E. japonicus, laying the foundation for its low-temperature tolerance mechanism and other molecular biology research.

Single Nematode Transcriptomic Analysis, Using Long-Read Technology, Reveals Two Novel Virulence Gene Candidates in the Soybean Cyst Nematode, Heterodera glycines

The soybean cyst nematode (Heterodera glycines, SCN), is the most damaging disease of soybean in North America. While management of this pest using resistant soybean is generally still effective, prolonged exposure to cultivars derived from the same source of resistance (PI 88788) has led to the emergence of virulence. Currently, the underlying mechanisms responsible for resistance breakdown remain unknown. In this study, we combined a single nematode transcriptomic profiling approach with long-read sequencing to reannotate the SCN genome. This resulted in the annotation of 1932 novel transcripts and 281 novel gene features. Using a transcript-level quantification approach, we identified eight novel effector candidates overexpressed in PI 88788 virulent nematodes in the late infection stage. Among these were the novel gene Hg-CPZ-1 and a pioneer effector transcript generated through the alternative splicing of the non-effector gene Hetgly21698. While our results demonstrate that alternative splicing in effectors does occur, we found limited evidence of direct involvement in the breakdown of resistance. However, our analysis highlighted a distinct pattern of effector upregulation in response to PI 88788 resistance indicative of a possible adaptation process by SCN to host resistance.

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.

Identification and validation of reference genes for real-time RT-PCR in Aphelenchoides besseyi

Fragments of four candidate reference genes of Aphelenchoides besseyi, including actin, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ubiquitin conjugating-3 enzyme (UBC) and alpha-tubulin (α-tubulin) were cloned from the transcriptome database of A. besseyi. The expression level of these four candidate reference genes and a commonly used reference gene of A. besseyi (18S rRNA) in three experimental conditions, including the four life stages (female, male, juvenile and egg) of two populations and the mixed-stage nematodes of four populations with different origins and hosts were analyzed by RT-qPCR. The expression stability of the five candidate reference genes under the three experimental conditions was analyzed by ΔCt, geNorm, NormFinder and RefFinder respectively. The analysis results of ΔCt, geNorm, NormFinder and RefFinder all indicated that UBC was the gene with the highest average ranking of stability. In conclusion, the expression stability of UBC was optimal under the three experimental conditions, indicating that UBC could be used as a suitable reference gene instead of 18S rRNA in the RT-qPCR analysis for A. besseyi.

Selection and Validation of Reference Genes for Gene Expression Studies in Codonopsis pilosula Based on Transcriptome Sequence Data

Relative gene expression analyses by RT-qPCR (reverse transcription-quantitative PCR) are highly dependent on the reference genes in normalizing the expression data of target genes. Therefore, inappropriate endogenous control genes will lead to inaccurate target gene expression profiles, and the selection and validation of suitable internal reference genes becomes essential. In this study, we retrieved the commonly used reference genes in transcriptome datasets of Codonopsis pilosula by RNA-Seq (unpublished data), and selected 15 candidate reference genes according to the coefficient of variation (CV) and fold change (FC) value of gene expression. The expression levels of candidate reference genes, which is at different growth stages, undergoing cold stress and drought stress, was determined by RT-qPCR. The expression stability of these genes was evaluated using software packages and algorithms including ΔCt, geNorm, NormFinder and Bestkeeper. Then appropriate reference genes were screened and validated by target gene-UDGPase (UDP glucose pyrophosphorylase). The optimal RGs combinations of C. pilosula, including PP2A59γ, CPY20-1, UBCE32, RPL5B and UBC18 for developmental stage, RPL5B, RPL13 and PP2A59γ for cold treatment, RPL13 and PP2A59γ for drought treatment, were found and proposed as reference genes for future work. This paper laid foundations for both the selection of reference genes and exploration in metabolic mechanism of C. pilosula.

Evaluation of duplicated reference genes for quantitative real-time PCR analysis in genome unknown hexaploid oat (Avena sativa L.)

BACKGROUND

Oat (Avena sativa L.), a hexaploid crop with unknown genome, has valuable nutritional, medicinal and pharmaceutical uses. However, no suitable RGs (reference genes) for qPCR (quantitative real-time PCR) has been documented for oat yet. Single-copy gene is often selected as RG, which is challengeable or impactable in unexplored polyploids.

RESULTS

In this study, eleven candidate RGs, including four duplicated genes, were selected from oat transcriptome. The stability and the optimal combination of these candidate RGs were assessed in 18 oat samples by using four statistical algorithms including the ΔCt method, geNorm, NormFinder and BestKeeper. The most stable RGs for "all samples", "shoots and roots of seedlings", "developing seeds" and "developing endosperms" were EIF4A (Eukaryotic initiation factor 4A-3), UBC21 (Ubiquitin-Conjugating Enzyme 21), EP (Expressed protein) and EIF4A respectively. Among these RGs, UBC21 was a four-copy duplicated gene. The reliability was validated by the expression patterns of four various genes normalized to the most and the least stable RGs in different sample sets.

CONCLUSIONS

Results provide a proof of concept that the duplicated RG is feasible for qPCR in polyploids. To our knowledge, this study is the first systematic research on the optimal RGs for accurate qPCR normalization of gene expression in different organs and tissues of oat.

"Piler Dirt" Survey for the Sampling and Detection of Potato Cyst Nematodes

Potato cyst nematodes are a significant threat to potato production worldwide and have important economic impacts due to yield losses but also because of the expenses associated with regulation procedures. In order to reduce the sampling labor, an alternative strategy named the "Piler Dirt" that collects the soil carried with potato tubers during their transfer to storage was proposed. The method showed a better sensitivity than the reference method to detect fields infested with G. rostochiensis. The quantification of the number of cysts per kilogram of soil was proportional between the two methods at low and moderate population densities (R² = 0.885) but no correlations were found at high density. However, the quantity of soil generated by the method was exceedingly large to be treated by diagnostic labs. It was shown that subsampling six aliquots, each equivalent to 5,000 cm³/ha, from the total quantity of soil generated by the Piler Dirt method, resulted in a probability of 97% to detect infested fields, 95% of the time in our dataset. Overall, Piler Dirt appears as a good compromise to reduce labor time and cost without significantly affecting sensitivity. However, it will be challenging to implement because it needs to be done simultaneously with harvest and will require the participation of farmers during a busy period.

What determines host specificity in hyperspecialized plant parasitic nematodes?

Background

In hyperspecialized parasites, the ability to grow on a particular host relies on specific virulence factors called effectors. These excreted proteins are involved in the molecular mechanisms of parasitism and distinguish virulent pathogens from non-virulent related species. The potato cyst nematodes (PCN) Globodera rostochiensis and G. pallida are major plant-parasitic nematodes developing on numerous solanaceous species including potato. Their close relatives, G. tabacum and G. mexicana are stimulated by potato root diffusate but unable to establish a feeding site on this plant host.

Results

RNA sequencing was used to characterize transcriptomic differences among these four Globodera species and to identify genes associated with host specificity. We identified seven transcripts that were unique to PCN species, including a protein involved in ubiquitination. We also found 545 genes that were differentially expressed between PCN and non-PCN species, including 78 genes coding for effector proteins, which represent more than a 6-fold enrichment compared to the whole transcriptome. Gene polymorphism analysis identified 359 homozygous non-synonymous variants showing a strong evidence for selection in PCN species.

Conclusions

Overall, we demonstrated that the determinant of host specificity resides in the regulation of essential effector gene expression that could be under the control of a single or of very few regulatory genes. Such genes are therefore promising targets for the development of novel and more sustainable resistances against potato cyst nematodes.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5853-4) contains supplementary material, which is available to authorized users.