Multi-Omic Investigation of Low-Nitrogen Conditional Resistance to Clubroot Reveals Brassica napus Genes Involved in Nitrate Assimilation.

Nitrogen fertilization has been reported to influence the development of clubroot, a root disease of Brassicaceae species, caused by the obligate protist Plasmodiophora brassicae. Our previous works highlighted that low-nitrogen fertilization induced a strong reduction of clubroot symptoms in some oilseed rape genotypes. To further understand the underlying mechanisms, the response to P. brassicae infection was investigated in two genotypes "Yudal" and HD018 harboring sharply contrasted nitrogen-driven modulation of resistance toward P. brassicae. Targeted hormone and metabolic profiling, as well as RNA-seq analysis, were performed in inoculated and non-inoculated roots at 14 and 27 days post-inoculation, under high and low-nitrogen conditions. Clubroot infection triggered a large increase of SA concentration and an induction of the SA gene markers expression whatever the genotype and nitrogen conditions. Overall, metabolic profiles suggested that N-driven induction of resistance was independent of SA signaling, soluble carbohydrate and amino acid concentrations. Low-nitrogen-driven resistance in "Yudal" was associated with the transcriptional regulation of a small set of genes, among which the induction of NRT2- and NR-encoding genes. Altogether, our results indicate a possible role of nitrate transporters and auxin signaling in the crosstalk between plant nutrition and partial resistance to pathogens.

Nitrogen Supply and Host-Plant Genotype Modulate the Transcriptomic Profile of Plasmodiophora brassicae.

Nitrogen fertilization can affect the susceptibility of Brassica napus to the telluric pathogen Plasmodiophora brassicae. Our previous works highlighted that the influence of nitrogen can strongly vary regarding plant cultivar/pathogen strain combinations, but the underlying mechanisms are unknown. The present work aims to explore how nitrogen supply can affect the molecular physiology of P. brassicae through its life epidemiological cycle. A time-course transcriptome experiment was conducted to study the interaction, under two conditions of nitrogen supply, between isolate eH and two B. napus genotypes (Yudal and HD-018), harboring (or not harboring) low nitrogen-conditional resistance toward this isolate (respectively). P. brassicae transcriptional patterns were modulated by nitrogen supply, these modulations being dependent on both host-plant genotype and kinetic time. Functional analysis allowed the identification of P. brassicae genes expressed during the secondary phase of infection, which may play a role in the reduction of Yudal disease symptoms in low-nitrogen conditions. Candidate genes included pathogenicity-related genes ("NUDIX," "carboxypeptidase," and "NEP-proteins") and genes associated to obligate biotrophic functions of P. brassicae. This work illustrates the importance of considering pathogen's physiological responses to get a better understanding of the influence of abiotic factors on clubroot resistance/susceptibility.

Soil microbiota influences clubroot disease by modulating Plasmodiophora brassicae and Brassica napus transcriptomes.

The contribution of surrounding plant microbiota to disease development has led to the 'pathobiome' concept, which represents the interaction between the pathogen, the host plant and the associated biotic microbial community, resulting or not in plant disease. The aim herein is to understand how the soil microbial environment may influence the functions of a pathogen and its pathogenesis, and the molecular response of the plant to the infection, with a dual-RNAseq transcriptomics approach. We address this question using Brassica napus and Plasmodiophora brassicae, the pathogen responsible for clubroot. A time-course experiment was conducted to study interactions between P. brassicae, two B. napus genotypes and three soils harbouring high, medium or low microbiota diversities and levels of richness. The soil microbial diversity levels had an impact on disease development (symptom levels and pathogen quantity). The P. brassicae and B. napus transcriptional patterns were modulated by these microbial diversities, these modulations being dependent on the host genotype plant and the kinetic time. The functional analysis of gene expressions allowed the identification of pathogen and plant host functions potentially involved in the change of plant disease level, such as pathogenicity-related genes (NUDIX effector) in P. brassicae and plant defence-related genes (glucosinolate metabolism) in B. napus.

pH effect on strain-specific transcriptomes of the take-all fungus.

The soilborne fungus Gaeumannomyces tritici (G. tritici) causes the take-all disease on wheat roots. Ambient pH has been shown to be critical in different steps of G. tritici life cycle such as survival in bulk soil, saprophytic growth, and pathogenicity on plants. There are however intra-specific variations and we previously found two types of G. tritici strains that grow preferentially either at acidic pH or at neutral/alkaline pH; gene expression involved in pH-signal transduction pathway and pathogenesis was differentially regulated in two strains representative of these types. To go deeper in the description of the genetic pathways and the understanding of this adaptative mechanism, transcriptome sequencing was achieved on two strains (PG6 and PG38) which displayed opposite growth profiles in two pH conditions (acidic and neutral). PG6, growing better at acidic pH, overexpressed in this condition genes related to cell proliferation. In contrast, PG38, which grew better at neutral pH, overexpressed in this condition genes involved in fatty acids and amino acid metabolisms, and genes potentially related to pathogenesis. This strain also expressed stress resistance mechanisms at both pH, to assert a convenient growth under various ambient pH conditions. These differences in metabolic pathway expression between strains at different pH might buffer the effect of field or soil variation in wheat fields, and explain the success of the pathogen.

Temporal dynamics of bacterial and fungal communities during the infection of Brassica rapa roots by the protist Plasmodiophora brassicae.

The temporal dynamics of rhizosphere and root microbiota composition was compared between healthy and infected Chinese cabbage plants by the pathogen Plasmodiophora brassicae. When inoculated with P. brassicae, disease was measured at five sampling dates from early root hair infection to late gall development. The first symptoms of clubroot disease appeared 14 days after inoculation (DAI) and increased drastically between 14 and 35 DAI. The structure of microbial communities associated to rhizosphere soil and root from healthy and inoculated plants was characterized through high-throughput DNA sequencing of bacterial (16S) and fungal (18S) molecular markers and compared at each sampling date. In healthy plants, Proteobacteria and Bacteroidetes bacterial phyla dominated the rhizosphere and root microbiota of Chinese cabbage. Rhizosphere bacterial communities contained higher abundances of Actinobacteria and Firmicutes compared to the roots. Moreover, a drastic shift of fungal communities of healthy plants occurred between the two last sampling dates, especially in plant roots, where most of Ascomycota fungi dominated until they were replaced by a fungus assigned to the Chytridiomycota phylum. Parasitic invasion by P. brassicae disrupted the rhizosphere and root-associated community assembly at a late step during the root secondary cortical infection stage of clubroot disease. At this stage, Flavisolibacter and Streptomyces in the rhizosphere, and Bacillus in the roots, were drastically less abundant upon parasite invasion. Rhizosphere of plants colonized by P. brassicae was significantly more invaded by the Chytridiomycota fungus, which could reflect a mutualistic relationship in this compartment between these two microorganisms.

Computational analysis of the Plasmodiophora brassicae genome: mitochondrial sequence description and metabolic pathway database design.

Plasmodiophora brassicae is an obligate biotrophic pathogenic protist responsible for clubroot, a root gall disease of Brassicaceae species. In addition to the reference genome of the P. brassicae European e3 isolate and the draft genomes of Canadian or Chinese isolates, we present the genome of eH, a second European isolate. Refinement of the annotation of the eH genome led to the identification of the mitochondrial genome sequence, which was found to be bigger than that of Spongospora subterranea, another plant parasitic Plasmodiophorid phylogenetically related to P. brassicae. New pathways were also predicted, such as those for the synthesis of spermidine, a polyamine up-regulated in clubbed regions of roots. A P. brassicae pathway genome database was created to facilitate the functional study of metabolic pathways in transcriptomics approaches. These available tools can help in our understanding of the regulation of P. brassicae metabolism during infection and in response to diverse constraints.

Strain-specific variation in a soilborne phytopathogenic fungus for the expression of genes involved in pH signal transduction pathway, pathogenesis and saprophytic survival in response to environmental pH changes.

The soilborne fungus Gaeumannomyces graminis var. tritici (Ggt) causes take-all, a wheat root disease. In an original strain-specific way, a previous study indicates that inside the Ggt species, some strains grow preferentially at acidic pH and other strains at neutral/alkaline pH. The most important mechanism for a fungal response to the environmental pH is the Pal pathway which integrates the products of the six pal genes and the transcription factor PacC. To evaluate whether the Ggt strain-specific growth in function of the ambient pH is mediated via the Pal pathway, a transcriptional study of the genes encoding this pathway was carried out. This study provided the first evidence that the pH signalling pathway similar to those described in other fungi operated in Ggt. The pacC gene was induced at neutral pH whatever the strain. In an original way, the expression of Ggt genes coding for the different Pal proteins depended on the strain and on the ambient pH. In the strain growing better at acidic pH, few pal genes were pH-regulated, and some were overexpressed at neutral pH when regulated. In the strain growing better at neutral pH, underexpression of most of the pal genes at neutral pH occurred. The strains displayed higher gene expression in the ambient pH that unfavoured their growth as if it was a compensation system. All pH taken together, a globally weaker Pal transcript level occurred in the strains that were less sensitive to acidic pH, and on the contrary, the strain growing better on neutral pH showed higher Pal mRNA levels. The expression of genes involved in pathogenesis and saprophytic growth was also regulated by the ambient pH and the strain: each gene displayed a specific pH-regulation that was similar between strains. But all pH taken together, the global transcript levels of four out of six genes were higher in the strain growing better on neutral pH. Altogether, for the first time, the results show that inside a species, conditions affecting environmental pH modulate the expression of genes in an original strain-specific way.

Genomic analysis of the biocontrol strain Pseudomonas fluorescens Pf29Arp with evidence of T3SS and T6SS gene expression on plant roots.

Several bacterial strains of the Pseudomonas genus provide plant growth stimulation, plant protection against pests or bioremediation. Among these bacteria, P. fluorescens Pf29Arp reduces the severity of take-all, a disease caused by the pathogenic fungus Gaeumannomyces graminis var. tritici (Ggt) on wheat roots. In this study, we obtained a draft genome of Pf29Arp and subsequent comparative genomic analyses have revealed that this bacterial strain is closely related to strains of the 'P. brassicacearum-like' subgroup including P. brassicacearum ssp. brassicacearum NFM421 and P. fluorescens F113. Despite an overall chromosomal organization similar to these strains, a number of features including antibiotic synthesis gene clusters from secondary metabolism are not found in the Pf29Arp genome. But Pf29Arp possesses different protein secretion systems including type III (T3SS) and type VI (T6SS) secretion systems. Pf29Arp is the first Pseudomonas sp. strain described with four T6SS clusters (cluster I, II, III and IV). In addition, some protein-coding genes involved in the assembly of these secretion systems are basally expressed during Pf29Arp colonization of healthy wheat roots and display different expression patterns on necrotized roots caused by Ggt. These data suggest a role of T3SS and T6SS in the Pf29Arp adaptation to different root environments.

The biocontrol bacterium Pseudomonas fluorescens Pf29Arp strain affects the pathogenesis-related gene expression of the take-all fungus Gaeumannomyces graminis var. tritici on wheat roots.

The main effects of antagonistic rhizobacteria on plant pathogenic fungi are antibiosis, fungistasis or an indirect constraint through the induction of a plant defence response. To explore different biocontrol mechanisms, an in vitro confrontation assay was conducted with the rhizobacterium Pseudomonas fluorescens Pf29Arp as a biocontrol agent of the fungus Gaeumannomyces graminis var. tritici (Ggt) on wheat roots. In parallel with the assessment of disease extension, together with the bacterial and fungal root colonization rates, the transcript levels of candidate fungal pathogenicity and plant-induced genes were monitored during the 10-day infection process. The bacterial inoculation of wheat roots with the Pf29Arp strain reduced the development of Ggt-induced disease expressed as attack frequency and necrosis length. The growth rates of Ggt and Pf29Arp, monitored through quantitative polymerase chain reaction of DNA amounts with a part of the Ggt 18S rDNA gene and a specific Pf29Arp strain detection probe, respectively, increased throughout the interactions. Bacterial antagonism and colonization had no significant effect on root colonization by Ggt. The expression of fungal and plant genes was quantified in planta by quantitative reverse transcription-polymerase chain reaction during the interactions thanks to the design of specific primers and an innovative universal reference system. During the early stages of the tripartite interaction, several of the fungal genes assayed were down-regulated by Pf29Arp, including two laccases, a ß-1,3-exoglucanase and a mitogen-activated protein kinase. The plant host glutathione-S-transferase gene was induced by Ggt alone and up-regulated by Pf29Arp bacteria in interaction with the pathogen. We conclude that Pf29Arp antagonism acts through the alteration of fungal pathogenesis and probably through the activation of host defences.

Several dystrophin-glycoprotein complex members are present in crude surface membranes but they are sodium dodecyl sulphate invisible in KCl-washed microsomes from mdx mouse muscle.

The dystrophin-glycoprotein complex (DGC) is a large trans-sarcolemmal complex that provides a linkage between the subsarcolemmal cytoskeleton and the extracellular matrix. In skeletal muscle, it consists of the dystroglycan, sarcoglycan and cytoplasmic complexes, with dystrophin forming the core protein. The DGC has been described as being absent or greatly reduced in dystrophin-deficient muscles, and this lack is considered to be involved in the dystrophic phenotype. Such a decrease in the DGC content was observed in dystrophin-deficient muscle from humans with muscular dystrophy and in mice with X-linked muscular dystrophy (mdx mice). These deficits were observed in total muscle homogenates and in partially membrane-purified muscle fractions, the so-called KCl-washed microsomes. Here, we report that most of the proteins of the DGC are actually present at normal levels in the mdx mouse muscle plasma membrane. The proteins are detected in dystrophic animal muscles when the immunoblot assay is performed with crude surface membrane fractions instead of the usually employed KCl-washed microsomes. We propose that these proteins form SDS-insoluble membrane complexes when dystrophin is absent.

The fatty acid desaturase 3 gene encodes for different FADS3 protein isoforms in mammalian tissues.

In 2000, Marquardt et al. (A. Marquardt, H. Stöhr, K. White, and B. H. F. Weber. 2000. cDNA cloning, genomic structure, and chromosomal localization of three members of the human fatty acid desaturase family. Genomics. 66: 176-183.) described the genomic structure of the fatty acid desaturase (FADS) cluster in humans. This cluster includes the FADS1 and FADS2 genes encoding, respectively, for the Delta 5- and Delta 6-desaturases involved in polyunsaturated fatty acid biosynthesis. A third gene, named FADS3, has recently been identified but no functional role has yet been attributed to the putative FADS3 protein. In this study, we investigated the FADS3 occurrence in rat tissues by using two specific polyclonal antibodies directed against the N-terminal and C-terminal ends of rat FADS3. Our results showed three potential protein isoforms of FADS3 (75 kDa, 51 kDa, and 37 kDa) present in a tissue-dependent manner. The occurrence of these FADS3 isoforms did not depend on the mRNA level determined by real-time PCR. In parallel, mouse tissues were also tested and showed the same three FADS3 isoforms but with a different tissue distribution. Finally, we reported the existence of FADS3 in human cells and tissues but different new isoforms were identified. To conclude, we showed in this study that FADS3 does exist under multiple protein isoforms depending on the mammalian tissues. These results will help further investigations to determine the physiological function of FADS3.

Identification and characterization of recombinant and native rat myristoyl-CoA: protein N-myristoyltransferases.

Compared to other species that possess a single functional myristoyl-CoA: protein N-myristoyltransferase gene copy, human, mouse and cow possess 2 NMT genes, and more than 2 protein isoforms. In mammals, the contribution of each gene transcript to multiple protein isoform expression and enzyme activity remains unclear. In order to get new insight on their respective physiological role, we have cloned and characterized the two rat NMT cDNAs. Rat NMT1 and NMT2 cDNAs contain 1491 and 1590 nucleotides, respectively, with high identity with their mouse homologues. Polypeptide sequences exhibited 68.1% identity between NMT1 and 2. Recombinant rat NMT1 and 2 showed major immunoreactive forms at 66 and 50 kDa, although NMT2 is 33-amino acid longer than NMT1. Both proteins exhibited functional myristoyltransferase activity but NMT2 appeared to be 4-time less active than NMT1. Studies of native protein expression revealed that the level and sizes of NMT proteins greatly vary among rat tissues although NMT1 and 2 did not display tissue specific expression at the mRNA level. Altogether, these results suggest that NMT2 may contribute little to total NMT activity levels in vivo.

Myristic acid increases delta6-desaturase activity in cultured rat hepatocytes.

In order to study the effects of saturated fatty acids on delta6-desaturase activity, rat hepatocytes in primary culture were incubated with lauric (C12:0), myristic (C14:0) or palmitic (C16:0) acids. After optimization, the standard in vitro conditions for the measurement of delta6-desaturase activity were as follows: 60 micromol x L(-1) alpha-linolenic acid (C18:3n-3), reaction time of 20 min and protein content of 0.4 mg. Data showed that cell treatment with 0.5 mmol x L(-1) myristic acid during 43 h specifically increased delta6-desaturase activity. This improvement, reproducible for three substrates of delta6-desaturase, i.e. oleic acid (C18:1n-9), linoleic acid (C18:2n-6) and alpha-linoleic acid (C18:3n-3) was dose-dependent in the range 0.1-0.5 mmol x L(-1) myristic acid concentration.

Genetic linkage and expression analysis of SREBP and lipogenic genes in fat and lean chicken.

To identify the genes directly responsible, through DNA polymorphism, for the difference in fatness observed between a lean and a fat chicken line, we studied five genes (ACL, ACC, FAS, ME, SCD1) encoding key enzymes involved in liver fatty acid synthesis and secretion. Genetic linkage was tested between polymorphic sites in the genes and the fatness trait segregating in an F2 design obtained by inter-crossing the two fat and lean lines. Despite a confirmation of a higher mRNA level in the fat birds, no genetic linkage of the gene alleles with the phenotype could be found. As a test of the implication of upstream regulatory transcription factors, SREPB genes were also studied. The lack of genetic linkage of SREBP genes with fatness shows that these genes are not directly responsible through polymorphism for fatness variability in our model. Moreover, the similar SREBP mRNA levels observed between the two lines led us to exclude also transcriptional factors regulating the two SREBP genes as being directly responsible for fatness variability. However, the genes involved in post-translational modifications of SREBPs remain candidates to investigate. These results emphasised the interest to perform expression and genetic linkage studies jointly, to progress in identifying the genetic origin of variability of a quantitative trait.

Although it is rapidly metabolized in cultured rat hepatocytes, lauric acid is used for protein acylation.

This study was designed to examine the metabolic fate of exogenous lauric acid in cultured rat hepatocytes, in terms of both lipid metabolism and acylation of proteins. Radiolabeled [14C]-lauric acid at 0.1 mM in the culture medium was rapidly taken up by the cells (94.8 +/- 2.2% of the initial radioactivity was cleared from the medium after a 4 h incubation) but its incorporation into cellular lipids was low (24.6 +/- 4.2% of initial radioactivity after 4 h), due to the high beta-oxidation of lauric acid in hepatocytes (38.7 +/- 4.4% after the same time). Among cellular lipids, lauric acid was preferentially incorporated into triglycerides (10.6 +/- 4.6% of initial radioactivity after 4 h). Lauric acid was also rapidly converted to palmitic acid by two successive elongations. Protein acylation was detected after metabolic labeling of the cells with [11,12-3H]-lauric acid. Two-dimensional electrophoresis separation of the cellular proteins and autoradiography evidenced the incorporation of radioactivity into 35 well-resolved proteins. Radiolabeling of several proteins resulted from covalent linkage to the precursor [11,12-3H]-lauric acid or to its elongation product, myristic acid. The covalent linkages between these proteins and lauric acid were broken by base hydrolysis, indicating that the linkage was of the thioester or ester-type. Endogenous myristic acid produced by lauric acid elongation was used for both protein N-myristoylation and protein S-acylation. Therefore, these results show for the first time that, although it is rapidly metabolized in hepatocytes, exogenous lauric acid is a substrate for the acylation of liver proteins.