Transcriptomic responses of Serratia liquefaciens cells grown under simulated Martian conditions of low temperature, low pressure, and CO2-enriched anoxic atmosphere

Results from previous experiments indicated that the Gram-negative α-proteobacterium Serratia liquefaciens strain ATCC 27592 was capable of growth under low temperature (0 °C), low pressure (0.7 kPa), and anoxic, CO₂-dominated atmosphere–conditions intended to simulate the near-subsurface environment of Mars. To probe the response of its transcriptome to this extreme environment, S. liquefaciens ATCC 27592 was cultivated under 4 different environmental simulations: 0 °C, 0.7 kPa, CO₂ atmosphere (Condition A); 0 °C, ~101.3 kPa, CO₂ atmosphere (Condition B); 0 °C, ~101.3 kPa, ambient N₂/O₂ atmosphere (Condition C); and 30 °C, ~101.3 kPa, N₂/O₂ atmosphere (Condition D; ambient laboratory conditions). RNA-seq was performed on ribosomal RNA-depleted total RNA isolated from triplicate cultures grown under Conditions A-D and the datasets generated were subjected to transcriptome analyses. The data from Conditions A, B, or C were compared to laboratory Condition D. Significantly differentially expressed transcripts were identified belonging to a number of KEGG pathway categories. Up-regulated genes under all Conditions A, B, and C included those encoding transporters (ABC and PTS transporters); genes involved in translation (ribosomes and their biogenesis, biosynthesis of both tRNAs and aminoacyl-tRNAs); DNA repair and recombination; and non-coding RNAs. Genes down-regulated under all Conditions A, B, and C included: transporters (mostly ABC transporters); flagellar and motility proteins; genes involved in phenylalanine metabolism; transcription factors; and two-component systems. The results are discussed in the context of Mars astrobiology and planetary protection.

Transcriptomic responses to extreme low salinity among locally adapted populations of Olympia oyster (Ostrea lurida)

The Olympia oyster (Ostrea lurida) is a foundation species inhabiting estuaries along the North American west coast. In California estuaries, O. lurida is adapted to local salinity regimes and populations differ in low salinity tolerance. In this study, oysters from three California populations were reared for two generations in a laboratory common garden and subsequently exposed to low salinity seawater. Comparative transcriptomics was then used to understand species-level responses to hyposmotic stress and population-level mechanisms underlying divergent salinity tolerances. Gene expression patterns indicate Olympia oysters are sensitive to hyposmotic stress: All populations respond to low salinity by up-regulating transcripts indicative of protein unfolding, DNA damage and cell cycle arrest after sub-lethal exposure. Among O. lurida populations, transcriptomic profiles differed constitutively and in response to low salinity. Despite two generations in common-garden conditions, transcripts encoding apoptosis modulators were constitutively expressed at significantly different levels in the most tolerant population. Expression of cell death regulators may facilitate cell fate decisions when salinity declines. Following low salinity exposure, oysters from the more tolerant population expressed a small number of mRNAs at significantly higher levels than less tolerant populations. Proteins encoded by these transcripts regulate ciliary activity within the mantle cavity and may function to prolong valve closure and reduce mortality in low salinity seawater. Collectively, gene expression patterns suggest sub-lethal impacts of hyposmotic stress in Olympia oysters are considerable and that even oysters with greater low salinity tolerance may be vulnerable to future freshwater flooding events.

Coexpression network revealing the plasticity and robustness of population transcriptome during the initial stage of domesticating energy crop Miscanthus lutarioriparius

Coexpression network revealing genes with Co-variation Expression pattern (CE) and those with Top rank of Expression fold change (TE) played different roles in responding to new environment of Miscanthus lutarioriparius. Variation in gene expression level, the product of genetic and/or environmental perturbation, determines the robustness-to-plasticity spectrum of a phenotype in plants. Understanding how expression variation of plant population response to a new field is crucial to domesticate energy crops. Weighted Gene Coexpression Network Analysis (WGCNA) was used to explore the patterns of expression variation based on 72 Miscanthus lutarioriparius transcriptomes from two contrasting environments, one near the native habitat and the other in one harsh domesticating region. The 932 genes with Co-variation Expression pattern (CE) and other 932 genes with Top rank of Expression fold change (TE) were identified and the former were strongly associated with the water use efficiency (r ≥ 0.55, P ≤ 10^-7). Functional enrichment of CE genes were related to three organelles, which well matched the annotation of twelve motifs identified from their conserved noncoding sequence; while TE genes were mostly related to biotic and/or abiotic stress. The expression robustness of CE genes with high genetic diversity kept relatively stable between environments while the harsh environment reduced the expression robustness of TE genes with low genetic diversity. The expression plasticity of CE genes was increased less than that of TE genes. These results suggested that expression variation of CE genes and TE genes could account for the robustness and plasticity of acclimation ability of Miscanthus, respectively. The patterns of expression variation revealed by transcriptomic network would shed new light on breeding and domestication of energy crops.

Evolution of heat-shock protein expression underlying adaptive responses to environmental stress

Heat-shock proteins (Hsps) and their cognates are primary mitigators of cell stress. With increasingly severe impacts of climate change and other human modifications of the biosphere, the ability of the heat-shock system to affect evolutionary fitness in environments outside the laboratory and to evolve in response is topic of growing importance. Since the last major reviews, several advances have occurred. First, demonstrations of the heat-shock response outside the laboratory now include many additional taxa and environments. Many of these demonstrations are only correlative, however. More importantly, technical advances in "omic" quantification of nucleic acids and proteins, genomewide association analysis, and manipulation of genes and their expression have enabled the field to move beyond correlation. Several consequent advances are already evident: The pathway from heat-shock gene expression to stress tolerance in nature can be extremely complex, mediated through multiple biological processes and systems, and even multiple species. The underlying genes are more numerous, diverse and variable than previously appreciated, especially with respect to their regulatory variation and epigenetic changes. The impacts and limitations (e.g., due to trade-offs) of natural selection on these genes have become more obvious and better established. At last, as evolutionary capacitors, Hsps may have distinctive impacts on the evolution of other genes and ecological consequences.

Plasticity of female reproductive resource allocation depends on the presence or absence of prior environmental sex determination in Ceratopteris richardii

Resource allocation plasticity enables individuals to alter patterns of nutrient use between reproductive and vegetative output to better fit their current environment. In sexually labile plant species, abiotic environmental factors can influence expression of dimorphic gender, resulting in environmental sex determination (ESD), which potentially reduces the need for plasticity of resource allocation by preemptively matching an individual's future nutrient demands to resource availability in its location. Ceratopteris richardii gametophytes exhibit gender-dependent differences in relative carbon and nitrogen content, and ESD in certain nutrient environments. This study examined whether prior ESD in C. richardii gametophyte populations reduced subsequent plasticity of reproductive allocation compared to instances where no ESD occurred, by quantifying phenotypic responses to reduced P, N, or CO 2 availabilities. All three nutrient-limited environments resulted in decreased size of egg-bearing (meristic) gametophytes compared to nonlimited environments, but gametophytes failed to respond to N and CO 2 limitation at the time of sex determination, resulting in no ESD. N limitation resulted in a predictable allometric re-allocation of resources based on small gametophyte size, whereas CO 2 limitation caused a change in reproductive output consistent with true plasticity. Withholding exogenous P caused ESD and had no effect on relative reproductive output of resultant meristic gametophytes because the size decrease was minor. Under P limitation, ESD matched the resource demands of gender phenotypes to their environment before the onset of developmental dimorphism, reducing the need for large allocation adjustments after sex determination.

Phenotypic and transcriptomic responses to salinity stress across genetically and geographically divergent Tigriopus californicus populations

Species inhabiting the North American west coast intertidal must tolerate an extremely variable environment, with large fluctuations in both temperature and salinity. Uncovering the mechanisms for this tolerance is key to understanding species' persistence. We tested for differences in salinity tolerance between populations of Tigriopus californicus copepods from locations in northern (Bodega Reserve) and southern (San Diego) California known to differ in temperature, precipitation and humidity. We also tested for differences between populations in their transcriptomic responses to salinity. Although these two populations have ~20% mtDNA sequence divergence and differ strongly in other phenotypic traits, we observed similarities in their phenotypic and transcriptomic responses to low and high salinity stress. Salinity significantly affected respiration rate (increased under low salinity and reduced under high salinity), but we found no significant effect of population on respiration or a population by salinity interaction. Under high salinity, there was no population difference in knock-down response, but northern copepods had a smaller knock-down under low salinity stress, corroborating previous results for T. californicus. Northern and southern populations had a similar transcriptomic response to salinity based on a principle components analysis, although differential gene expression under high salinity stress was three times lower in the northern population compared to the southern population. Transcripts differentially regulated under salinity stress were enriched for "amino acid transport" and "ion transport" annotation categories, supporting previous work demonstrating that the accumulation of free amino acids is important for osmotic regulation in T. californicus.

Comparative transcriptomics of cyprinid minnows and carp in a common wild setting: a resource for ecological genomics in freshwater communities

Comparative transcriptomics can now be conducted on organisms in natural settings, which has greatly enhanced understanding of genome–environment interactions. Here, we demonstrate the utility and potential pitfalls of comparative transcriptomics of wild organisms, with an example from three cyprinid fish species (Teleostei:Cypriniformes). We present extensively filtered and annotated transcriptome assemblies that provide a valuable resource for studies of genome evolution (e.g. polyploidy), ecological and morphological diversification, speciation, and shared and unique responses to environmental variation in cyprinid fishes. Our results and analyses address the following points: (i) ‘essential developmental genes’ are shown to be ubiquitously expressed in a diverse suite of tissues across later ontogenetic stages (i.e. juveniles and adults), making these genes are useful for assessing the quality of transcriptome assemblies, (ii) the influence of microbiomes and other exogenous DNA, (iii) potentially novel, species-specific genes, and (iv) genomic rearrangements (e.g. whole genome duplication). The data we present provide a resource for future comparative work in cypriniform fishes and other taxa across a variety of sub-disciplines, including stress response, morphological diversification, community ecology, ecotoxicology, and climate change.

Modularity Facilitates Flexible Tuning of Plastic and Evolutionary Gene Expression Responses during Early Divergence

Gene expression changes have been recognized as important drivers of adaptation to changing environmental conditions. Little is known about the relative roles of plastic and evolutionary responses in complex gene expression networks during the early stages of divergence. Large gene expression data sets coupled with in silico methods for identifying coexpressed modules now enable systems genetics approaches also in nonmodel species for better understanding of gene expression responses during early divergence. Here, we combined gene coexpression analyses with population genetics to separate plastic and population (evolutionary) effects in expression networks using small salmonid populations as a model system. We show that plastic and population effects were highly variable among the six identified modules and that the plastic effects explained larger proportion of the total eigengene expression than population effects. A more detailed analysis of the population effects using a QST - FST comparison across 16,622 annotated transcripts revealed that gene expression followed neutral expectations within modules and at the global level. Furthermore, two modules showed enrichment for genes coding for early developmental traits that have been previously identified as important phenotypic traits in thermal responses in the same model system indicating that coexpression analysis can capture expression patterns underlying ecologically important traits. We suggest that module-specific responses may facilitate the flexible tuning of expression levels to local thermal conditions. Overall, our study indicates that plasticity and neutral evolution are the main drivers of gene expression variance in the early stages of thermal adaptation in this system.

Functional genomic analysis of corals from natural CO2 -seeps reveals core molecular responses involved in acclimatization to ocean acidification

Little is known about the potential for acclimatization or adaptation of corals to ocean acidification and even less about the molecular mechanisms underpinning these processes. Here, we examine global gene expression patterns in corals and their intracellular algal symbionts from two replicate population pairs in Papua New Guinea that have undergone long-term acclimatization to natural variation in pCO₂ . In the coral host, only 61 genes were differentially expressed in response to pCO₂ environment, but the pattern of change was highly consistent between replicate populations, likely reflecting the core expression homeostasis response to ocean acidification. Functional annotations highlight lipid metabolism and a change in the stress response capacity of corals as key parts of this process. Specifically, constitutive downregulation of molecular chaperones was observed, which may impact response to combined climate change-related stressors. Elevated CO₂ has been hypothesized to benefit photosynthetic organisms but expression changes of in hospite Symbiodinium in response to acidification were greater and less consistent among reef populations. This population-specific response suggests hosts may need to adapt not only to an acidified environment, but also to changes in their Symbiodinium populations that may not be consistent among environments, adding another challenging dimension to the physiological process of coping with climate change.

Transcriptomic response of the Antarctic pteropod Limacina helicina antarctica to ocean acidification

Background

Ocean acidification (OA), a change in ocean chemistry due to the absorption of atmospheric CO₂ into surface oceans, challenges biogenic calcification in many marine organisms. Ocean acidification is expected to rapidly progress in polar seas, with regions of the Southern Ocean expected to experience severe OA within decades. Biologically, the consequences of OA challenge calcification processes and impose an energetic cost.

Results

In order to better characterize the response of a polar calcifier to conditions of OA, we assessed differential gene expression in the Antarctic pteropod, Limacina helicina antarctica. Experimental levels of pCO₂ were chosen to create both contemporary pH conditions, and to mimic future pH expected in OA scenarios. Significant changes in the transcriptome were observed when juvenile L. h. antarctica were acclimated for 21 days to low-pH (7.71), mid-pH (7.9) or high-pH (8.13) conditions. Differential gene expression analysis of individuals maintained in the low-pH treatment identified down-regulation of genes involved in cytoskeletal structure, lipid transport, and metabolism. High pH exposure led to increased expression and enrichment for genes involved in shell formation, calcium ion binding, and DNA binding. Significant differential gene expression was observed in four major cellular and physiological processes: shell formation, the cellular stress response, metabolism, and neural function. Across these functional groups, exposure to conditions that mimic ocean acidification led to rapid suppression of gene expression.

Conclusions

Results of this study demonstrated that the transcriptome of the juvenile pteropod, L. h. antarctica, was dynamic and changed in response to different levels of pCO₂. In a global change context, exposure of L. h. antarctica to the low pH, high pCO₂ OA conditions resulted in a suppression of transcripts for genes involved in key physiological processes: calcification, metabolism, and the cellular stress response. The transcriptomic response at both acute and longer-term acclimation time frames indicated that contemporary L. h. antarctica may not have the physiological plasticity necessary for adaptation to OA conditions expected in future decades. Lastly, the differential gene expression results further support the role of shelled pteropods such as L. h. antarctica as sentinel organisms for the impacts of ocean acidification.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4161-0) contains supplementary material, which is available to authorized users.

Conserved Transcriptional Responses to Nutrient Stress in Bloom-Forming Algae

The concentration and composition of bioavailable nitrogen (N) and phosphorus (P) in the upper ocean shape eukaryotic phytoplankton communities and influence their physiological responses. Phytoplankton are known to exhibit similar physiological responses to limiting N and P conditions such as decreased growth rates, chlorosis, and increased assimilation of N and P. Are these responses similar at the molecular level across multiple species? To interrogate this question, five species from biogeochemically important, bloom-forming taxa (Bacillariophyta, Dinophyta, and Haptophyta) were grown under similar low N, low P, and replete nutrient conditions to identify transcriptional patterns and associated changes in biochemical pools related to N and P stress. Metabolic profiles, revealed through the transcriptomes of these taxa, clustered together based on species rather than nutrient stressor, suggesting that the global metabolic response to nutrient stresses was largely, but not exclusively, species-specific. Nutrient stress led to few transcriptional changes in the two dinoflagellates, consistent with other research. An orthologous group analysis examined functionally conserved (i.e., similarly changed) responses to nutrient stress and therefore focused on the diatom and haptophytes. Most conserved ortholog changes were specific to a single nutrient treatment, but a small number of orthologs were similarly changed under both N and P stress in 2 or more species. Many of these orthologs were related to photosynthesis and may represent generalized stress responses. A greater number of orthologs were conserved across more than one species under low P compared to low N. Screening the conserved orthologs for functions related to N and P metabolism revealed increased relative abundance of orthologs for nitrate, nitrite, ammonium, and amino acid transporters under N stress, and increased relative abundance of orthologs related to acquisition of inorganic and organic P substrates under P stress. Although the global transcriptional responses were dominated by species-specific changes, the analysis of conserved responses revealed functional similarities in resource acquisition pathways among different phytoplankton taxa. This overlap in nutrient stress responses observed among species may be useful for tracking the physiological ecology of phytoplankton field populations.

Assessing elements of an extended evolutionary synthesis for plant domestication and agricultural origin research

The development of agricultural societies, one of the most transformative events in human and ecological history, was made possible by plant and animal domestication. Plant domestication began 12,000-10,000 y ago in a number of major world areas, including the New World tropics, Southwest Asia, and China, during a period of profound global environmental perturbations as the Pleistocene epoch ended and transitioned into the Holocene. Domestication is at its heart an evolutionary process, and for many prehistorians evolutionary theory has been foundational in investigating agricultural origins. Similarly, geneticists working largely with modern crops and their living wild progenitors have documented some of the mechanisms that underwrote phenotypic transformations from wild to domesticated species. Ever-improving analytic methods for retrieval of empirical data from archaeological sites, together with advances in genetic, genomic, epigenetic, and experimental research on living crop plants and wild progenitors, suggest that three fields of study currently little applied to plant domestication processes may be necessary to understand these transformations across a range of species important in early prehistoric agriculture. These fields are phenotypic (developmental) plasticity, niche construction theory, and epigenetics with transgenerational epigenetic inheritance. All are central in a controversy about whether an Extended Evolutionary Synthesis is needed to reconceptualize how evolutionary change occurs. An exploration of their present and potential utility in domestication study shows that all three fields have considerable promise in elucidating important issues in plant domestication and in agricultural origin and dispersal research and should be increasingly applied to these issues.

Evolutionary Transcriptomics and Proteomics: Insight into Plant Adaptation

Comparative transcriptomics and proteomics (T&P) have brought biological insight into development, gene function, and physiological stress responses. However, RNA-seq and high-throughput proteomics remain underutilised in studies of plant adaptation. These methodologies have created discovery tools with the potential to significantly advance our understanding of adaptive diversification. We outline experimental recommendations for their application. We discuss analysis models and approaches that accelerate the identification of adaptive gene sets and integrate transcriptome, proteome, phenotypic, and environmental data. Finally, we encourage widespread uptake and future developments in T&P that will advance our understanding of evolution and adaptation.

Locally adapted populations of a copepod can evolve different gene expression patterns under the same environmental pressures

As populations diverge in allopatry, but under similar thermal conditions, do similar thermal performance phenotypes evolve by maintaining similar gene expression patterns, or does genetic divergence lead to divergent patterns of gene expression between these populations? We used genetically divergent populations of the copepod Tigriopus californicus, whose performance at different thermal conditions is well characterized, to investigate transcriptome-wide expression responses under two different thermal regimes: (1) a nonvariable temperature regime and (2) a regime with variable temperature. Our results show the expression profiles of the response to these regimes differed substantially among populations, even for populations that are geographically close. This pattern was accentuated when populations were raised in the variable temperature environment. Less heat-tolerant populations mounted strong but divergent responses to the different thermal regimes, with a large heat-shock response observed in one population, and an apparent reduction in the expression of genes involved in basic cellular processes in the other. Our results suggest that as populations diverge in allopatry, they may evolve starkly different responses to changes in temperature, at the gene expression level, while maintaining similar thermal performance phenotypes.

Transcriptomic resources for an endemic Neotropical plant lineage (Gesneriaceae)

PREMISE OF THE STUDY

Despite the extensive phenotypic variation that characterizes the Gesneriaceae family, there is a lack of genomic resources to investigate the molecular basis of their diversity. We developed and compared the transcriptomes for two species of the Neotropical lineage of the Gesneriaceae.

METHODS AND RESULTS

Illumina sequencing and de novo assembly of floral and leaf samples were used to generate multigene sequence data for Sinningia eumorpha and S. magnifica, two species endemic to the Brazilian Atlantic Forest. A total of 300 million reads were used to assemble the transcriptomes, with an average of 92,038 transcripts and 43,506 genes per species. The transcriptomes showed good quality metrics, with the presence of all eukaryotic core genes, and an equal representation of clusters of orthologous groups (COG) classifications between species. The orthologous search produced 8602 groups, with 15-20% of them annotated using BLAST tools.

DISCUSSION

This study provides the first step toward a comprehensive multispecies transcriptome characterization of the Gesneriaceae family. These resources are the basis for comparative analyses in this species-rich Neotropical plant group; they will also allow the investigation of the evolutionary importance of multiple metabolic pathways and phenotypic diversity, as well as developmental programs in these nonmodel species.

Aquatic insects in a multistress environment: cross-tolerance to salinity and desiccation

Exposing organims to a particular stressor may enhance tolerance to a subsequent stress, when protective mechanisms against both stressors are shared. Such cross-tolerance is a common adaptive response in dynamic multivariate environments and often indicates potential co-evolution of stress traits. Many aquatic insects in inland saline waters from Mediterranean-climate regions are sequentially challenged with salinity and desiccation stress. Thus, cross-tolerance to these physiologically similar stressors could have been positively selected in insects of these regions. We used adults of the saline water beetles Enochrus jesusarribasi (Hydrophilidae) and Nebrioporus baeticus (Dytiscidae) to test cross-tolerance responses to desiccation and salinity. In independent laboratory experiments, we evaluated the effects of i) salinity stress on the subsequent resistance to desiccation and ii) desiccation stress (rapid and slow dehydration) on the subsequent tolerance to salinity. Survival, water loss and haemolymph osmolality were measured. Exposure to stressful salinity improved water control under subsequent desiccation stress in both species, with a clear cross-tolerance (enhanced performance) in N. baeticus. In contrast, general negative effects on performance were found under the inverse stress sequence. The rapid and slow dehydration produced different water loss and haemolymph osmolality dynamics that were reflected in different survival patterns. Our finding of cross-tolerance to salinity and desiccation in ecologically similar species from distant lineages, together with parallel responses between salinity and thermal stress previously found in several aquatic taxa, highlights the central role of adaption to salinity and co-occurring stressors in arid inland waters, having important implications for the species' persistence under climate change.

Adaptation to climate change: trade-offs among responses to multiple stressors in an intertidal crustacean

Trade-offs may influence both physiological and evolutionary responses to co-occurring stressors, but their effects on both plastic and adaptive responses to climate change are poorly understood. To test for genetic and physiological trade-offs incurred in tolerating multiple stressors, we hybridized two populations of the intertidal copepod Tigriopus californicus that were divergent for both heat and salinity tolerance. Starting in the F₂ generation, we selected for increased tolerance of heat, low salinity, and high salinity in replicate lines. After five generations of selection, heat-selected lines had greater heat tolerance but lower fecundity, indicating an energetic cost to tolerance. Lines selected for increased salinity tolerance did not show evidence of adaptation to their respective environments; however, hypo-osmotic selection lines showed substantial loss of tolerance to hyperosmotic stress. Neither of the salinity selection regimes resulted in diminished heat tolerance at ambient salinity; however, simultaneous exposure to heat and hypo-osmotic stress led to decreased heat tolerance, implying a physiological trade-off in tolerance to the two stressors. When we quantified the transcriptomic response to heat and salinity stress via RNA sequencing, we observed little overlap in the stress responses, suggesting the observed synergistic effects of heat and salinity stress were driven by competing energetic demands, rather than shared stress response pathways.

Killing two birds with one stone: Model plant systems as a tool to teach the fundamental concepts of gene expression while analyzing biological data

Plants are ideal systems to teach core biology concepts due to their unique physiological and developmental features. Advances in DNA sequencing technology and genomics have allowed scientists to generate genome sequences and transcriptomics data for numerous model plant species. This information is publicly available and presents a valuable tool to introduce undergraduate students to the fundamental concepts of gene expression in the context of modern quantitative biology and bioinformatics. Modern biology classrooms must provide authentic research experiences to allow developing core competencies such as scientific inquiry, critical interpretation of experimental results, and quantitative analyses of large dataset using computational approaches. Recent educational research has shown that undergraduate students struggle when connecting gene expression concepts to classic genetics, phenotypic analyses, and overall flow of biological information in living organisms, suggesting that novel approaches are necessary to enhance learning of gene expression and regulation. This review describes different strategies and resources available to instructors willing to incorporate authentic research experiences, genomic tools, and bioinformatics analyses when teaching transcriptional regulation and gene expression in undergraduate courses. A variety of laboratory exercises and pedagogy materials developed to teach gene expression using plants are discussed. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.