Adaptive evolution of the enigmatic Takakia now facing climate change in Tibet.

The most extreme environments are the most vulnerable to transformation under a rapidly changing climate. These ecosystems harbor some of the most specialized species, which will likely suffer the highest extinction rates. We document the steepest temperature increase (2010-2021) on record at altitudes of above 4,000 m, triggering a decline of the relictual and highly adapted moss Takakia lepidozioides. Its de-novo-sequenced genome with 27,467 protein-coding genes includes distinct adaptations to abiotic stresses and comprises the largest number of fast-evolving genes under positive selection. The uplift of the study site in the last 65 million years has resulted in life-threatening UV-B radiation and drastically reduced temperatures, and we detected several of the molecular adaptations of Takakia to these environmental changes. Surprisingly, specific morphological features likely occurred earlier than 165 mya in much warmer environments. Following nearly 400 million years of evolution and resilience, this species is now facing extinction.

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target.

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

Multilevel analysis between Physcomitrium patens and Mortierellaceae endophytes explores potential long-standing interaction among land plants and fungi.

The model moss species Physcomitrium patens has long been used for studying divergence of land plants spanning from bryophytes to angiosperms. In addition to its phylogenetic relationships, the limited number of differential tissues, and comparable morphology to the earliest embryophytes provide a system to represent basic plant architecture. Based on plant-fungal interactions today, it is hypothesized these kingdoms have a long-standing relationship, predating plant terrestrialization. Mortierellaceae have origins diverging from other land fungi paralleling bryophyte divergence, are related to arbuscular mycorrhizal fungi but are free-living, observed to interact with plants, and can be found in moss microbiomes globally. Due to their parallel origins, we assess here how two Mortierellaceae species, Linnemannia elongata and Benniella erionia, interact with P. patens in coculture. We also assess how Mollicute-related or Burkholderia-related endobacterial symbionts (MRE or BRE) of these fungi impact plant response. Coculture interactions are investigated through high-throughput phenomics, microscopy, RNA-sequencing, differential expression profiling, gene ontology enrichment, and comparisons among 99 other P. patens transcriptomic studies. Here we present new high-throughput approaches for measuring P. patens growth, identify novel expression of over 800 genes that are not expressed on traditional agar media, identify subtle interactions between P. patens and Mortierellaceae, and observe changes to plant-fungal interactions dependent on whether MRE or BRE are present. Our study provides insights into how plants and fungal partners may have interacted based on their communications observed today as well as identifying L. elongata and B. erionia as modern fungal endophytes with P. patens.

Drying without dying: A genome database for desiccation-tolerant plants and evolution of desiccation tolerance.

Desiccation is typically fatal, but a small number of land plants have evolved vegetative desiccation tolerance (VDT), allowing them to dry without dying through a process called anhydrobiosis. Advances in sequencing technologies have enabled the investigation of genomes for desiccation-tolerant plants over the past decade. However, a dedicated and integrated database for these valuable genomic resources has been lacking. Our prolonged interest in VDT plant genomes motivated us to create the "Drying without Dying" database, which contains a total of 16 VDT-related plant genomes (including 10 mosses) and incorporates 10 genomes that are closely related to VDT plants. The database features bioinformatic tools, such as blast and homologous cluster search, sequence retrieval, Gene Ontology term and metabolic pathway enrichment statistics, expression profiling, co-expression network extraction, and JBrowser exploration for each genome. To demonstrate its utility, we conducted tailored PFAM family statistical analyses, and we discovered that the drought-responsive ABA transporter AWPM-19 family is significantly tandemly duplicated in all bryophytes but rarely so in tracheophytes. Transcriptomic investigations also revealed that response patterns following desiccation diverged between bryophytes and angiosperms. Combined, the analyses provided genomic and transcriptomic evidence supporting a possible divergence and lineage-specific evolution of VDT in plants. The database can be accessed at http://desiccation.novogene.com. We expect this initial release of the "Drying without Dying" plant genome database will facilitate future discovery of VDT genetic resources.

Decline in biological soil crust N-fixing lichens linked to increasing summertime temperatures.

Biological soil crusts (biocrusts), comprised of mosses, lichens, and cyanobacteria, are key components to many dryland communities. Climate change and other anthropogenic disturbances are thought to cause a decline in mosses and lichens, yet few long-term studies exist to track potential shifts in these sensitive soil-surface communities. Using a unique long-term observational dataset from a temperate dryland with initial observations dating back to 1967, we examine the effects of 53 y of observed environmental variation and Bromus tectorum invasion on biocrust communities in a grassland never grazed by domestic livestock. Annual observations show a steep decline in N-fixing lichen cover (dominated by Collema species) from 1996 to 2002, coinciding with a period of extended drought, with Collema communities never able to recover. Declines in other lichen species were also observed, both in number of species present and by total cover, which were attributed to increasing summertime temperatures. Conversely, moss species gradually gained in cover over the survey years, especially following a large Bromus tectorum invasion at the study onset (ca. 1996 to 2001). These results support a growing body of studies that suggests climate change is a key driver in changes to certain components of late-successional biocrust communities. Results here suggest that warming may partially negate decades of protection from disturbance, with biocrust communities reaching a vital tipping point. The accelerated rate of ongoing warming observed in this study may have resulted in the loss of lichen cover and diversity, which could have long-term implications for global temperate dryland ecosystems.

An optimized transformation protocol for Anthoceros agrestis and three more hornwort species.

Land plants comprise two large monophyletic lineages, the vascular plants and the bryophytes, which diverged from their most recent common ancestor approximately 480 million years ago. Of the three lineages of bryophytes, only the mosses and the liverworts are systematically investigated, while the hornworts are understudied. Despite their importance for understanding fundamental questions of land plant evolution, they only recently became amenable to experimental investigation, with Anthoceros agrestis being developed as a hornwort model system. Availability of a high-quality genome assembly and a recently developed genetic transformation technique makes A. agrestis an attractive model species for hornworts. Here we describe an updated and optimized transformation protocol for A. agrestis, which can be successfully used to genetically modify one more strain of A. agrestis and three more hornwort species, Anthoceros punctatus, Leiosporoceros dussii, and Phaeoceros carolinianus. The new transformation method is less laborious, faster, and results in the generation of greatly increased numbers of transformants compared with the previous method. We have also developed a new selection marker for transformation. Finally, we report the development of a set of different cellular localization signal peptides for hornworts providing new tools to better understand the hornwort cell biology.

Horizontally acquired fungal killer protein genes affect cell development in mosses.

The moss Physcomitrium patens is crucial for studying plant development and evolution. Although the P. patens genome includes genes acquired from bacteria, fungi and viruses, the functions and evolutionary significance of these acquired genes remain largely unclear. Killer protein 4 (KP4) is a toxin secreted by the phytopathogenic fungus Ustilago maydis that inhibits the growth of sensitive target strains by blocking their calcium uptake. Here, we show that KP4 genes in mosses were acquired from fungi through at least three independent events of horizontal gene transfer. Two paralogous copies of KP4 (PpKP4-1 and PpKP4-2) exist in P. patens. Knockout mutants ppkp4-1 and ppkp4-2 showed cell death at the protonemal stage, and ppkp4-2 also exhibited defects in tip growth. We provide experimental evidence indicating that PpKP4-1/2 affects P. patens protonemal cell development by mediating cytoplasmic calcium and that KP4 genes are functionally conserved between P. patens and fungi. The present study provides additional insights into the role of horizontal gene transfer in land plant development and evolution.

The ScAPD1-like gene from the desert moss Syntrichia caninervis enhances resistance to Verticillium dahliae via phenylpropanoid gene regulation.

Soloist is a member of a distinct and small subfamily within the AP2/ERF transcriptional factor family that play important roles in plant biotic and abiotic stress responses. There are limited studies of Soloist genes and their functions are poorly understood. We characterized the abiotic and biotic stress tolerance function of the ScSoloist gene (designated as ScAPD1-like) from the desert moss Syntrichia caninervis. ScAPD1-like responded to multiple abiotic, biotic stresses and plant hormone treatments. ScAPD1-like protein located to the nucleus and bound to several DNA elements. Overexpression of ScAPD1-like in Arabidopsis did not alter abiotic stress resistance or inhibit Pseudomonas syringae pv. tomato (Pst) DC3000 infection. However, overexpression of ScAPD1-like significantly increased the resistance of transgenic Arabidopsis and S. caninervis to Verticillium dahliae infection, decreased reactive oxygen species accumulation and improved reactive oxygen species scavenging activity. ScAPD1-like overexpression plants altered the abundance of transcripts for lignin synthesis and promoted lignin accumulation in Arabidopsis. ScAPD1-like directly bind to RAV1, AC elements, and TATA-box in the promoters of AtPAL1 and AtC4H genes, respectively, in vitro. Chromatin immunoprecipitation-quantitative polymerase chain reaction assays demonstrated ScAPD1-like directly bound to PAL and C4H genes promoters in Arabidopsis and their homologs in S. caninervis. In S. caninervis, ScAPD1-like overexpression and RNAi directly regulated the abundance of ScPAL and ScC4H transcripts and modified the metabolites of phenylpropanoid pathway. We provide insight into the function of Soloist in plant defense mechanisms that likely occurs through activation of the phenylpropanoid biosynthesis pathway. ScAPD1-like is a promising candidate gene for breeding strategies to improve resistance to Verticillium wilt.

PpSCARECROW1 (PpSCR1) regulates leaf blade and mid-vein development in Physcomitrium patens.

In plants, asymmetric cell divisions result in distinct cell fates forming large and small daughter cells, adding to the cellular diversity in an organ. SCARECROW (SCR), a GRAS domain-containing transcription factor controls asymmetric periclinal cell divisions in flowering plants by governing radial patterning of ground tissue in roots and cell proliferation in leaves. Though SCR homologs are present across land plant lineages, the current understanding of their role in cellular patterning and leaf development is mostly limited to flowering plants. Our phylogenetic analysis identified three SCR homologs in moss Physcomitrium patens, amongst which PpSCR1 showed highest expression in gametophores and its promoter activity was prominent at the mid-vein and the flanking leaf blade cells pointing towards its role in leaf development. Notably, out of the three SCR homologs, only the ppscr1 knock-out lines developed slender leaves with four times narrower leaf blade and three times thicker mid-vein. Detailed histology studies revealed that slender leaf phenotype is either due to the loss of anticlinal cell divisions or failure of periclinal division suppression in the leaf blade. RNA-Seq analyses revealed that genes responsible for cell division and differentiation are expressed differentially in the mutant. PpSCR1 overexpression lines exhibited significantly wider leaf lamina, further reconfirming the role in leaf development. Together, our data suggests that PpSCR1 is involved in the leaf blade and mid-vein development of moss and that its role in the regulation of cell division and proliferation is ancient and conserved among flowering plants and mosses.

Biological nitrogen fixation, diversity and community structure of diazotrophs in two mosses in 25 temperate forests.

Many moss species are associated with nitrogen (N)-fixing bacteria (diazotrophs) that support the N supply of mosses. Our knowledge relates primarily to pristine ecosystems with low atmospheric N input, but knowledge of biological N fixation (BNF) and diazotrophic communities in mosses in temperate forests with high N deposition is limited. We measured BNF rates using the direct stable isotope method and studied the total and potentially active diazotrophic communities in two abundant mosses, Brachythecium rutabulum and Hypnum cupressiforme, both growing on lying deadwood trunks in 25 temperate forest sites. BNF rates in both mosses were similar to those observed in moss species of pristine ecosystems. H. cupressiforme fixed three times more N2 and exhibited lower diazotrophic richness than B. rutabulum. Frankia was the most prominent diazotroph followed by cyanobacteria Nostoc. Manganese, iron, and molybdenum contents in mosses were positively correlated with BNF and diazotrophic communities. Frankia maintained high BNF rates in H. cupressiforme and B. rutabulum even under high chronic N deposition in Central European forests. Moss N concentration and 15 N abundance indicate a rather minor contribution of BNF to the N nutrition of these mosses.

Habitat conditions and not moss composition mediate microbial community structure in Swiss peatlands.

Peatlands, one of the oldest ecosystems, globally store significant amounts of carbon and freshwater. However, they are under severe threat from human activities, leading to changes in water, nutrient and temperature regimes in these delicate systems. Such shifts can trigger a substantial carbon flux into the atmosphere and diminish the water-holding capacity of peatlands. Microbes associated with moss in peatlands play a crucial role in providing these ecosystem services, which are at risk due to global change. Therefore, understanding the factors influencing microbial composition and function is vital. Our study focused on five peatlands along an altitudinal gradient in Switzerland, where we sampled moss on hummocks containing Sarracenia purpurea. Structural equation modelling revealed that habitat condition was the primary predictor of community structure and directly influenced other environmental variables. Interestingly, the microbial composition was not linked to the local moss species identity. Instead, microbial communities varied significantly between sites due to differences in acidity levels and nitrogen availability. This finding was also mirrored in a co-occurrence network analysis, which displayed a distinct distribution of indicator species for acidity and nitrogen availability. Therefore, peatland conservation should take into account the critical habitat characteristics of moss-associated microbial communities.

Evolution in the Cycles of Life.

The life cycles of eukaryotes alternate between haploid and diploid phases, which are initiated by meiosis and gamete fusion, respectively. In both ascomycete and basidiomycete fungi and chlorophyte algae, the haploid-to-diploid transition is regulated by a pair of paralogous homeodomain protein encoding genes. That a common genetic program controls the haploid-to-diploid transition in phylogenetically disparate eukaryotic lineages suggests this may be the ancestral function for homeodomain proteins. Multicellularity has evolved independently in many eukaryotic lineages in either one or both phases of the life cycle. Organisms, such as land plants, exhibiting a life cycle whereby multicellular bodies develop in both the haploid and diploid phases are often referred to as possessing an alternation of generations. We review recent progress on understanding the genetic basis for the land plant alternation of generations and highlight the roles that homeodomain-encoding genes may have played in the evolution of complex multicellularity in this lineage.

Evolution of parental care in haploid-diploid plants.

In bryophytes that alternate between haploid gametophytes and diploid sporophytes through sexual reproduction, sporophytes are often attached to and nurtured on the female gametophyte. A similar phenomenon is seen in Florideophyceae (a group of red algae). These systems in which a gametophyte (mother) invests nutrients in sporophytes (offspring) are ideal for studying the evolution of 'parental care' in non-animal organisms. Here, we propose a model of a haploid-diploid life cycle and examine the evolution of maternal investment in sporophytes focusing on two effects: the degree of paternal or maternal control of investment and the number of sporophytes. We find that when the female dominantly controls the investment, the evolutionarily stable level of investment is that which maximizes the expected reproductive success of the female gametophyte. The genomic conflict between maternal and paternal alleles complicates the evolutionary outcome; however, a greater male allelic effect and a higher number of sporophytes favour a higher energy investment, which may lead to evolutionary branching or run-away escalation of the investment level. This suggests that the selfishness of the paternal gene is the evolutionary driver of parental care and that complex structures such as fusion cells in red algae may have evolved to suppress it.

Plasmodesmata dynamics in bryophyte model organisms: secondary formation and developmental modifications of structure and function.

MAIN CONCLUSION: Developing bryophytes differentially modify their plasmodesmata structure and function. Secondary plasmodesmata formation via twinning appears to be an ancestral trait. Plasmodesmata networks in hornwort sporophyte meristems resemble those of angiosperms. All land-plant taxa use plasmodesmata (PD) cell connections for symplasmic communication. In angiosperm development, PD networks undergo an extensive remodeling by structural and functional PD modifications, and by postcytokinetic formation of additional secondary PD (secPD). Since comparable information on PD dynamics is scarce for the embryophyte sister groups, we investigated maturating tissues of Anthoceros agrestis (hornwort), Physcomitrium patens (moss), and Marchantia polymorpha (liverwort). As in angiosperms, quantitative electron microscopy revealed secPD formation via twinning in gametophytes of all model bryophytes, which gives rise to laterally adjacent PD pairs or to complex branched PD. This finding suggests that PD twinning is an ancient evolutionary mechanism to adjust PD numbers during wall expansion. Moreover, all bryophyte gametophytes modify their existing PD via taxon-specific strategies resembling those of angiosperms. Development of type II-like PD morphotypes with enlarged diameters or formation of pit pairs might be required to maintain PD transport rates during wall thickening. Similar to angiosperm leaves, fluorescence redistribution after photobleaching revealed a considerable reduction of the PD permeability in maturating P. patens phyllids. In contrast to previous reports on monoplex meristems of bryophyte gametophytes with single initials, we observed targeted secPD formation in the multi-initial basal meristems of A. agrestis sporophytes. Their PD networks share typical features of multi-initial angiosperm meristems, which may hint at a putative homologous origin. We also discuss that monoplex and multi-initial meristems may require distinct types of PD networks, with or without secPD formation, to control maintenance of initial identity and positional signaling.

Characterization of the fiber-like cortical cells in moss gametophytes.

MAIN CONCLUSION: Fiber-like cells with thickened cell walls of specific structure and polymer composition that includes (1 → 4)-ß-galactans develop in the outer stem cortex of several moss species gametophytes. The early land plants evolved several specialized cell types and tissues that did not exist in their aquatic ancestors. Of these, water-conducting elements and reproductive organs have received most of the research attention. The evolution of tissues specialized to fulfill a mechanical function is by far less studied despite their wide distribution in land plants. For vascular plants following a homoiohydric trajectory, the evolutionary emergence of mechanical tissues is mainly discussed starting with the fern-like plants with their hypodermal sterome or sclerified fibers that have xylan and lignin-based cell walls. However, mechanical challenges were also faced by bryophytes, which lack lignified cell-walls. To characterize mechanical tissues in the bryophyte lineage, following a poikilohydric trajectory, we used six wild moss species (Polytrichum juniperinum, Dicranum sp., Rhodobryum roseum, Eurhynchiadelphus sp., Climacium dendroides, and Hylocomium splendens) and analyzed the structure and composition of their cell walls. In all of them, the outer stem cortex of the leafy gametophytic generation had fiber-like cells with a thickened but non-lignified cell wall. Such cells have a spindle-like shape with pointed tips. The additional thick cell wall layer in those fiber-like cells is composed of sublayers with structural evidence for different cellulose microfibril orientation, and with specific polymer composition that includes (1 → 4)-ß-galactans. Thus, the basic cellular characters of the cells that provide mechanical support in vascular plant taxa (elongated cell shape, location at the periphery of a primary organ, the thickened cell wall and its peculiar composition and structure) also exist in mosses.

Paternal imprinting in Marchantia polymorpha.

We are becoming aware of a growing number of organisms that do not express genetic information equally from both parents as a result of an epigenetic phenomenon called genomic imprinting. Recently, it was shown that the entire paternal genome is repressed during the diploid phase of the life cycle of the liverwort Marchantia polymorpha. The deposition of the repressive epigenetic mark H3K27me3 on the male pronucleus is responsible for the imprinted state, which is reset by the end of meiosis. Here, we put these recent reports in perspective of other forms of imprinting and discuss the potential mechanisms of imprinting in bryophytes and the causes of its evolution.

Impact of changing climate on bryophyte contributions to terrestrial water, carbon, and nitrogen cycles.

Bryophytes, including the lineages of mosses, liverworts, and hornworts, are the second-largest photoautotroph group on Earth. Recent work across terrestrial ecosystems has highlighted how bryophytes retain and control water, fix substantial amounts of carbon (C), and contribute to nitrogen (N) cycles in forests (boreal, temperate, and tropical), tundra, peatlands, grasslands, and deserts. Understanding how changing climate affects bryophyte contributions to global cycles in different ecosystems is of primary importance. However, because of their small physical size, bryophytes have been largely ignored in research on water, C, and N cycles at global scales. Here, we review the literature on how bryophytes influence global biogeochemical cycles, and we highlight that while some aspects of global change represent critical tipping points for survival, bryophytes may also buffer many ecosystems from change due to their capacity for water, C, and N uptake and storage. However, as the thresholds of resistance of bryophytes to temperature and precipitation regime changes are mostly unknown, it is challenging to predict how long this buffering capacity will remain functional. Furthermore, as ecosystems shift their global distribution in response to changing climate, the size of different bryophyte-influenced biomes will change, resulting in shifts in the magnitude of bryophyte impacts on global ecosystem functions.

Divergent evolution of the alcohol-forming pathway of wax biosynthesis among bryophytes.

The plant cuticle is a hydrophobic barrier, which seals the epidermal surface of most aboveground organs. While the cuticle biosynthesis of angiosperms has been intensively studied, knowledge about its existence and composition in nonvascular plants is scarce. Here, we identified and characterized homologs of Arabidopsis thaliana fatty acyl-CoA reductase (FAR) ECERIFERUM 4 (AtCER4) and bifunctional wax ester synthase/acyl-CoA:diacylglycerol acyltransferase 1 (AtWSD1) in the liverwort Marchantia polymorpha (MpFAR2 and MpWSD1) and the moss Physcomitrium patens (PpFAR2A, PpFAR2B, and PpWSD1). Although bryophyte harbor similar compound classes as described for angiosperm cuticles, their biosynthesis may not be fully conserved between the bryophytes M. polymorpha and P. patens or between these bryophytes and angiosperms. While PpFAR2A and PpFAR2B contribute to the production of primary alcohols in P. patens, loss of MpFAR2 function does not affect the wax profile of M. polymorpha. By contrast, MpWSD1 acts as the major wax ester-producing enzyme in M. polymorpha, whereas mutations of PpWSD1 do not affect the wax ester levels of P. patens. Our results suggest that the biosynthetic enzymes involved in primary alcohol and wax ester formation in land plants have either evolved multiple times independently or undergone pronounced radiation followed by the formation of lineage-specific toolkits.

Life under the snow: A year-round transcriptome analysis of Antarctic mosses in natural habitats provides insight into the molecular adaptation of plants under extreme environment.

Mosses are vital components of ecosystems, exhibiting remarkable adaptability across diverse habitats from deserts to polar ice caps. Sanionia uncinata (Hedw.) Loeske, a dominant Antarctic moss survives extreme environmental condition through perennial lifecycles involving growth and dormancy alternation. This study explores genetic controls and molecular mechanisms enabling S. uncinata to cope with seasonality of the Antarctic environment. We analysed the seasonal transcriptome dynamics of S. uncinata collected monthly from February 2015 to January 2016 in King George Island, Antarctica. Findings indicate that genes involved in plant growth were predominantly upregulated in Antarctic summer, while those associated with protein synthesis and cell cycle showed marked expression during the winter-to-summer transition. Genes implicated in cellular stress and abscisic acid signalling were highly expressed in winter. Further, validation included a comparison of the Antarctic field transcriptome data with controlled environment simulation of Antarctic summer and winter temperatures, which revealed consistent gene expression patterns in both datasets. This proposes a seasonal gene regulatory model of S. uncinate to understand moss adaptation to extreme environments. Additionally, this data set is a valuable resource for predicting genetic responses to climatic fluctuations, enhancing our knowledge of Antarctic flora's resilience to global climate change.

Enhanced nucleotide analysis enables the quantification of deoxynucleotides in plants and algae revealing connections between nucleoside and deoxynucleoside metabolism.

Detecting and quantifying low-abundance (deoxy)ribonucleotides and (deoxy)ribonucleosides in plants remains difficult; this is a major roadblock for the investigation of plant nucleotide (NT) metabolism. Here, we present a method that overcomes this limitation, allowing the detection of all deoxy- and ribonucleotides as well as the corresponding nucleosides from the same plant sample. The method is characterized by high sensitivity and robustness enabling the reproducible detection and absolute quantification of these metabolites even if they are of low abundance. Employing the new method, we analyzed Arabidopsis thaliana null mutants of CYTIDINE DEAMINASE, GUANOSINE DEAMINASE, and NUCLEOSIDE HYDROLASE 1, demonstrating that the deoxyribonucleotide (dNT) metabolism is intricately interwoven with the catabolism of ribonucleosides (rNs). In addition, we discovered a function of rN catabolic enzymes in the degradation of deoxyribonucleosides in vivo. We also determined the concentrations of dNTs in several mono- and dicotyledonous plants, a bryophyte, and three algae, revealing a correlation of GC to AT dNT ratios with genomic GC contents. This suggests a link between the genome and the metabolome previously discussed but not experimentally addressed. Together, these findings demonstrate the potential of this new method to provide insight into plant NT metabolism.