Multivariable regulation of gene expression plasticity in metazoans

Gene Expression Shifts in Emperor Penguin Adaptation to the Extreme Antarctic Environment

Gene expression can accelerate ecological divergence by rapidly tweaking the response of an organism to novel environments, with more divergent environments exerting stronger selection and supposedly, requiring faster adaptive responses. Organisms adapted to extreme environments provide ideal systems to test this hypothesis, particularly when compared to related species with milder ecological niches. The Emperor penguin (Aptenodytes forsteri) is the only endothermic vertebrate breeding in the harsh Antarctic winter, in stark contrast with the less cold-adapted sister species, the King penguin (A. patagonicus). Assembling the first de novo transcriptomes and analysing multi-tissue (brain, kidney, liver, muscle, skin) RNA-Seq data from natural populations of both species, we quantified the shifts in tissue-enhanced genes, co-expression gene networks, and differentially expressed genes characterising Emperor penguin adaptation to the extreme Antarctic. Our analyses revealed the crucial role played by muscle and liver in temperature homeostasis, fasting, and whole-body energy metabolism (glucose/insulin regulation, lipid metabolism, fatty acid beta-oxidation, and blood coagulation). Repatterning at the regulatory level appears as more important in the brain of the Emperor penguin, showing the lowest signature of differential gene expression, but the largest co-expression gene network shift. Nevertheless, over-expressed genes related to mTOR signalling in the brain and the liver support their central role in cold and fasting responses. Besides contributing to understanding the genetics underlying complex traits, like body energy reservoir management, our results provide a first insight into the role of gene expression in adaptation to one of the most extreme environmental conditions endured by an endotherm.

Plasticity for the win: Flexible transcriptional response to host plant switches in the comma butterfly (Polygonia c-album)

Generalist plant-feeding insects are characterised by a broad host repertoire that can comprise several families or even different orders of plants. The genetic and physiological mechanisms underlying the use of such a wide host range are still not fully understood. Earlier studies indicate that the consumption of different host plants is associated with host-specific gene expression profiles. It remained, however, unclear if and how larvae can alter these profiles in the case of a changing host environment. Using the polyphagous comma butterfly (Polygonia c-album) we show that larvae can adjust their transcriptional profiles in response to a new host plant. The switch to some of the host plants, however, resulted in a larger transcriptional response and, thus, seems to be more challenging. At a physiological level, no correspondence for these patterns could be found in larval performance. This suggests that a high transcriptional but also phenotypic flexibility are essential for the use of a broad and diverse host range. We furthermore propose that host switch tests in the laboratory followed by transcriptomic investigations can be a valuable tool to examine not only plasticity in host use but also subtle and/or transient trade-offs in the evolution of host plant repertoires.

Patterns of methylation and transcriptional plasticity during thermal acclimation in a reef-building coral

Phenotypic plasticity can buffer organisms against short-term environmental fluctuations. For example, previous exposure to increased temperatures can increase thermal tolerance in many species. Prior studies have found that acclimation to higher temperature can influence the magnitude of transcriptional response to subsequent acute thermal stress (hereafter, "transcriptional response modulation"). However, mechanisms mediating this gene expression response and, ultimately, phenotypic plasticity remain largely unknown. Epigenetic modifications are good candidates for modulating transcriptional response, as they broadly correlate with gene expression. Here, we investigate changes in DNA methylation as a possible mechanism controlling shifts in gene expression plasticity and thermal acclimation in the reef-building coral Acropora nana. We find that gene expression response to acute stress is altered in corals acclimated to different temperatures, with many genes exhibiting a dampened response to heat stress in corals pre-conditioned to higher temperatures. At the same time, we observe shifts in methylation during both acclimation (11 days) and acute heat stress (24 h). We observed that the acute heat stress results in shifts in gene-level methylation and elicits an acute transcriptional response in distinct gene sets. Further, acclimation-induced shifts in gene expression plasticity and differential methylation also largely occur in separate sets of genes. Counter to our initial hypothesis no overall correlation between the magnitude of differential methylation and the change in gene expression plasticity. We do find a small but statistically significant overlap in genes exhibiting both dampened expression response and shifts in methylation (14 genes), which could be candidates for further inquiry. Overall, our results suggest transcriptional response modulation occurs independently from methylation changes induced by thermal acclimation.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) - the variability of gene expression in response to environmental changes - remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also find that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Genome-wide mapping of regulatory variants for temperature- and salinity-adaptive genes reveals genetic basis of genotype-by-environment interaction in Crassostrea ariakensis

Regulatory variants in gene expression serve as bridges linking genetic variation and phenotypic plasticity. Environmental conditions typically influence the effects of regulatory variants on phenotypic plasticity; however, such genotype-by-environment interactions (G × E) are poorly understood. This study aimed to investigate the genetic basis of G × E in estuarine oyster (Crassostrea ariakensis), which is an important model animal for studying environmental adaption owing to its high plasticity and large intraspecific divergence. Genome-wide mapping of expression quantitative trait loci (eQTLs) for 23 environmental adaptive genes was performed for 256 estuarine oysters. We identified 1194 eQTL single nucleotide polymorphisms (eSNPs), including 433 cis-eSNPs in four genes and 722 trans-eSNPs in eight genes. The expression variation explanation of cis-eSNPs (9.95%) was significantly higher than that of trans-eSNPs (9.15%). We specifically showed cis- and trans-eSNPs with high linkage disequilibrium (LD) for Traf7, Slc6a5, Ggt, and Dap3. For example, we identified a cis-regulatory LD block containing 68 cis-eSNP and a trans-regulatory LD block, including 20 trans-eSNPs in Traf7. A high proportion (85%) of 40 vital eSNPs exhibited significant G × E effects. We identified crossing and nonparallel interactions of G × E, with the tag cis-eSNPs of Baat and Slc6a5 as representatives. Our results indicated that cis-eQTLs are highly conserved. This study provides insights into the understanding of adaptive evolutionary mechanisms and phenotypic response prediction to variable environments, as well as the genetic improvement for superior adaptive traits for genetic resource conservation and aquaculture.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Genetically encoded transcriptional plasticity underlies stress adaptation in Mycobacterium tuberculosis

Transcriptional regulation is a critical adaptive mechanism that allows bacteria to respond to changing environments, yet the concept of transcriptional plasticity (TP) remains largely unexplored. In this study, we investigate the genome-wide TP profiles of Mycobacterium tuberculosis (Mtb) genes by analyzing 894 RNA sequencing samples derived from 73 different environmental conditions. Our data reveal that Mtb genes exhibit significant TP variation that correlates with gene function and gene essentiality. We also found that critical genetic features, such as gene length, GC content, and operon size independently impose constraints on TP, beyond trans-regulation. By extending our analysis to include two other Mycobacterium species -- M. smegmatis and M. abscessus -- we demonstrate a striking conservation of the TP landscape. This study provides a comprehensive understanding of the TP exhibited by mycobacteria genes, shedding light on this significant, yet understudied, genetic feature encoded in bacterial genomes.

Weak gene-gene interaction facilitates the evolution of gene expression plasticity

BACKGROUND

Individual organisms may exhibit phenotypic plasticity when they acclimate to different conditions. Such plastic responses may facilitate or constrain the adaptation of their descendant populations to new environments, complicating their evolutionary trajectories beyond the genetic blueprint. Intriguingly, phenotypic plasticity itself can evolve in terms of its direction and magnitude during adaptation. However, we know little about what determines the evolution of phenotypic plasticity, including gene expression plasticity. Recent laboratory-based studies suggest dominance of reversing gene expression plasticity-plastic responses that move the levels of gene expression away from the new optima. Nevertheless, evidence from natural populations is still limited.

RESULTS

Here, we studied gene expression plasticity and its evolution in the montane and lowland populations of an elevationally widespread songbird-the Rufous-capped Babbler (Cyanoderma ruficeps)-with reciprocal transplant experiments and transcriptomic analyses; we set common gardens at altitudes close to these populations' native ranges. We confirmed the prevalence of reversing plasticity in genes associated with altitudinal adaptation. Interestingly, we found a positive relationship between magnitude and degree of evolution in gene expression plasticity, which was pertinent to not only adaptation-associated genes but also the whole transcriptomes from multiple tissues. Furthermore, we revealed that genes with weaker expressional interactions with other genes tended to exhibit stronger plasticity and higher degree of plasticity evolution, which explains the positive magnitude-evolution relationship.

CONCLUSIONS

Our experimental evidence demonstrates that species may initiate their adaptation to new habitats with genes exhibiting strong expression plasticity. We also highlight the role of expression interdependence among genes in regulating the magnitude and evolution of expression plasticity. This study illuminates how the evolution of phenotypic plasticity in gene expression facilitates the adaptation of species to challenging environments in nature.

DOT-1.1 (DOT1L) deficiency in C. elegans leads to small RNA-dependent gene activation

Methylation of histone H3 at lysine 79 (H3K79) is conserved from yeast to humans and is accomplished by Dot1 (disruptor of telomeric silencing-1) methyltransferases. The C. elegans enzyme DOT-1.1 and its interacting partners are similar to the mammalian DOT1L (Dot1-like) complex. The C. elegans DOT-1.1 complex has been functionally connected to RNA interference. Specifically, we have previously shown that embryonic and larval lethality of dot-1.1 mutant worms deficient in H3K79 methylation was suppressed by mutations in the RNAi pathway genes responsible for generation (rde-4) and function (rde-1) of primary small interfering RNAs (siRNAs). This suggests that dot-1.1 mutant lethality is dependent on the enhanced production of some siRNAs. We have also found that this lethality is suppressed by a loss-of-function of CED-3, a conserved apoptotic protease. Here, we describe a comparison of gene expression and primary siRNA production changes between control and dot-1.1 deletion mutant embryos. We found that elevated antisense siRNA production occurred more often at upregulated than downregulated genes. Importantly, gene expression changes were dependent on RDE-4 in both instances. Moreover, the upregulated group, which is potentially activated by ectopic siRNAs, was enriched in protease-coding genes. Our findings are consistent with a model where in the absence of H3K79 methylation there is a small RNA-dependent activation of protease genes, which leads to embryonic and larval lethality. DOT1 enzymes' conservation suggests that the interplay between H3K79 methylation and small RNA pathways may exist in higher organisms.

Cis- and Trans-variations of Stearoyl-CoA Desaturase Provide New Insights into the Mechanisms of Diverged Pattern of Phenotypic Plasticity for Temperature Adaptation in Two Congeneric Oyster Species

The evolution of phenotypic plasticity plays an essential role in adaptive responses to climate change; however, its regulatory mechanisms in marine organisms which exhibit high phenotypic plasticity still remain poorly understood. The temperature-responsive trait oleic acid content and its major gene stearoyl-CoA desaturase (Scd) expression have diverged in two allopatric congeneric oyster species, cold-adapted Crassostrea gigas and warm-adapted Crassostrea angulata. In this study, genetic and molecular methods were used to characterize fatty acid desaturation and membrane fluidity regulated by oyster Scd. Sixteen causative single-nucleotide polymorphisms (SNPs) were identified in the promoter/cis-region of the Scd between wild C. gigas and C. angulata. Further functional experiments showed that an SNP (g.-333C [C. gigas allele] >T [C. angulata allele]) may influence Scd transcription by creating/disrupting the binding motif of the positive trans-factor Y-box factor in C. gigas/C. angulata, which mediates the higher/lower constitutive expression of Scd in C. gigas/C. angulata. Additionally, the positive trans-factor sterol-regulatory element-binding proteins (Srebp) were identified to specifically bind to the promoter of Scd in both species, and were downregulated during cold stress in C. gigas compared to upregulated in C. angulata. This partly explains the relatively lower environmental sensitivity (plasticity) of Scd in C. gigas. This study serves as an experimental case to reveal that both cis- and trans-variations shape the diverged pattern of phenotypic plasticity, which provides new insights into the formation of adaptive traits and the prediction of the adaptive potential of marine organisms to future climate change.

Data-driven identification of inherent features of eukaryotic stress-responsive genes

Living organisms are continuously challenged by changes in their environment that can propagate to stresses at the cellular level, such as rapid changes in osmolarity or oxygen tension. To survive these sudden changes, cells have developed stress-responsive mechanisms that tune cellular processes. The response of Saccharomyces cerevisiae to osmostress includes a massive reprogramming of gene expression. Identifying the inherent features of stress-responsive genes is of significant interest for understanding the basic principles underlying the rewiring of gene expression upon stress. Here, we generated a comprehensive catalog of osmostress-responsive genes from 5 independent RNA-seq experiments. We explored 30 features of yeast genes and found that 25 (83%) were distinct in osmostress-responsive genes. We then identified 13 non-redundant minimal osmostress gene traits and used statistical modeling to rank the most stress-predictive features. Intriguingly, the most relevant features of osmostress-responsive genes are the number of transcription factors targeting them and gene conservation. Using data on HeLa samples, we showed that the same features that define yeast osmostress-responsive genes can predict osmostress-responsive genes in humans, but with changes in the rank-ordering of feature-importance. Our study provides a holistic understanding of the basic principles of the regulation of stress-responsive gene expression across eukaryotes.

Flavonoids Modulate the Accumulation of Toxins From Aspergillus flavus in Maize Kernels

Aspergillus flavus is an opportunistic fungal pathogen capable of producing aflatoxins, potent carcinogenic toxins that accumulate in maize kernels after infection. To better understand the molecular mechanisms of maize resistance to A. flavus growth and aflatoxin accumulation, we performed a high-throughput transcriptomic study in situ using maize kernels infected with A. flavus strain 3357. Three maize lines were evaluated: aflatoxin-contamination resistant line TZAR102, semi-resistant MI82, and susceptible line Va35. A modified genotype-environment association method (GEA) used to detect loci under selection via redundancy analysis (RDA) was used with the transcriptomic data to detect genes significantly influenced by maize line, fungal treatment, and duration of infection. Gene ontology enrichment analysis of genes highly expressed in infected kernels identified molecular pathways associated with defense responses to fungi and other microbes such as production of pathogenesis-related (PR) proteins and lipid bilayer formation. To further identify novel genes of interest, we incorporated genomic and phenotypic field data from a genome wide association analysis with gene expression data, allowing us to detect significantly expressed quantitative trait loci (eQTL). These results identified significant association between flavonoid biosynthetic pathway genes and infection by A. flavus. In planta fungal infections showed that the resistant line, TZAR102, has a higher fold increase of the metabolites naringenin and luteolin than the susceptible line, Va35, when comparing untreated and fungal infected plants. These results suggest flavonoids contribute to plant resistance mechanisms against aflatoxin contamination through modulation of toxin accumulation in maize kernels.

Inter-embryo gene expression variability recapitulates the hourglass pattern of evo-devo

BACKGROUND

The evolution of embryological development has long been characterized by deep conservation. In animal development, the phylotypic stage in mid-embryogenesis is more conserved than either early or late stages among species within the same phylum. Hypotheses to explain this hourglass pattern have focused on purifying the selection of gene regulation. Here, we propose an alternative-genes are regulated in different ways at different stages and have different intrinsic capacities to respond to perturbations on gene expression.

RESULTS

To eliminate the influence of natural selection, we quantified the expression variability of isogenetic single embryo transcriptomes throughout fly Drosophila melanogaster embryogenesis. We found that the expression variability is lower at the phylotypic stage, supporting that the underlying regulatory architecture in this stage is more robust to stochastic variation on gene expression. We present evidence that the phylotypic stage is also robust to genetic variations on gene expression. Moreover, chromatin regulation appears to play a key role in the variation and evolution of gene expression.

CONCLUSIONS

We suggest that a phylum-level pattern of embryonic conservation can be explained by the intrinsic difference of gene regulatory mechanisms in different stages.

Caenorhabditis elegans Deficient in DOT-1.1 Exhibit Increases in H3K9me2 at Enhancer and Certain RNAi-Regulated Regions

The methylation of histone H3 at lysine 79 is a feature of open chromatin. It is deposited by the conserved histone methyltransferase DOT1. Recently, DOT1 localization and H3K79 methylation (H3K79me) have been correlated with enhancers in C. elegans and mammalian cells. Since earlier research implicated H3K79me in preventing heterochromatin formation both in yeast and leukemic cells, we sought to inquire whether a H3K79me deficiency would lead to higher levels of heterochromatic histone modifications, specifically H3K9me2, at developmental enhancers in C. elegans. Therefore, we used H3K9me2 ChIP-seq to compare its abundance in control and dot-1.1 loss-of-function mutant worms, as well as in rde-4; dot-1.1 and rde-1; dot-1.1 double mutants. The rde-1 and rde-4 genes are components of the RNAi pathway in C. elegans, and RNAi is known to initiate H3K9 methylation in many organisms, including C. elegans. We have previously shown that dot-1.1(-) lethality is rescued by rde-1 and rde-4 loss-of-function. Here we found that H3K9me2 was elevated in enhancer, but not promoter, regions bound by the DOT-1.1/ZFP-1 complex in dot-1.1(-) worms. We also found increased H3K9me2 at genes targeted by the ALG-3/4-dependent small RNAs and repeat regions. Our results suggest that ectopic H3K9me2 in dot-1.1(-) could, in some cases, be induced by small RNAs.