An official website of the United States government
Abstract The regulation of gene expression is central to many biological processes. Gene regulatory networks (GRNs) link transcription factors (TFs) to their target genes and represent maps of potential transcriptional regulation. Here, we analyzed a large number of publically available maize (Zea mays) transcriptome data sets including >6000 RNA sequencing samples to generate 45 coexpression-based GRNs that represent potential regulatory relationships between TFs and other genes in different populations of samples (cross-tissue, cross-genotype, and tissue-and-genotype samples). While these networks are all enriched for biologically relevant interactions, different networks capture distinct TF-target associations and biological processes. By examining the power of our coexpression-based GRNs to accurately predict covarying TF-target relationships in natural variation data sets, we found that presence/absence changes rather than quantitative changes in TF gene expression are more likely associated with changes in target gene expression. Integrating information from our TF-target predictions and previous expression quantitative trait loci (eQTL) mapping results provided support for 68 TFs underlying 74 previously identified trans-eQTL hotspots spanning a variety of metabolic pathways. This study highlights the utility of developing multiple GRNs within a species to detect putative regulators of important plant pathways and provides potential targets for breeding or biotechnological applications.
more »
« less
Polygenic pathogen networks influence transcriptional plasticity in the Arabidopsis–Botrytis pathosystem
Abstract Bidirectional flow of information shapes the outcome of the host–pathogen interactions and depends on the genetics of each organism. Recent work has begun to use co-transcriptomic studies to shed light on this bidirectional flow, but it is unclear how plastic the co-transcriptome is in response to genetic variation in both the host and pathogen. To study co-transcriptome plasticity, we conducted transcriptomics using natural genetic variation in the pathogen, Botrytis cinerea, and large-effect genetic variation abolishing defense signaling pathways within the host, Arabidopsis thaliana. We show that genetic variation in the pathogen has a greater influence on the co-transcriptome than mutations that abolish defense signaling pathways in the host. Genome-wide association mapping using the pathogens’ genetic variation and both organisms’ transcriptomes allowed an assessment of how the pathogen modulates plasticity in response to the host. This showed that the differences in both organism's responses were linked to trans-expression quantitative trait loci (eQTL) hotspots within the pathogen's genome. These hotspots control gene sets in either the host or pathogen and show differential allele sensitivity to the host’s genetic variation rather than qualitative host specificity. Interestingly, nearly all the trans-eQTL hotspots were unique to the host or pathogen transcriptomes. In this system of differential plasticity, the pathogen mediates the shift in the co-transcriptome more than the host.
more »
« less
- Award ID(s):
- 2020754
- PAR ID:
- 10432805
- Editor(s):
- Birchler, J
- Date Published:
- Journal Name:
- GENETICS
- Volume:
- 224
- Issue:
- 3
- ISSN:
- 1943-2631
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Trans-Acting Genotypes Associated with mRNA Expression Affect Metabolic and Thermal Tolerance TraitsBetancourt, Andrea (Ed.)Abstract Evolutionary processes driving physiological trait variation depend on the underlying genomic mechanisms. Evolution of these mechanisms depends on the genetic complexity (involving many genes) and how gene expression impacting the traits is converted to phenotype. Yet, genomic mechanisms that impact physiological traits are diverse and context dependent (e.g., vary by environment and tissues), making them difficult to discern. We examine the relationships between genotype, mRNA expression, and physiological traits to discern the genetic complexity and whether the gene expression affecting the physiological traits is primarily cis- or trans-acting. We use low-coverage whole genome sequencing and heart- or brain-specific mRNA expression to identify polymorphisms directly associated with physiological traits and expressed quantitative trait loci (eQTL) indirectly associated with variation in six temperature specific physiological traits (standard metabolic rate, thermal tolerance, and four substrate specific cardiac metabolic rates). Focusing on a select set of mRNAs belonging to co-expression modules that explain up to 82% of temperature specific traits, we identified hundreds of significant eQTL for mRNA whose expression affects physiological traits. Surprisingly, most eQTL (97.4% for heart and 96.7% for brain) were trans-acting. This could be due to higher effect size of trans- versus cis-acting eQTL for mRNAs that are central to co-expression modules. That is, we may have enhanced the identification of trans-acting factors by looking for single nucleotide polymorphisms associated with mRNAs in co-expression modules that broadly influence gene expression patterns. Overall, these data indicate that the genomic mechanism driving physiological variation across environments is driven by trans-acting heart- or brain-specific mRNA expression.more » « less
-
IntroductionSeagrass meadows serve as an integral component of coastal ecosystems but are declining rapidly due to numerous anthropogenic stressors including climate change. Eelgrass wasting disease, caused by opportunisticLabyrinthulaspp., is an increasing concern with rising seawater temperature. To better understand the host-pathogen interaction, we paired whole organism physiological assays with dual transcriptomic analysis of the infected host and parasite. MethodsEelgrass (Zostera marina) shoots were placed in one of two temperature treatments, 11° C or 18° C, acclimated for 10 days, and exposed to a waterborne inoculation containing infectiousLabyrinthula zosterae(Lz) or sterile seawater. At two- and five-days post-exposure, pathogen load, visible disease signs, whole leaf phenolic content, and both host- and pathogen- transcriptomes were characterized. ResultsTwo days after exposure, more than 90% of plants had visible lesions andLzDNA was detectable in 100% percent of sampled plants in theLzexposed treatment. Concentrations of total phenolic compounds were lower after 5 days of combined exposure to warmer temperatures andLz, but were unaffected in other treatments. Concentrations of condensed tannins were not affected byLzor temperature, and did not change over time. Analysis of the eelgrass transcriptome revealed 540 differentially expressed genes in response toLzexposure, but not temperature.Lz-exposed plants had gene expression patterns consistent with increased defense responses through altered regulation of phytohormone biosynthesis, stress response, and immune function pathways. Analysis of the pathogen transcriptome revealed up-regulation of genes potentially involved in breakdown of host defense, chemotaxis, phagocytosis, and metabolism. DiscussionThe lack of a significant temperature signal was unexpected but suggests a more pronounced physiological response toLzinfection as compared to temperature. Pre-acclimation of eelgrass plants to the temperature treatments may have contributed to the limited physiological responses to temperature. Collectively, these data characterize a widespread physiological response to pathogen attack and demonstrate the value of paired transcriptomics to understand infections in a host-pathogen system.more » « less
-
Abstract The evolutionary fate of multi-strain pathogens is shaped by host-pathogen ecological interactions. In bacterial pathogens of plants, enhanced strain characterization and advances in our understanding of molecular mechanisms underlying defense pathways open the door for revisiting the role of negative frequency-dependent selection (NFDS) in strain structure, including its interplay with genetic exchange. NFDS arising from specific defense is one potential mechanism for generating, maintaining, and structuring pathogen diversity. In plants, specific protection against microbial pathogens involves Resistance proteins (R-proteins) that recognize virulence factors (effectors) secreted by pathogens, typically to subvert the initial line of host defense. Here we formulate a stochastic computational co-evolution model that explicitly incorporates variable length R-gene and effector repertoires, and migration from their regional pools. We use this model to understand potential mechanisms shaping effector repertoire structure and associated strain coexistence in the generalist plant pathogenP. syringae. The demonstration of a modular structure in our numerical simulations motivates the analysis of genome sequences from 76 strains collected in the Midwestern US and 1104 strains from global sources. We find that effector repertories both locally and globally exhibit a modular structure, with higher similarity within than between clusters. The observed modules are consistent with the core genome phylogeny and are unexplained by plant host species, location of isolation, and genetic linkage between effectors. An extension of the model is needed to take into account the evidence for genetic exchange and the phylogenetic congruence of effector modules. We initialize the system with a phylogenetically congruent modular structure and include recombination rates decreasing as a function of phylogenetic distance. We show that NFDS can counter-balance the effects of mixing due to recombination and in so doing, contributes to the maintenance of strain structure. These findings indicate that the observed similarity clusters may constitute, in part, emergent niches arising from eco-evolutionary dynamics that contribute to strain coexistence.more » « less
-
Plants are continuously exposed to beneficial and pathogenic microbes, but how plants recognize and respond to friends ver- sus foes remains poorly understood. Here, we compared the molecular response of Arabidopsis thaliana independently chal- lenged with a Fusarium oxysporum endophyte Fo47 versus a pathogen Fo5176. These two F. oxysporum strains share a core genome of about 46 Mb, in addition to 1,229 and 5,415 unique accessory genes. Metatranscriptomic data reveal a shared pattern of expression for most plant genes (about 80%) in responding to both fungal inoculums at all timepoints from 12 to 96 h postinoculation (HPI). However, the distinct responding genes depict transcriptional plasticity, as the path- ogenic interaction activates plant stress responses and sup- presses functions related to plant growth and development, while the endophytic interaction attenuates host immunity but activates plant nitrogen assimilation. The differences in reprogramming of the plant transcriptome are most obvious in 12 HPI, the earliest timepoint sampled, and are linked to accessory genes in both fungal genomes. Collectively, our results indicate that the A. thaliana and F. oxysporum interac- tion displays both transcriptome conservation and plasticity in the early stages of infection, providing insights into the fine- tuning of gene regulation underlying plant differential responses to fungal endophytes and pathogens.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1093/genetics/iyad099
