Circadian clock- and temperature-associated genes contribute to overall genomic differentiation along elevation in lichenized fungi.

Circadian regulation is linked to local environmental adaptation, and many species with broad climatic niches display variation in circadian genes. Here, we hypothesize that lichenizing fungi occupying different climate zones tune their metabolism to local environmental conditions with the help of their circadian systems. We study two species of the genus Umbilicaria occupying similar climatic niches (Mediterranean and the cold temperate) in different continents. Using homology to Neurospora crassa genes, we identify gene sets associated with circadian rhythms (11 core, 39 peripheral genes) as well as temperature response (37 genes). Nucleotide diversity of these genes is significantly correlated with mean annual temperature, minimum temperature of the coldest month and mean temperature of the coldest quarter. Furthermore, we identify altitudinal clines in allele frequencies in several non-synonymous substitutions in core clock components, for example, white collar-like, frh-like and various ccg-like genes. A dN/dS approach revealed a few significant peripheral clock- and temperature-associated genes (e.g. ras-1-like, gna-1-like) that may play a role in fine-tuning the circadian clock and temperature-response machinery. An analysis of allele frequency changes demonstrated the strongest evidence for differentiation above the genomic background in the clock-associated genes in U. pustulata. These results highlight the likely relevance of the circadian clock in environmental adaptation, particularly frost tolerance, of lichens. Whether or not the fungal clock modulates the symbiotic interaction within the lichen consortium remains to be investigated. We corroborate the finding of genetic variation in clock components along altitude-not only latitude-as has been reported in other species.

Gene abundance linked to climate zone: Parallel evolution of gene content along elevation gradients in lichenized fungi.

Introduction: Intraspecific genomic variability affects a species' adaptive potential toward climatic conditions. Variation in gene content across populations and environments may point at genomic adaptations to specific environments. The lichen symbiosis, a stable association of fungal and photobiont partners, offers an excellent system to study environmentally driven gene content variation. Many of these species have remarkable environmental tolerances, and often form populations across different climate zones. Here, we combine comparative and population genomics to assess the presence and absence of genes in high and low elevation genomes of two lichenized fungi of the genus Umbilicaria. Methods: The two species have non-overlapping ranges, but occupy similar climatic niches in North America (U. phaea) and Europe (U. pustulata): high elevation populations are located in the cold temperate zone and low elevation populations in the Mediterranean zone. We assessed gene content variation along replicated elevation gradients in each of the two species, based on a total of 2050 individuals across 26 populations. Specifically, we assessed shared orthologs across species within the same climate zone, and tracked, which genes increase or decrease in abundance within populations along elevation. Results: In total, we found 16 orthogroups with shared orthologous genes in genomes at low elevation and 13 at high elevation. Coverage analysis revealed one ortholog that is exclusive to genomes at low elevation. Conserved domain search revealed domains common to the protein kinase superfamily. We traced the discovered ortholog in populations along five replicated elevation gradients on both continents and found that the number of this protein kinase gene linearly declined in abundance with increasing elevation, and was absent in the highest populations. Discussion: We consider the parallel loss of an ortholog in two species and in two geographic settings a rare find, and a step forward in understanding the genomic underpinnings of climatic tolerances in lichenized fungi. In addition, the tracking of gene content variation provides a widely applicable framework for retrieving biogeographical determinants of gene presence/absence patterns. Our work provides insights into gene content variation of lichenized fungi in relation to climatic gradients, suggesting a new research direction with implications for understanding evolutionary trajectories of complex symbioses in relation to climatic change.

Identification and expression of functionally conserved circadian clock genes in lichen-forming fungi.

Lichen-forming fungi establish stable symbioses with green algae or cyanobacteria. Many species have broad distributions, both in geographic and ecological space, making them ideal subjects to study organism-environment interactions. However, little is known about the specific mechanisms that contribute to environmental adaptation in lichen-forming fungi. The circadian clock provides a well-described mechanism that contributes to regional adaptation across a variety of species, including fungi. Here, we identify the putative circadian clock components in phylogenetically divergent lichen-forming fungi. The core circadian genes (frq, wc-1, wc-2, frh) are present across the Fungi, including 31 lichen-forming species, and their evolutionary trajectories mirror overall fungal evolution. Comparative analyses of the clock genes indicate conserved domain architecture among lichen- and non-lichen-forming taxa. We used RT-qPCR to examine the core circadian loop of two unrelated lichen-forming fungi, Umbilicaria pustulata (Lecanoromycetes) and Dermatocarpon miniatum (Eurotiomycetes), to determine that the putative frq gene is activated in a light-dependent manner similar to the model fungus Neurospora crassa. Together, these results demonstrate that lichen-forming fungi retain functional light-responsive mechanisms, including a functioning circadian clock. Our findings provide a stepping stone into investigating the circadian clock in the lichen symbiosis, e.g. its role in adaptation, and in synchronizing the symbiotic interaction.

A Candidate Gene Cluster for the Bioactive Natural Product Gyrophoric Acid in Lichen-Forming Fungi.

Natural products of lichen-forming fungi are structurally diverse and have a variety of medicinal properties. Despite this, they have limited implementation in industry mostly because the corresponding genes are unknown for most of their natural products. Here, we implement a long-read sequencing and bioinformatic approach to identify the putative biosynthetic gene cluster of the bioactive natural product gyrophoric acid (GA). Using 15 high-quality genomes representing nine GA-producing species of the lichen-forming fungal genus Umbilicaria, we identify the most likely GA cluster and investigate the cluster gene organization and composition across the nine species. Our results show that GA clusters are promiscuous within Umbilicaria, and only three genes are conserved across species, including the polyketide synthase (PKS) gene. In addition, our results suggest that the same cluster codes for different, but structurally similar compounds, namely, GA, umbilicaric-, and hiascic acid, bringing new evidence that lichen metabolite diversity is also generated through regulatory mechanisms at the molecular level. Ours is the first study to identify the most likely GA cluster and, thus, provides essential information to open new avenues for biotechnological approaches to producing and modifying GA and similar lichen-derived compounds. GA PKS is the first tridepside PKS to be identified. IMPORTANCE The implementation of natural products in the pharmaceutical industry relies on the possibility of modifying the natural product (NP) pathway to optimize yields and pharmacological effects. Characterization of genes and pathways underlying natural product biosynthesis is a major bottleneck for exploiting the medicinal properties of the natural products. Genome mining is a promising and relatively cost- and time-effective approach to utilize unexplored NP resources for drug discovery. In this study, we identify the most likely gene cluster for the lichen-forming fungal depside gyrophoric acid in nine Umbilicaria species. This compound shows cytotoxic and antiproliferative properties against several cancer cell lines and is also a broad-spectrum antimicrobial agent. This information paves the way for generating GA analogs with modified properties by selective activation/deactivation of genes.

TOC1 in Nicotiana attenuata regulates efficient allocation of nitrogen to defense metabolites under herbivory stress.

The circadian clock contextualizes plant responses to environmental signals. Plants use temporal information to respond to herbivory, but many of the functional roles of circadian clock components in these responses, and their contribution to fitness, remain unknown. We investigate the role of the central clock regulator TIMING OF CAB EXPRESSION 1 (TOC1) in Nicotiana attenuata's defense responses to the specialist herbivore Manduca sexta under both field and glasshouse conditions. We utilize 15 N pulse-labeling to quantify nitrogen incorporation into pools of three defense compounds: caffeoylputrescine (CP), dicaffeoyl spermidine (DCS) and nicotine. Nitrogen incorporation was decreased in CP and DCS and increased in nicotine pools in irTOC1 plants compared to empty vector (EV) under control conditions, but these differences were abolished after simulated herbivory. Differences between EV and irTOC1 plants in nicotine, but not phenolamide production, were abolished by treatment with the ethylene agonist 1-methylcyclopropene. Using micrografting, TOC1's effect on nicotine was isolated to the root and did not affect the fitness of heterografts under field conditions. These results suggest that the circadian clock contributes to plant fitness by balancing production of metabolically expensive nitrogen-rich defense compounds and mediating the allocation of resources between vegetative biomass and reproduction.

Determining the scale at which variation in a single gene changes population yields.

Plant trait diversity is known to influence population yield, but the scale at which this happens remains unknown: divergent individuals might change yields of immediate neighbors (neighbor scale) or of plants across a population (population scale). We use Nicotiana attenuata plants silenced in mitogen-activated protein kinase 4 (irMPK4) - with low water-use efficiency (WUE) - to study the scale at which water-use traits alter intraspecific population yields. In the field and glasshouse, we observed overyielding in populations with low percentages of irMPK4 plants, unrelated to water-use phenotypes. Paired-plant experiments excluded the occurrence of overyielding effects at the neighbor scale. Experimentally altering field arbuscular mycorrhizal fungal associations by silencing the Sym-pathway gene NaCCaMK did not affect reproductive overyielding, implicating an effect independent of belowground AMF interactions. Additionally, micro-grafting experiments revealed dependence on shoot-expressed MPK4 for N. attenuata to vary its yield per neighbor presence. We find that variation in a single gene, MPK4, is responsible for population overyielding through a mechanism, independent of irMPK4's WUE phenotype, at the aboveground, population scale.

Whether on farmland or in a forest, plants do not grow in isolation. Plants compete with their neighbors over limited space and resources, and individual plants respond to this competition in different ways by changing how much they grow and how they use resources. The efficiency with which crop plants use water, for example, is one trait that is dramatically influenced by neighboring plants and is of increasing concern given the warming climate. Understanding the effects of interactions between individual plants in a population as a whole is complicated, especially in natural plant communities where neighbors are often from different species. For this reason, McGale et al. took a different approach and looked at neighbors that were all from the same species and differed only in the activity of a single gene. The species in question was coyote tobacco, a plant that is native to western North America. McGale et al. used genetic engineering to silence a gene called MPK4, which was known from previous studies to have the effect of reducing water-use efficiency. Some of these 'water-inefficient' plants were then grown in mixed populations with plants that had normal levels of MPK4. In experiments conducted both in a glasshouse and at a field station in the Utah desert, McGale et al. found that populations with a low percentage of the MPK4-silenced plants were actually more productive than 'monocultures' that were all one type or the other. Further analysis showed that the increase in productivity did not depend on the different soil nutrient or water use of the different populations, or even the density of the plants in the populations. Pairs of plants grown in single pots essentially ruled out any interactions between immediate neighbors being responsible for the increased productivity, suggesting that that effect must instead emerge at the level of the population. Perhaps unexpectedly, McGale et al. also found that the MPK4-silenced plants and control plants did not actually differ in how they used water when grown in the field (previous studies had all been conducted in glasshouses), indicating that this trait also could not explain the observed population-level effect. Finally, experiments that involved grafting the shoots of one plant onto the roots of another suggested that the effect most likely comes from the aboveground parts of the plant. Ecologists have previously noted that more diverse populations typically have higher productivity. This new finding that a small percentage of slightly different plants in an otherwise uniform population can increase overall productivity will likely to be of special interest to researchers looking to boost the efficiency of agricultural ecosystems. Also, since MPK4 is highly conserved, and thus likely to be found in many plant species, this could be an interesting trait with which to study the interactions of natural plant communities.

The Clock Gene TOC1 in Shoots, Not Roots, Determines Fitness of Nicotiana attenuata under Drought.

The highly conserved core circadian clock component TIMING OF CAB EXPRESSION1 (TOC1) contextualizes environmental stress responses in plants, for example by gating abscisic acid signaling and suppressing thermoresponsive growth. Selective interaction of TOC1 with PHYTOCHROME B under far-red-enriched light suggests a connection between circadian gating of light responses and sensitivity to ABA, an important regulator of growth and stress responses, including under drought. However, the fitness consequences of TOC1 function, particularly in the root, are poorly understood. Here, we used the desert annual, Nicotiana attenuata, to investigate the function of TOC1 in shoots and roots for maintaining fitness under drought, in both field and glasshouse experiments. Despite marked decreases in leaf water loss, TOC1-deficient lines failed to maintain fitness in response to drought stress as measured by total seed capsule production. Restoring TOC1 transcript levels in shoots via micrografting was sufficient to restore wild-type drought responses under field conditions. Microarrays identified a coexpression module in leaves strongly linking red and far-red light signaling to drought responses in a TOC1-dependent manner, but experiments with phytochrome-deficient lines revealed that the effects of TOC1 deficiency under drought cannot be attributed to changes in red/far-red light perception alone. Taken together, these results elucidate the sophisticated, tissue-dependent role of the circadian clock in maintaining fitness in the face of long-term abiotic stresses such as drought.