PubMed HealthSearch

SEARCH · PubMed Health

Results for “Environmental gradients”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

Environmental Gradients as a Dominant Force in the Macroevolution of a Host-Associated Marine Bacterium.

Natural selection is imposed by both abiotic environmental filtering and biotic interactions, yet their relative roles in shaping the deep phylogeny of widespread, generalist host-associated bacteria remain unclear. Here, we integrate large-scale phylogenomics, environmental sequencing, functional genomics, and global metagenomic analysis to demonstrate that tidal zonation overrides host association as the dominant macroevolutionary force structuring the marine bacterial genus Ruegeria. Analysis of 533 genomes and 74 global coastal metagenomes reveals that the intertidal-subtidal boundary structures the deepest phylogenetic splits, driving the repeated evolution of distinct ecotypes through independent zonation transitions across global coastlines. These ecotypes possess divergent genomic toolkits: intertidal strains are enriched for genes coding for stress resistance and anaerobic metabolism, whereas subtidal strains specialize in high-affinity nutrient scavenging. Our findings establish that predictable physicochemical gradients act as filters that generate foundational diversity from which specialized host symbionts subsequently emerge, reframing how environmental gradients shape microbial evolution at the eco-evolutionary interface.

Journal Article

Global Patterns of Net Ecosystem Exchange in peatlands: A Systematic Review and Meta-analysis of Drivers Across Land Use and Environmental Gradients.

Peatlands play an essential role in the global carbon cycle, storing approximately one-third of the world's soil carbon despite covering less than 3% of the land surface. Peatland degradation from anthropogenic activities and climate change can convert peatlands from net carbon sinks to sources by altering carbon cycling. Net Ecosystem Exchange (NEE), the balance between CO2 uptake and emission, is a critical indicator for assessing peatland condition and restoration efforts. We conducted a systematic quantitative literature review to investigate global patterns of NEE in peatlands and identify key environmental and anthropogenic drivers of CO2 flux variability. Annual NEE values from 120 globally distributed sites reported in peer-reviewed literature were analyzed in relation to climatic zone, land use, vegetation type, peatland condition, and water table depth. Our synthesis revealed significant geographic gaps, with peatland NEE studies substantially underrepresented in the Tropics, Africa, and Oceania. Agricultural peatlands emitted significantly more CO2 than sites under natural land uses or peat extraction, while degraded peatlands were significantly greater net CO2 sources than intact and restored systems. Restored peatlands remained net CO2 sources on average, emphasizing the importance of long-term monitoring and adaptive management following restoration interventions. Water table depth significantly affected NEE variability, with CO2 emissions increasing approximately 7.2 gCO2-C m-2yr-1 for every centimeter of water table drawdown. A substantial variability in measurement methods, data processing software, and protocols highlighted the critical need for methodological standardization. Our findings provide evidence-based targets for peatland conservation and restoration monitoring as nature-based climate solutions.

Ecosystem

Hidden genomic structure and widespread structural polymorphism across environmental gradients in the spiny sea star Marthasterias glacialis.

Genomic regions of reduced recombination can preserve linkage among co-adapted alleles, facilitating local adaptation despite high connectivity. Such regions-often generated by chromosomal inversions-may be especially important in highly dispersive marine taxa yet remain poorly documented in echinoderms. Here, we combined a chromosome-level reference genome with genome-wide ddRAD-seq from 296 Marthasterias glacialis individuals across 19 Atlantic-Mediterranean locations to quantify population structure and scan for recombination-suppressed haploblocks. Genome-wide neutral markers showed significant population differentiation together with evidence of high connectivity, revealed by the presence of inter-ecoregion migrants. Additionally, we identified 16 polymorphic haploblocks with patterns consistent with putative chromosomal inversions spanning 18.6% of the genome. Haploblock haplotypes were strongly environmentally and geographically structured and contained genes with key functions in stress response, osmoregulation and thermal tolerance. Haplotype distributions also paralleled previously described mitochondrial lineages despite nuclear gene flow, consistent with a model of ancient divergence followed by secondary contact. Overall, our results suggest a role for widespread structural polymorphism in adaptive differentiation in Echinodermata, providing a framework for linking echinoderm genome rearrangements to ecological divergence. Marthasterias glacialis thus emerges as a promising system to explore how structural variation contributes to adaptation and genome evolution in highly dispersive organisms.

Animals

Genomic Insights Into Local Adaptation Across Heterogeneous Understory Habitats and Climate Change Vulnerability.

Understanding adaptive evolution and survival risks in understory herbs is crucial for the effective conservation of biodiversity. How environmental gradients shape species local adaptation patterns is not well understood, nor is how populations of understory herbs respond to a changing climate. In this study, we conducted population genomic analyses of Adenocaulon himalaicum (Asteraceae) with a pan-East Asian distribution, representing a good model for dominant understory herbs to elucidate adaptation mechanisms in heterogeneous forest ecosystems. Based on 34,398 putatively neutral single nucleotide polymorphisms (SNPs) across 27 populations, we identified three genetic lineages accompanied by high levels of genetic differentiation between populations. Our isolation by environment results (IBE) indicated a significant effect of environmental gradients on genomic variation of A. himalaicum (r = 0.18, p = 0.03). To decompose the relative contributions of climate, geography and population structure in explaining genetic variance, our partial RDA found that the prominent contribution of environmental effects (climatic and soil variables) explained 29% and 36% of the neutral and adaptive genetic variation, respectively. Using two genotype-environment association (GEA) methods, we identified 13 SNPs as candidates for core climate-related adaptation loci, with two of these loci further validated by qRT-PCR experiments. Projections of spatiotemporal genomic vulnerability under different future climate scenarios revealed that populations in the southeastern edge of the Himalayas, near the Sichuan Basin, the southernmost region of Northeast China and the northern Korean Peninsula, as well as northern Japan, were identified as the most vulnerable and should be prioritised for conservation. Therefore, our current study provides the genomic foundations for conservation and management strategies to elucidate how these understory herbs cope with future climate changes.

Climate Change

Climate and soil shape Daqu wheat quality and seed microbiome via rhizosphere taxa and microbial assembly.

The grain quality and seed microbiome of Daqu wheat are fundamental determinants of Daqu fermentation performance; however, the mechanisms by which cultivation environments influence these traits via rhizosphere microbial communities remain unclear. Bacterial and fungal communities across the bulk soil-rhizosphere-seed continuum of three wheat cultivars grown in four ecoregions were characterized using absolute quantitative amplicon sequencing. The rhizosphere microbiome was treated as a central intermediary, while the response variables were seed microbial diversity and grain-quality traits, including starch content, protein content, and grain hardness. Twelve physicochemical properties of soil and 11 climatic factors were integrated into a multidimensional association framework. Environmental conditions exerted stronger influences on both seed quality traits and microbial diversity than cultivar identity. Distinct regional signatures were also evident in rhizosphere microbiomes, with environmental gradients explaining community variation more effectively than geographic distance. Bacterial communities exhibited greater sensitivity to environmental fluctuations than fungi. Mantel analyses identified available nitrogen, precipitation, and atmospheric pressure as significant drivers of core rhizosphere taxa (P&#xa0;<&#xa0;0.05). iCAMP revealed that stochastic processes predominantly governed rhizosphere bacterial assembly, whereas stochastic and deterministic mechanisms jointly shaped fungal assembly. Partial least squares path modeling further uncovered a rhizosphere-mediated environment-seed cascade, wherein sunlight intensity and duration, atmospheric pressure, and soil nitrogen directly or indirectly affected seed wet gluten content, grain hardness, and seed microbial diversity through their influences on rhizosphere microbiota. Rhizosphere bacterial diversity was negatively associated with seed bacterial diversity (path coefficient&#xa0;=&#xa0;-0.118, P&#xa0;<&#xa0;0.05), indicating that rhizosphere communities may shape seed endophytic bacterial assemblages via environmental filtering and competitive interactions. Collectively, these findings elucidate how environments shape the quality and seed microbiomes of Daqu wheat, providing scientific guidance for optimal site selection and the standardized production of high-quality brewing wheat for industrial Baijiu.

Triticum

A trait-based ecological perspective on the soil microbial antibiotic-related genetic machinery.

Antibiotic resistance crisis dictates the need for resistance monitoring and the search for new antibiotics. The development of monitoring protocols is hindered by the great diversity of resistance factors, while the "streetlight effect" denies the possibility of discovering novel drugs based on existing databases. In this study, we address these challenges using high-throughput environmental screening viewed from a trait-based ecological perspective. Through an in-depth analysis of the metagenomes of 658 topsoil samples spanning Europe, we explored the distribution of 241 prokaryotic and fungal genes responsible for producing metabolites with antibiotic properties and 485 antibiotic resistance genes. We analyzed the diversity of these gene collections at different levels and modeled the distribution of each gene across environmental gradients. Our analyses revealed several nonparallel distribution patterns of the genes encoding sequential steps of enzymatic pathways synthesizing large antibiotic groups, pointing to gaps in existing databases and suggesting potential for discovering new analogues of known antibiotics. We show that agricultural activity caused a continental-scale homogenization of microbial antibiotic-related machinery, emphasizing the importance of maintaining indigenous ecosystems within the landscape mosaic. Based on the relationships between the proportion of the genes in the metagenomes with the main predictors (soil pH, land cover type, climate temperature and humidity), we illustrate how the properties of chemical structures dictate the distribution of the genes responsible for their synthesis across environments. With this understanding, we propose general principles to facilitate the discovery of antibiotics, including principally new ones, establish abundance baselines for antibiotic resistance genes, and predict their dissemination.

Soil Microbiology

Annotated genome of the Atlantic dog whelk, Nucella lapillus.

Nucella lapillus is an important player in rocky shore food chains and has been a focal organism of ecological and evolutionary studies for decades. Despite poor dispersal, they have a broad geographic range, which makes them an ideal species to examine isolation by distance and selection across environmental gradients. Here we present the fully annotated genome of N. lapillus generated with Oxford Nanopore Techonology sequencing at &#x223c;37&#xd7; coverage. The genome assembly is 2.32 Gbp and consists of 2,525 contigs, with an N50 length of 2 Mbp. Repeat annotation identified 2,491 families that cover 67.56% of the genome, which is similar to other gastropods. Despite its large size and high proportion of repeats, the genome is of high quality. Benchmarking Universal Single-Copy Ortholog (BUSCO) analysis revealed a score of 96.8%. Functional annotation of the genome produced 45,848 protein-coding genes with a 96.6% BUSCO score. Genomic resources for mollusks lag behind that of other phyla, perhaps because many of their innate characteristics complicate DNA extraction, sequencing, and assembly. This new N. lapillus genome will increase our genomic understanding of the second largest phylum (and the most diverse class within said phylum) and serve as a key resource to advance studies on the organismal biology and population genetics of this iconic species as well as the connection between genomic variation and community-level processes.

Animals

Comparisons Between Large-Scale Genomic Variants and SNPs in Driving Population Divergence and Local Adaptation.

Genomic variations, such as indels (2-49 bp) and structural variants (SVs, &#x2265;50 bp), are larger-scale mutations than single nucleotide polymorphisms (SNPs) and can substantially impact evolutionary processes, including speciation, adaptation, and phenotypes. Despite their functional importance, integrative population genetic analyses that jointly consider genome-wide SNPs, indels, and SVs remain under-explored. The ground tit (Pseudopodoces humilis), an endemic species to the Qinghai-Tibet Plateau (QTP), exhibits divergence across distinct glacial refugia, accompanied by habitat and morphological divergence, making it an excellent example for investigating how different types of genomic variants contribute to population divergence and local adaptation. Here, by retrieving 81 whole-genome sequence data, over 13 million SNPs, 2 million indels, and 22,101 SVs were identified. Variants were unevenly distributed across the genome, characterized by distinct hotspot regions. Indels and SVs revealed four genetic clusters consistent with previous SNP-based results, thereby validating the reliability of our variant datasets. FST and genotype-environment association (GEA) analyses independently revealed numerous candidate indels and SVs; each showed minimal overlap with previously identified SNPs, and were enriched in similar functional pathways such as signal transduction, skeletal muscle development, water transport, DNA repair, reproduction, nervous system development, and immunity. Collectively, our results demonstrated that indels and SVs could capture additional signatures besides SNPs. Furthermore, similar but distinct gene functions among different types of genomic variants collectively and complementarily drive genomic divergence across environmental gradients in such a high-elevation endemic species, underscoring its evolutionary relevance in local adaptation.

indels

Experimental Validation of Genome-Environment Associations in Arabidopsis.

Identifying the genetic basis of local adaptation is a key goal in evolutionary biology. Allele frequency clines along environmental gradients, known as genotype-environment associations (GEA), are often used to detect potential loci causing local adaptation but are rarely followed by experimental validation. Here, we tested loci identified in three moisture-related GEA studies on Arabidopsis. We studied 42 GEA-identified genes using t-DNA knockout lines under drought and tested effects on flowering time, an adaptive trait, and genotype-by-environment (GxE) interactions for performance and fitness. In total, 16/42 genes had significant effects on traits involved in local adaptation or performance responses to the environment. We found that wrky38 mutants had significant GxE effects for fitness; lsd1 plants had a significant GxE effect for flowering time, and 11 genes showed flowering time effects with no drought interaction. However, most GEA candidates did not exhibit GxE. In the follow-up experiments, wrky38 caused decreased stomatal conductance and specific leaf area under drought, indicating potentially adaptive drought avoidance. Additionally, GEA identified natural putative LoF variants of WRKY38 associated with dry environments, as well as alleles associated with variation in LSD1 expression. While only a few GEA-identified genes were validated for GxE interactions for fitness, we likely overlooked some genes because experiments might not well represent natural environments and t-DNA insertions might not well represent natural alleles. Nevertheless, GEAs apparently identified some genes contributing to local adaptation. GEA and follow-up experiments are straightforward to implement in model systems and demonstrate prospects for GEA discovery of new local adaptations.

Arabidopsis

Interface-dependent V. parahaemolyticus biofilm under varying temperatures, media, and oxygen conditions: implications for seafood safety.

Vibrio parahaemolyticus biofilms play a critical role in pathogen persistence in marine and seafood-processing environments, where oxygen availability, temperature, and surface interfaces vary widely. This study investigated biofilm development by three strains on partially submerged stainless-steel coupons under gas-liquid-wall (GLW) and fully submerged (SM) interfaces. Viable cell counts (log&#x2081;&#x2080;CFU/cm2) along with normalized protein concentration per viable cell (nProt) and normalized polysaccharide concentration per viable cell (nPol) were measured, under aerobic and anaerobic conditions across a temperature range of 15-30&#xa0;&#xb0;C, using tryptic soy broth with 3% NaCl (TSB) and seawater-based medium (SW). GLW biofilms consistently exhibited higher cell counts (6.4-7.3 log&#x2081;&#x2080;CFU/cm2) compared to SM biofilms (5.9-6.3 log&#x2081;&#x2080;CFU/cm2), suggesting that enhanced oxygen diffusion promotes bacterial proliferation. Conversely, SM biofilms exhibited significantly higher nProt and nPol levels (p&#xa0;<&#xa0;0.001), indicating increased production of the extracellular polymeric substance (EPS) matrix under low-oxygen, high-nutrient conditions. Microscopy and three-dimensional surface plot analyses revealed relatively uniform biofilm layers at the GLW interface, whereas SM biofilms formed heterogeneous, tower-like structures. EPS production was further influenced by medium composition, oxygen, and temperature. SM biofilms grown in SW exhibited significantly higher nProt and nPol than those in TSB under aerobic conditions (p&#xa0;<&#xa0;0.001), indicating enhanced matrix stabilization. Under anaerobic conditions at 15&#xa0;&#xb0;C, nProt and nPol were higher, whereas under aerobic conditions, peak nProt and nPol occurred at elevated temperatures. These findings highlight a trade-off between bacterial growth and matrix production and provide insight into biofilm adaptation and persistence in seafood-processing environments. These insights may help develop improved biofilm control and seafood safety management.

Biofilms

Evolutionary Reorganization of Transcriptomic Architecture Across a UVB Tolerance Gradient in Fish.

Environmental stressors such as ultraviolet radiation impose strong selective pressures on organisms, yet how adaptation to such stressors shapes transcriptomic responses at the network level remains poorly understood. Although stratospheric ozone is recovering globally, substantial regional variation in UV exposure persists, particularly in high-altitude environments where extreme UV levels can occur. Here, we compared three fish models representing distinct biological responses to UVB exposure: wild-type zebrafish (Danio rerio), a melanin-deficient zebrafish mutant (nacre) lacking a major protective mechanism against UVB damage, and the high-altitude Andean killifish Orestias ascotanensis, a species naturally exposed to extreme UVB radiation. Together, these models define a gradient spanning physiological protection, impaired protection, and evolutionary adaptation to UVB stress. Using RNA-seq and protein-protein interaction networks, we show that transcriptomic responses differ markedly across this gradient. Wild-type and nacre zebrafish exhibited relatively limited transcriptomic changes (&#x223c;2%-2.4% of genes changing), whereas O. ascotanensis displayed a large-scale and highly coordinated response (&#x223c;21.6% of genes changing) characterized by functionally specialized networks enriched in DNA repair pathways. These differences involved not only transcriptomic magnitude but also marked reorganization of transcriptomic architecture. Integration with positive selection analyses revealed that positively selected genes were concentrated within highly interconnected regions of transcriptomic networks, consistent with adaptation involving network reorganization. Furthermore, ortholog-based analyses suggest that adaptive responses involve differential reorganization of a conserved functional background. Together, our results support a model in which adaptation to environmental stress is associated with the reorganization of conserved transcriptomic networks across physiological and evolutionary contexts, providing a systems-level perspective on the molecular basis of adaptation.

UVB radiation

Isolation of plasmid deoxyribonucleic acid from Pseudomonas putida.

Conditions suitable for reproducible recovery of covalently closed circular deoxyribonucleic acid from strains of Pseudomonas putida containing degradative plasmids (CAM, SAL, OCT, etc.) have been defined. These degradative plasmids could not be isolated by the usual procedure, whereas RP1, an R factor of the P group, present in the isogenic strain of P. putida, was isolated equally well by either the usual procedure or the modified procedure. Characterization by electron microscopy of RP1 deoxyribonucleic acid confirmed the molecular weight (about 40 X 10(6)) previously determined by sucrose gradient centrifugation.

Biodegradation, Environmental

[Alveolar arterial O2 gradient in patients with cardiopulmonary pathology. Its study at rest with respiration of environmental air].

The arterial oxygen pressure (PAO2) and the arterial carbon dioxide pressure (PACO2) are the mirror of the whole stage in alveolar ventilation, because there is a numerical correlation between them, in the alveolar air equation. In our material no difference was found when the respiratory cocient is used to calculate the equation modifying the PACO2 value. On the other hand, the PaO depends on a great amount of variables, i.e., the rationship V/Q. Qs/Qt and the arteriovenous oxygen difference in volume percentage, which reflect the functional stage of the gaseous interchange; other variables depend essentially of technical factors as the methodology used in the obtention and management of the samples and the measurement of the PAO2 at the laboratory. Thus, the alveolo-arterial oxigen Difference delta (A-a)O2 is considered as a mirror of the fluctuations in gaseous interchange, only if the alveolar ventilation the cardiac output, the systemic arterial-vein oxygen difference and in minor importance the respiratory cocient (RP remain constant. The delta (A-a)O2 is not always correlated with other parameters far from the gaseous interchange at lung level in the critically ill patient, for this, it has not a pronostic mecaning by itself. Nevertheless, we believe that delta (A-a)O2 continues being a useful measurement to evaluate the gaseous interchange if at the same time all the factors which may have influence in it are analysed. We consider, for trying to be simple in the management of the acute respiratory failure at the bedside that the (A-a)O2 must not be linked to only one parameter of the respiratory function, i.e. the intrapulmonar veno-arterial shunt, without considering the cardiac output and the systemic arterial-vein oxygen difference.

Blood Pressure

Soil erosion and landscape elevation as unnoticed determinants of environmental antibiotic resistance distribution.

Climate change is reshaping the global antibiotic resistance gene (ARG) landscape through geomorphological processes that remain largely overlooked in the One Health framework. This critical review synthesises evidence on how soil erosion and landscape elevation gradients redistribute, select for, and disseminate ARGs across terrestrial and aquatic ecosystems. Erosion physically removes and transports ARG-bearing microbes, depletes nutrients, and co-selects for resistance via heavy metal exposure and horizontal gene transfer, creating source-sink dynamics that connect eroding hillslopes to downstream water bodies and food systems. Elevation gradients impose abiotic stressors-declining temperature, elevated UV radiation, and shifting pH-that drive microbial community reassembly through environmental selection and dispersal limitation, with emerging evidence linking bacterial competition at high altitude to enhanced multidrug efflux and resistome complexity. The review identifies critical knowledge gaps, including unquantified ARG mass fluxes across erosion-deposition gradients, unresolved dispersal-versus-selection mechanisms along elevation transects, and the absence of integrated One Health surveillance linking environmental ARG reservoirs to clinical outcomes. A synthesis of global case studies illustrates how these processes converge across diverse landscapes. The review concludes with a mechanistic research agenda-including reciprocal transplant experiments, landscape connectivity modelling, and cross-sectoral surveillance-needed to translate these emerging drivers into actionable climate-AMR mitigation policy.

Drug Resistance, Microbial

Environmental Stresses Constrain Soil Microbial Community Functions by Regulating Deterministic Assembly and Niche Width.

Increasing evidence indicates that the loss of soil microbial &#x3b1;-diversity triggered by environmental stress negatively impacts microbial functions; however, the effects of microbial &#x3b1;-diversity on community functions under environmental stress are poorly understood. Here, we investigated the changes in bacterial and fungal &#x3b1;- diversity along gradients of five natural stressors (temperature, precipitation, plant diversity, soil organic C and pH) across 45 grasslands in China and evaluated their connection with microbial functional traits. By quantifying the five environmental stresses into an integrated stress index, we found that the bacterial and fungal &#x3b1;-diversity declined under high environmental stress across three soil layers (0-20&#x2009;cm, 20-40&#x2009;cm and 40-60&#x2009;cm). Metagenomic-based analyses showed that the diversity of functional genes decreased along the stress gradients. High stress enhanced the abundance of genes associated with broad functional categories (e.g., glycolysis/gluconeogenesis, TCA cycle, DNA replication/repair and cell growth/death) but reduced the abundance of genes linked to specialised functional categories (e.g., C, N, S and methane metabolism). Phylogenetic null models and niche analyses indicated that stochastic assembly processes predominated in high-diversity communities, in which bacterial and fungal taxa had a narrow ecological niche. However, in low-diversity communities, deterministic assembly processes were dominant, and taxa had wide niches, correlating with the reduction in gene abundance observed for broad and specialised functional categories. Given the essential role of the microbiome in regulating ecosystem functions, our findings suggest that low-diversity-induced deterministic community assembly processes and a wide niche under high environmental stress may regulate microbial functions. These findings emphasise the ecological mechanisms through which microbial biodiversity regulates terrestrial ecosystem functioning.

Soil Microbiology

A single-nucleus and spatial transcriptomic atlas of poplar leaves reveals the regulation of leaf polarity and cuticle deposition.

Leaf adaxial-abaxial polarity is fundamental for plant morphogenesis and environmental adaptation through asymmetric cell differentiation. Emerging evidence reveals dorsoventral metabolic gradients act downstream of transcriptional networks to fine-tune cellular specialization. While conserved transcription factors (e.g., HD-ZIP III and KANADI) establish initial polarity, the molecular networks driving position-specific cellular differentiation and their integration with metabolic adaptation remain unclear. Leveraging single-nucleus and spatial transcriptomics, we resolve major cell classes (mesophyll, epidermal, and vascular-associated) and their adaxial-abaxial subtypes, revealing dorsoventral polarity in transcriptional profiles and metabolic pathways. Adaxial cells are enriched in phenylpropanoid/flavonoid biosynthesis, while abaxial cells show preferential activation of stress and hormone signaling. Notably, we identify MYC2 as a key regulator of adaxial cuticle biosynthesis, binding to promoters of lipid biosynthetic and transport genes (e.g., CER10 and LTPG1) and promoting cuticle thickening. Our study uncovers how positional identity shapes transcriptional and metabolic polarity in leaves, with MYC2 emerging as a central regulator coordinating organ-specific adaptations. These findings provide insights into the spatial regulation of plant development and stress resilience, offering potential strategies for engineering stress-tolerant woody crops.

Plant Leaves

Fungal and algal lichen symbionts show different transcriptional expression patterns in two climate zones.

In the lichen symbiosis, the fungal and algal partners constitute a closely integrated system. The combination of fungal and algal partners changes along climate gradients in many species, and is expected to be adaptive. However, the functional mechanisms behind this symbiosis-mediated environmental adaptation are unknown. We investigated which transcriptional profiles are associated with specific fungal-algal symbiont pairings found in lichens from high-elevation (Lower Supratemperate) and low-elevation (Lower Mesomediterranean) sites at two extremes of a climatic gradient on Mount Limbara, Sardinia. Using laboratory-acclimatized thalli, we found that lichen fungal and algal symbionts show variable expression profiles between high- and low-elevation individuals: circadian- and temperature-associated genes for fungi and light-responsive genes for algae show climate-specific patterns. High- and low-elevation individuals differentially express sugar transporters in both symbionts, pointing to symmetrical and climate-dependent sugar transport mechanisms between them. A light pulse treatment identified asymmetries between fungal and algal light responses, with high- and low-elevation fungal symbionts but only low-elevation algal symbionts showing a response. Together, these results tie previously observed genomic variation along climatic gradients in a lichen species to functional differences in transcription for the fungal and algal symbionts, contributing to our understanding of environmental specialization and niche-specific partner combinations in lichens.

Lichens