Development of a novel target-based cell assay, reporter of the activity of Mycobacterium tuberculosis protein-O-mannosyltransferase.

The Protein-O-mannosyltransferase is crucial for the virulence of Mycobacterium tuberculosis, the etiological agent of tuberculosis. This enzyme, called MtPMT (Rv1002c), is responsible for the post-translational O-mannosylation of mycobacterial proteins. It catalyzes the transfer of a single mannose residue from a polyprenol phospho-mannosyl lipidic donor to the hydroxyl groups of selected Ser/Thr residues in acceptor proteins during their translocation across the membrane. Previously, we provided evidence that the loss of MtPMT activity causes the absence of mannoproteins in Mycobacterium tuberculosis, severely impacting its intracellular growth, as well as a strong attenuation of its pathogenicity in immunocompromised mice. Therefore, it is of interest to develop specific inhibitors of this enzyme to better understand mycobacterial infectious diseases. Here we report the development of a "target-based" phenotypic assay for this enzyme, assessing its O-mannosyltransferase activity in bacteria, in the non-pathogenic Mycobacterium smegmatis strain. Robustness of the quantitative contribution of this assay was evaluated by intact protein mass spectrometry, using a panel of control strains, overexpressing the MtPMT gene, carrying different key point-mutations. Then, screening of a limited library of 30 compounds rationally chosen allowed us to identify 2 compounds containing pyrrole analogous rings, as significant inhibitors of MtPMT activity, affecting neither the growth of the mycobacterium nor its secretion of mannoproteins. These molecular cores could therefore serve as scaffold for the design of new pharmaceutical agents that could improve treatment of mycobacterial diseases. We report here the implementation of a miniaturized phenotypic activity assay for a glycosyltransferase of the C superfamily.

Natural mutations in the sensor kinase of the PhoPR two-component regulatory system modulate virulence of ancestor-like tuberculosis bacilli.

The molecular factors and genetic adaptations that contributed to the emergence of Mycobacterium tuberculosis (MTB) from an environmental Mycobacterium canettii-like ancestor, remain poorly investigated. In MTB, the PhoPR two-component regulatory system controls production and secretion of proteins and lipid virulence effectors. Here, we describe that several mutations, present in phoR of M. canettii relative to MTB, impact the expression of the PhoP regulon and the pathogenicity of the strains. First, we establish a molecular model of PhoR and show that some substitutions found in PhoR of M. canettii are likely to impact the structure and activity of this protein. Second, we show that STB-K, the most attenuated available M. canettii strain, displays lower expression of PhoP-induced genes than MTB. Third, we demonstrate that genetic swapping of the phoPR allele from STB-K with the ortholog from MTB H37Rv enhances expression of PhoP-controlled functions and the capacities of the recombinant strain to colonize human macrophages, the MTB target cells, as well as to cause disease in several mouse infection models. Fourth, we extended these observations to other M. canettii strains and confirm that PhoP-controlled functions are expressed at lower levels in most M. canettii strains than in M. tuberculosis. Our findings suggest that distinct PhoR variants have been selected during the evolution of tuberculosis bacilli, contributing to higher pathogenicity and persistence of MTB in the mammalian host.

A Mycobacterium tuberculosis fingerprint in human breath allows tuberculosis detection.

An estimated one-third of tuberculosis (TB) cases go undiagnosed or unreported. Sputum samples, widely used for TB diagnosis, are inefficient at detecting infection in children and paucibacillary patients. Indeed, developing point-of-care biomarker-based diagnostics that are not sputum-based is a major priority for the WHO. Here, in a proof-of-concept study, we tested whether pulmonary TB can be detected by analyzing patient exhaled breath condensate (EBC) samples. We find that the presence of Mycobacterium tuberculosis (Mtb)-specific lipids, lipoarabinomannan lipoglycan, and proteins in EBCs can efficiently differentiate baseline TB patients from controls. We used EBCs to track the longitudinal effects of antibiotic treatment in pediatric TB patients. In addition, Mtb lipoarabinomannan and lipids were structurally distinct in EBCs compared to ex vivo cultured bacteria, revealing specific metabolic and biochemical states of Mtb in the human lung. This provides essential information for the rational development or improvement of diagnostic antibodies, vaccines and therapeutic drugs. Our data collectively indicate that EBC analysis can potentially facilitate clinical diagnosis of TB across patient populations and monitor treatment efficacy. This affordable, rapid and non-invasive approach seems superior to sputum assays and has the potential to be implemented at point-of-care.

Robust transcriptional indicators of immune cell death revealed by spatiotemporal transcriptome analyses.

Recognition of a pathogen by the plant immune system often triggers a form of regulated cell death traditionally known as the hypersensitive response (HR). This type of cell death occurs precisely at the site of pathogen recognition, and it is restricted to a few cells. Extensive research has shed light on how plant immune receptors are mechanistically activated. However, two central key questions remain largely unresolved: how does cell death zonation take place, and what are the mechanisms that underpin this phenomenon? Consequently, bona fide transcriptional indicators of HR are lacking, which prevents deeper insight into its mechanisms before cell death becomes macroscopic and precludes early or live observation. In this study, to identify the transcriptional indicators of HR we used the paradigmatic Arabidopsis thaliana-Pseudomonas syringae pathosystem and performed a spatiotemporally resolved gene expression analysis that compared infected cells that will undergo HR upon pathogen recognition with bystander cells that will stay alive and activate immunity. Our data revealed unique and time-dependent differences in the repertoire of differentially expressed genes, expression profiles, and biological processes derived from tissue undergoing HR and that of its surroundings. Furthermore, we generated a pipeline based on concatenated pairwise comparisons between time, zone, and treatment that enabled us to define 13 robust transcriptional HR markers. Among these genes, the promoter of an uncharacterized AAA-ATPase was used to obtain a fluorescent reporter transgenic line that displays a strong spatiotemporally resolved signal specifically in cells that will later undergo pathogen-triggered cell death. This valuable set of genes can be used to define cells that are destined to die upon infection with HR-triggering bacteria, opening new avenues for specific and/or high-throughput techniques to study HR processes at a single-cell level.

Bioinformatic Mining and Structure-Activity Profiling of Baeyer-Villiger Monooxygenases from Mycobacterium tuberculosis.

Mycobacterium tuberculosis is the etiological agent of tuberculosis (TB), one of the deadliest infectious diseases. The alarming health context coupled with the emergence of resistant M. tuberculosis strains highlights the urgent need to expand the range of anti-TB antibiotics. A subset of anti-TB drugs in use are prodrugs that require bioactivation by a class of M. tuberculosis enzymes called Baeyer-Villiger monooxygenases (BVMOs), which remain understudied. To examine the prevalence and the molecular function of BVMOs in mycobacteria, we applied a comprehensive bioinformatic analysis that identified six BVMOs in M. tuberculosis, including Rv3083 (MymA), Rv3854c (EthA), Rv0565c, and Rv0892, which were selected for further characterization. Homology modeling and substrate docking analysis, performed on this subset, suggested that Rv0892 is closer to the cyclohexanone BVMO, while Rv0565c and EthA are structurally and functionally similar to MymA, which is by far the most prominent type I BVMO enzyme. Thanks to an unprecedented purification and assay optimization, biochemical studies confirmed that all four BVMOs display BV-oxygenation activity. We also showed that MymA displays a distinctive substrate preference that we further investigated by kinetic parameter determination and that correlates with in silico modeling. We provide insights into distribution of BVMOs and the structural basis of their substrate profiling, and we discuss their possible redundancy in M. tuberculosis, raising questions about their versatility in prodrug activation and their role in physiology and infection. IMPORTANCE Tuberculosis (TB), caused by Mycobacterium tuberculosis, is one of the leading causes of death worldwide. The rise in drug resistance highlights the urgent need for innovation in anti-TB drug development. Many anti-TB drugs require bioactivation by Baeyer-Villiger monooxygenases (BVMOs). Despite their emerging importance, BVMO structural and functional features remain enigmatic. We applied a comprehensive bioinformatic analysis and confirmed the presence of six BVMOs in M. tuberculosis, including MymA, EthA, and Rv0565c-activators of the second-line prodrug ethionamide-and the novel BVMO Rv0892. Combining in silico characterization with in vitro validation, we outlined their structural framework and substrate preference. Markedly, MymA displayed an enhanced capacity and a distinct selectivity profile toward ligands, in agreement with its catalytic site topology. These features ground the molecular basis for structure-function comprehension of the specificity in these enzymes and expand the repertoire of BVMOs with selective and/or overlapping activity for application in the context of improving anti-TB therapy.

Developmental Modulation of Root Cell Wall Architecture Confers Resistance to an Oomycete Pathogen.

The cell wall is the primary interface between plant cells and their immediate environment and must balance multiple functionalities, including the regulation of growth, the entry of beneficial microbes, and protection against pathogens. Here, we demonstrate how API, a SCAR2 protein component of the SCAR/WAVE complex, controls the root cell wall architecture important for pathogenic oomycete and symbiotic bacterial interactions in legumes. A mutation in API results in root resistance to the pathogen Phytophthora palmivora and colonization defects by symbiotic rhizobia. Although api mutant plants do not exhibit significant overall growth and development defects, their root cells display delayed actin and endomembrane trafficking dynamics and selectively secrete less of the cell wall polysaccharide xyloglucan. Changes associated with a loss of API establish a cell wall architecture with altered biochemical properties that hinder P. palmivora infection progress. Thus, developmental stage-dependent modifications of the cell wall, driven by SCAR/WAVE, are important in balancing cell wall developmental functions and microbial invasion.

HpaP Sequesters HrpJ, an Essential Component of Ralstonia solanacearum Virulence That Triggers Necrosis in Arabidopsis.

The Gram-negative bacterium Ralstonia solanacearum, the causal agent of bacterial wilt, is a worldwide major crop pathogen whose virulence strongly relies on a type III secretion system (T3SS). This extracellular apparatus allows the translocation of proteins, called type III effectors (T3Es), directly into the host cells. To date, very few data are available in plant-pathogenic bacteria concerning the role played by type III secretion (T3S) regulators at the posttranslational level. We have demonstrated that HpaP, a putative T3S substrate specificity switch protein of R. solanacearum, controls T3E secretion. To better understand the role of HpaP on T3S control, we analyzed the secretomes of the GMI1000 wild-type strain as well as the hpaP mutant using a mass spectrometry experiment (liquid chromatography tandem mass spectrometry). The secretomes of both strains appeared to be very similar and highlighted the modulation of the secretion of few type III substrates. Interestingly, only one type III-associated protein, HrpJ, was identified as specifically secreted by the hpaP mutant. HrpJ appeared to be an essential component of the T3SS, essential for T3S and pathogenicity. We further showed that HrpJ is specifically translocated in planta by the hpaP mutant and that HrpJ can physically interact with HpaP. Moreover, confocal microscopy experiments demonstrated a cytoplasmic localization for HrpJ once in planta. When injected into Arabidopsis thaliana leaves, HrpJ is able to trigger a necrosis on 16 natural accessions. A genome-wide association mapping revealed a major association peak with 12 highly significant single-nucleotide polymorphisms located on a plant acyl-transferase.

A secreted WY-domain-containing protein present in European isolates of the oomycete Plasmopara viticola induces cell death in grapevine and tobacco species.

Plasmopara viticola is a biotrophic oomycete pathogen causing grapevine downy mildew. We characterized the repertoire of P. viticola effector proteins which may be translocated into plants to support the disease. We found several secreted proteins that contain canonical dEER motifs and conserved WY-domains but lack the characteristic RXLR motif reported previously from oomycete effectors. We cloned four candidates and showed that one of them, Pv33, induces plant cell death in grapevine and Nicotiana species. This activity is dependent on the nuclear localization of the protein. Sequence similar effectors were present in seven European, but in none of the tested American isolates. Together our work contributes a new type of conserved P. viticola effector candidates.

Spontaneous mutations in a regulatory gene induce phenotypic heterogeneity and adaptation of Ralstonia solanacearum to changing environments.

An evolution experiment with the bacterial plant pathogen Ralstonia solanacearum revealed that several adaptive mutations conferring enhanced fitness in plants arose in the efpR gene encoding a regulator of virulence and metabolic functions. In this study, we found that an efpR mutant systematically displays colonies with two morphotypes: the type S ('smooth', similar to the wild type) and the type EV ('efpR variant'). We demonstrated that the efpH gene, a homologue of efpR, plays a key role in the control of phenotypic heterogeneity, the ΔefpR-ΔefpH double mutant being stably locked into the EV type. Using mixed infection assays, we demonstrated that the type EV is metabolically more proficient than the type S and displays fitness gain in specific environments, whereas the type S has a better fitness into the plant environment. We provide evidence that this efpR-dependent phenotypic heterogeneity is a general feature of strains of the R. solanacearum species complex and could occur in natural conditions. This study highlights the potential role of phenotypic heterogeneity in this plant pathogen as an adaptive trait to changing environments.

Endosymbiotic Sinorhizobium meliloti modulate Medicago root susceptibility to secondary infection via ethylene.

A complex network of pathways coordinates nodulation and epidermal root hair infection in the symbiotic interaction between rhizobia and legume plants. Whereas nodule formation was known to be autoregulated, it was so far unclear whether a similar control is exerted on the infection process. We assessed the capacity of Medicago plants nodulated by Sinorhizobium meliloti to modulate root susceptibility to secondary bacterial infection or to purified Nod factors in split-root and volatile assays using bacterial and plant mutant combinations. Ethylene implication in this process emerged from gas production measurements, use of a chemical inhibitor of ethylene biosynthesis and of a Medicago mutant affected in ethylene signal transduction. We identified a feedback mechanism that we named AOI (for Autoregulation Of Infection) by which endosymbiotic bacteria control secondary infection thread formation by their rhizospheric peers. AOI involves activation of a cyclic adenosine 3',5'-monophosphate (cAMP) cascade in endosymbiotic bacteria, which decreases both root infectiveness and root susceptibility to bacterial Nod factors. These latter two effects are mediated by ethylene. AOI is a novel component of the complex regulatory network controlling the interaction between Sinorhizobium meliloti and its host plants that emphasizes the implication of endosymbiotic bacteria in fine-tuning the interaction.

Sunflower resistance to multiple downy mildew pathotypes revealed by recognition of conserved effectors of the oomycete Plasmopara halstedii.

Over the last 40 years, new sunflower downy mildew isolates (Plasmopara halstedii) have overcome major gene resistances in sunflower, requiring the identification of additional and possibly more durable broad-spectrum resistances. Here, 354 RXLR effectors defined in silico from our new genomic data were classified in a network of 40 connected components sharing conserved protein domains. Among 205 RXLR effector genes encoding conserved proteins in 17 P. halstedii pathotypes of varying virulence, we selected 30 effectors that were expressed during plant infection as potentially essential genes to target broad-spectrum resistance in sunflower. The transient expression of the 30 core effectors in sunflower and in Nicotiana benthamiana leaves revealed a wide diversity of targeted subcellular compartments, including organelles not so far shown to be targeted by oomycete effectors such as chloroplasts and processing bodies. More than half of the 30 core effectors were able to suppress pattern-triggered immunity in N. benthamiana, and five of these induced hypersensitive responses (HR) in sunflower broad-spectrum resistant lines. HR triggered by PhRXLRC01 co-segregated with Pl22 resistance in F3 populations and both traits localized in 1.7 Mb on chromosome 13 of the sunflower genome. Pl22 resistance was physically mapped on the sunflower genome recently sequenced, unlike all the other downy mildew resistances published so far. PhRXLRC01 and Pl22 are proposed as an avirulence/resistance gene couple not previously described in sunflower. Core effector recognition is a successful strategy to accelerate broad-spectrum resistance gene identification in complex crop genomes such as sunflower.

Comparative transcriptomic studies identify specific expression patterns of virulence factors under the control of the master regulator PhcA in the Ralstonia solanacearum species complex.

The global regulator PhcA controls numerous traits associated to virulence and bacterial proliferation in strains of the plant pathogen Ralstonia solanacearum species complex. Here, we conducted a genome-wide RNA sequencing study of the GMI1000 wild-type strain and a derived phcA mutant grown in complete medium. The PhcA regulon we identified is the largest regulon described to date in the R. solanacearum species complex with 1581 regulated genes, representing about 30% of the bacterial genome. Among these genes, 166 transcription regulators were identified including known regulators controlling major cellular functions such as the Type 3 secretion system and 27 novel regulators that were not identified in previous transcriptomic studies. This study highlights that PhcA controls other functions beside pathogenicity stricto sensu which participate to the global cell homeostasis (metabolism, energy storage). We then compared the PhcA regulon identified in complete medium to the recently published PhcA regulon obtained in planta. This comparison of the set of GMI1000 genes subjected to PhcA regulation in both conditions revealed 383 common genes. Among them, 326 (85%) had a similar PhcA dependent regulation pattern in complete medium and in planta, and 57 (15%) displayed an opposite regulation pattern. A large majority of the genes repressed by PhcA in complete medium but activated in planta belong to the HrpG-HrpB regulon, which represents a set of key genes required for R. solanacearum pathogenesis. This latter class of genes appears to be specifically induced by PhcA in the plant environment whereas PhcA represses their expression in complete medium. The large set of direct and indirect targets identified in this study will contribute to enrich our knowledge of the intricate regulatory network coordinating the expression of virulence and metabolic functions in the model plant pathogen R. solanacearum.

Nod factors potentiate auxin signaling for transcriptional regulation and lateral root formation in Medicago truncatula.

Nodulation (Nod) factors (NFs) are symbiotic molecules produced by rhizobia that are essential for establishment of the rhizobium-legume endosymbiosis. Purified NFs can stimulate lateral root formation (LRF) in Medicago truncatula, but little is known about the molecular mechanisms involved. Using a combination of reporter constructs, pharmacological and genetic approaches, we show that NFs act on early steps of LRF in M. truncatula, independently of the ethylene signaling pathway and of the cytokinin receptor MtCRE1, but in interaction with auxin. We conducted a whole-genome transcriptomic study upon NF and/or auxin treatments, using a lateral root inducible system adapted for M. truncatula. This revealed a large overlap between NF and auxin signaling and, more interestingly, synergistic interactions between these molecules. Three groups showing interaction effects were defined: group 1 contained more than 1500 genes responding specifically to the combinatorial treatment of NFs and auxin; group 2 comprised auxin-regulated genes whose expression was enhanced or antagonized by NFs; and in group 3 the expression of NF regulated genes was antagonized by auxin. Groups 1 and 2 were enriched in signaling and metabolic functions, which highlights important crosstalk between NF and auxin signaling for both developmental and symbiotic processes.

Enhanced in planta Fitness through Adaptive Mutations in EfpR, a Dual Regulator of Virulence and Metabolic Functions in the Plant Pathogen Ralstonia solanacearum.

Experimental evolution of the plant pathogen Ralstonia solanacearum, where bacteria were maintained on plant lineages for more than 300 generations, revealed that several independent single mutations in the efpR gene from populations propagated on beans were associated with fitness gain on bean. In the present work, novel allelic efpR variants were isolated from populations propagated on other plant species, thus suggesting that mutations in efpR were not solely associated to a fitness gain on bean, but also on additional hosts. A transcriptomic profiling and phenotypic characterization of the efpR deleted mutant showed that EfpR acts as a global catabolic repressor, directly or indirectly down-regulating the expression of multiple metabolic pathways. EfpR also controls virulence traits such as exopolysaccharide production, swimming and twitching motilities and deletion of efpR leads to reduced virulence on tomato plants after soil drenching inoculation. We studied the impact of the single mutations that occurred in efpR during experimental evolution and found that these allelic mutants displayed phenotypic characteristics similar to the deletion mutant, although not behaving as complete loss-of-function mutants. These adaptive mutations therefore strongly affected the function of efpR, leading to an expanded metabolic versatility that should benefit to the evolved clones. Altogether, these results indicated that EfpR is a novel central player of the R. solanacearum virulence regulatory network. Independent mutations therefore appeared during experimental evolution in the evolved clones, on a crucial node of this network, to favor adaptation to host vascular tissues through regulatory and metabolic rewiring.

Single-cell analysis reveals new subset markers of murine peritoneal macrophages and highlights macrophage dynamics upon Staphylococcus aureus peritonitis.

Resident macrophages play a central role in maintaining tissue homeostasis and immune surveillance. Here, we used single cell-based qPCR coupled with flow cytometry analysis to further define the phenotypes of large and small resident peritoneal macrophages (LPMs and SPMs, respectively) in mice. We demonstrated that the expression of Cxcl13, IfngR1, Fizz-1 and Mrc-1 clearly distinguished between LPMs and SPMs subsets. Using these markers, the dynamics of peritoneal macrophages in a Staphylococcus aureus-induced peritonitis model were analyzed. We found that S. aureus infection triggers a massive macrophage disappearance reaction in both subsets. Thereafter, inflammatory monocytes rapidly infiltrated the cavity and differentiated to replenish the SPMs. Although phenotypically indistinguishable from resident SPMs by flow cytometry, newly recruited SPMs had a different pattern of gene expression dominated by M2 markers combined with M1 associated features (inos expression). Interestingly, S. aureus elicited SPMs showed a robust expression of Cxcl13, suggesting that these cells may endorse the role of depleted LPMs and contribute to restoring peritoneal homeostasis. These data provide information on both resident and recruited macrophages dynamics upon S. aureus infection and demonstrate that single-cell phenotyping is a promising and highly valuable approach to unraveling macrophage diversity and plasticity.

Comparative Secretome Analysis of Ralstonia solanacearum Type 3 Secretion-Associated Mutants Reveals a Fine Control of Effector Delivery, Essential for Bacterial Pathogenicity.

Ralstonia solanacearum, the causal agent of bacterial wilt, exerts its pathogenicity through more than a hundred secreted proteins, many of them depending directly on the functionality of a type 3 secretion system. To date, only few type 3 effectors have been identified as required for bacterial pathogenicity, notably because of redundancy among the large R. solanacearum effector repertoire. In order to identify groups of effectors collectively promoting disease on susceptible hosts, we investigated the role of putative post-translational regulators in the control of type 3 secretion. A shotgun secretome analysis with label-free quantification using tandem mass spectrometry was performed on the R. solanacearum GMI1000 strain. There were 228 proteins identified, among which a large proportion of type 3 effectors, called Rip (Ralstonia injected proteins). Thanks to this proteomic approach, RipBJ was identified as a new effector specifically secreted through type 3 secretion system and translocated into plant cells. A focused Rip secretome analysis using hpa (hypersensitive response and pathogenicity associated) mutants revealed a fine secretion regulation and specific subsets of Rips with different secretion patterns. We showed that a set of Rips (RipF1, RipW, RipX, RipAB, and RipAM) are secreted in an Hpa-independent manner. We hypothesize that these Rips could be preferentially involved in the first stages of type 3 secretion. In addition, the secretion of about thirty other Rips is controlled by HpaB and HpaG. HpaB, a candidate chaperone was shown to positively control secretion of numerous Rips, whereas HpaG was shown to act as a negative regulator of secretion. To evaluate the impact of altered type 3 effectors secretion on plant pathogenesis, the hpa mutants were assayed on several host plants. HpaB was required for bacterial pathogenicity on multiple hosts whereas HpaG was found to be specifically required for full R. solanacearum pathogenicity on the legume plant Medicago truncatula.

Combined genetic and transcriptomic analysis reveals three major signalling pathways activated by Myc-LCOs in Medicago truncatula.

Myc-LCOs are newly identified symbiotic signals produced by arbuscular mycorrhizal (AM) fungi. Like rhizobial Nod factors, they are lipo-chitooligosaccharides that activate the common symbiotic signalling pathway (CSSP) in plants. To increase our limited understanding of the roles of Myc-LCOs we aimed to analyse Myc-LCO-induced transcriptional changes and their genetic control. Whole genome RNA sequencing (RNA-seq) was performed on roots of Medicago truncatula wild-type plants, and dmi3 and nsp1 symbiotic mutants affected in nodulation and mycorrhizal signalling. Plants were treated separately with the two major types of Myc-LCOs, sulphated and nonsulphated. Generalized linear model analysis identified 2201 differentially expressed genes and classified them according to genotype and/or treatment effects. Three genetic pathways for Myc-LCO-regulation of transcriptomic reprogramming were highlighted: DMI3- and NSP1-dependent; DMI3-dependent and NSP1-independent; and DMI3- and NSP1-independent. Comprehensive analysis revealed overlaps with previous AM studies, and highlighted certain functions, especially signalling components and transcription factors. These data provide new insights into mycorrhizal signalling mechanisms, supporting a role for NSP1, and specialisation for NSP1-dependent and -independent pathways downstream of DMI3. Our data also indicate significant Myc-LCO-activated signalling upstream of DMI3 and/or parallel to the CSSP and some constitutive activity of the CSSP.

Bridging physiological and evolutionary time-scales in a gene regulatory network.

Gene regulatory networks (GRNs) govern phenotypic adaptations and reflect the trade-offs between physiological responses and evolutionary adaptation that act at different time-scales. To identify patterns of molecular function and genetic diversity in GRNs, we studied the drought response of the common sunflower, Helianthus annuus, and how the underlying GRN is related to its evolution. We examined the responses of 32,423 expressed sequences to drought and to abscisic acid (ABA) and selected 145 co-expressed transcripts. We characterized their regulatory relationships in nine kinetic studies based on different hormones. From this, we inferred a GRN by meta-analyses of a Gaussian graphical model and a random forest algorithm and studied the genetic differentiation among populations (FST ) at nodes. We identified two main hubs in the network that transport nitrate in guard cells. This suggests that nitrate transport is a critical aspect of the sunflower physiological response to drought. We observed that differentiation of the network genes in elite sunflower cultivars is correlated with their position and connectivity. This systems biology approach combined molecular data at different time-scales and identified important physiological processes. At the evolutionary level, we propose that network topology could influence responses to human selection and possibly adaptation to dry environments.

Combined linkage and association mapping of flowering time in Sunflower (Helianthus annuus L.).

Association mapping and linkage mapping were used to identify quantitative trait loci (QTL) and/or causative mutations involved in the control of flowering time in cultivated sunflower Helianthus annuus. A panel of 384 inbred lines was phenotyped through testcrosses with two tester inbred lines across 15 location × year combinations. A recombinant inbred line (RIL) population comprising 273 lines was phenotyped both per se and through testcrosses with one or two testers in 16 location × year combinations. In the association mapping approach, kinship estimation using 5,923 single nucleotide polymorphisms was found to be the best covariate to correct for effects of panel structure. Linkage disequilibrium decay ranged from 0.08 to 0.26 cM for a threshold of 0.20, after correcting for structure effects, depending on the linkage group (LG) and the ancestry of inbred lines. A possible hitchhiking effect is hypothesized for LG10 and LG08. A total of 11 regions across 10 LGs were found to be associated with flowering time, and QTLs were mapped on 11 LGs in the RIL population. Whereas eight regions were demonstrated to be common between the two approaches, the linkage disequilibrium approach did not detect a documented QTL that was confirmed using the linkage mapping approach.

A gene-phenotype network based on genetic variability for drought responses reveals key physiological processes in controlled and natural environments.

Identifying the connections between molecular and physiological processes underlying the diversity of drought stress responses in plants is key for basic and applied science. Drought stress response involves a large number of molecular pathways and subsequent physiological processes. Therefore, it constitutes an archetypical systems biology model. We first inferred a gene-phenotype network exploiting differences in drought responses of eight sunflower (Helianthus annuus) genotypes to two drought stress scenarios. Large transcriptomic data were obtained with the sunflower Affymetrix microarray, comprising 32423 probesets, and were associated to nine morpho-physiological traits (integrated transpired water, leaf transpiration rate, osmotic potential, relative water content, leaf mass per area, carbon isotope discrimination, plant height, number of leaves and collar diameter) using sPLS regression. Overall, we could associate the expression patterns of 1263 probesets to six phenotypic traits and identify if correlations were due to treatment, genotype and/or their interaction. We also identified genes whose expression is affected at moderate and/or intense drought stress together with genes whose expression variation could explain phenotypic and drought tolerance variability among our genetic material. We then used the network model to study phenotypic changes in less tractable agronomical conditions, i.e. sunflower hybrids subjected to different watering regimes in field trials. Mapping this new dataset in the gene-phenotype network allowed us to identify genes whose expression was robustly affected by water deprivation in both controlled and field conditions. The enrichment in genes correlated to relative water content and osmotic potential provides evidence of the importance of these traits in agronomical conditions.