The rate and molecular spectrum of spontaneous mutations in Arabidopsis thaliana

Population Genomics Reveals Demographic History and Climate Adaptation in Japanese Arabidopsis halleri

Climate oscillations in the Quaternary forced species to major latitudinal or altitudinal range shifts. It has been suggested that adaptation concomitant with range shifts plays key roles in species responses during climate oscillations, but the role of selection for local adaptation to climatic changes remains largely unexplored. Here, we investigated population structure, demographic history and signatures of climate-driven selection based on genome-wide polymorphism data of 141 Japanese Arabidopsis halleri individuals, with European ones as outgroups. Coalescent-based analyses suggested a genetic differentiation between Japanese subpopulations since the Last Glacial Period (LGP), which would have contributed to shaping the current pattern of population structure. Population demographic analysis revealed the population size fluctuations in the LGP, which were particularly prominent since the subpopulations started to diverge (∼50, 000 years ago). The ecological niche modeling predicted the geographic or distribution range shifts from southern coastal regions to northern coastal and mountainous areas, possibly in association with the population size fluctuations. Through genome-wide association analyses of bioclimatic variables and selection scans, we investigated whether climate-associated loci are enriched in the extreme tails of selection scans, and demonstrated the prevailing signatures of selection, particularly toward a warmer climate in southern subpopulations and a drier environment in northern subpopulations, which may have taken place during or after the LGP. Our study highlights the importance of integrating climate associations, selection scans and population demographic analyses for identifying genomic signatures of population-specific adaptation, which would also help us predict the evolutionary responses to future climate changes.

Two haplotype-resolved telomere-to-telomere genome assemblies of Xanthoceras sorbifolium

Yellowhorn (Xanthoceras sorbifolium) is widely used in northern China for landscaping, desertification control, and oil production. However, the lack of high-quality genomes has hindered breeding and evolutionary studies. Here, we present the first haplotype-resolved, telomere-to-telomere (T2T) yellowhorn genomes of PBN-43 (white single-flowered) and PBN-126 (white double-flowered) using PacBio HiFi and Hi-C data. These assemblies range from 464.34 Mb to 468.97 Mb and include all centromeres and telomeres. Genome annotation revealed that an average of 67.99% (317.09 Mb) of yellowhorn genomic regions consist of repetitive elements across all haplotypes. The number of protein-coding genes ranges from 35,039 to 35,174 among assemblies, representing an average 50.16% increase over the first published yellowhorn genome. Additionally, 93.90% of the annotated genes have functional annotations. We found yellowhorn experienced an LTR-RT burst during the last 0.45-0.48 Mya. These data provide a resource for investigating genomic variations, phylogenetic relationships, duplication modes, and the distribution of nucleotide-binding leucine-rich repeat (NLR) genes, and support further research into yellowhorn breeding.

Genomic signature and evolutionary history of completely cleistogamous lineages in the non-photosynthetic orchid Gastrodia

Despite a long-standing interest since Darwin's time, the genomic implications of obligate self-fertilization remain elusive. Complete cleistogamy-the obligate production of closed, self-pollinating flowers-represents an extreme reproductive strategy. Here, we present the genomic profiles and evolutionary history of two lineages of the mycoheterotrophic orchid Gastrodia, both of which independently acquired complete cleistogamy, based on detailed sampling and a combination of simple sequence repeat (SSR), multiplexed ISSR genotyping by sequencing (MIG-seq) and RNA-seq data. Our analysis reveals clear species delimitation, with no evidence of introgression between the completely cleistogamous species and their co-occurring allogamous sisters. Intriguingly, all analyses indicate that both the completely cleistogamous Gastrodia species and their allogamous sisters exhibit genetic profiles typical of self-pollinating plants. This pattern suggests that their ancestors, probably bearing allogamous flowers, had already evolved mechanisms to mitigate the deleterious effects of selfing, potentially facilitating the emergence of complete cleistogamy through benefits such as reproductive assurance, enhanced colonization ability and species reinforcement. Meanwhile, further analyses suggest that complete cleistogamy evolved very recently (possibly within the last 1000-2000 years) in these two Gastrodia lineages. Combined with the scant evidence of complete cleistogamy outside Gastrodia, our findings imply a limited and ephemeral role for complete cleistogamy in plant speciation.

Electronic insights into the role of nuclear quantum effects in proton transfer reactions of nucleobase pairs

Double proton transfer in nucleobase pairs leads to point mutations in nucleic acids. A series of constrained nuclear-electronic orbital calculations combined with natural bond orbital and non-covalent interaction analyses, and kinetic studies have quantitatively revealed the importance of nuclear quantum effects (NQEs) in the reaction. Compared with the classical treatment of the nuclei, the probability of forming the tautomeric isomers of Cytosine-Guanine, when explicitly accounting for NQEs, increased by a factor of 8.0. This outcome can be attributed to enhancing the interaction between the orbitals at the reactive site due to NQEs, which increased the number of electrons occupying the antibonding orbitals.

Ectopic Gene Conversion Causing Quantitative Trait Variation

Why is there so much non-neutral genetic variation segregating in natural populations? We dissect function and evolution of a near-cryptic quantitative trait locus (QTL) for defense metabolites in Arabidopsis using the CRISPR/Cas9 system and nucleotide polymorphism patterns. The QTL is explained by genetic variation in a family of 4 tightly linked indole-glucosinolate O-methyltransferase genes. Some of this variation appears to be maintained by balancing selection, some appears to be generated by non-reciprocal transfer of sequence, also known as ectopic gene conversion (EGC), between functionally diverged gene copies. Here, we elucidate how EGC, as an inevitable consequence of gene duplication, could be a general mechanism for generating genetic variation for fitness traits.

Evolutionary History and Rhizosphere Microbial Community Composition in Domesticated Hops (Humulus lupulus L.)

Humulus lupulus L., commonly known as hop, is a perennial crop grown worldwide and is well known for its pharmacological, commercial, and most importantly brewing applications. For hundreds of years, hop has undergone intense artificial selection, with over 250 cultivated varieties being developed worldwide, all displaying differences in key characteristics such as bitter acid concentrations, flavour and aroma profiles, changes in photoperiod, growth, and pathogen/pest resistances. Previous studies have individually explored differences between cultivars, aiming to identify markers that can quickly and cost-effectively differentiate between cultivars. However, little is known about their evolutionary history and the variability in their associated rhizospheric microbial communities. Coupling phenotypic, genomic, and soil metagenomic data, our study explores the global population structure and domestication history of 98 hop cultivars. We assessed differences in growth rates, rates of viral infection, usage of dissolvable nitrogen, and soil microbial community compositions between US and non-US based cultivars. Our study revealed that worldwide hop cultivars cluster into four subpopulations: Central European, English, and American ancestry as previously reported, and one new group, the Nobles, revealing further substructure amongst Central European cultivars. Modelling the evolutionary history of domesticated hop reveals divergence of the common ancestors of modern US cultivars around 2800 years before present (ybp), and more recent divergences with gene flow across English, Central European, and Noble cultivars, reconciled with key events in human history and migrations. Furthermore, cultivars of US origin were shown to overall outperform non-US cultivars in both growth rates and usage of dissolvable nitrogen and display novel microbial composition under common-garden settings in the United States.

A more elaborate genetic clock for clonal species

The genetic clock is a well-established tool used in evolutionary biology for estimating divergence times between species, individuals, or cells based on DNA sequence changes. Yu et al. have revisited the clock to make it applicable to clonal multicellular organisms that expand through asexual reproduction mechanisms, enabling more comprehensive evolutionary tracking.

Incorporating genetic load contributes to predicting Arabidopsis thaliana's response to climate change

Understanding how species respond to climate change can facilitate species conservation and crop breeding. Current prediction frameworks about population vulnerability focused on predicting range shifts or local adaptation but ignored genetic load, which is also crucial for adaptation. By analyzing 1115 globally distributed Arabidopsis thaliana natural accessions, we find that effective population size (Ne) is the major contributor of genetic load variation, both along genome and among populations, and can explain 74-94% genetic load variation in natural populations. Intriguingly, Ne affects genetic load by changing both effectiveness of purifying selection and GC biased gene conversion strength. In particular, by incorporating genetic load, genetic offset and species distribution models (SDM), we predict that, the populations at species' range edge are generally at higher risk. The populations at the eastern range perform poorer in all aspects, southern range have higher genetic offset and lower SDM suitability, while northern range have higher genetic load. Among the diverse natural populations, the Yangtze River basin population is the most vulnerable population under future climate change. Overall, here we deciphered the driving forces of genetic load in A. thaliana, and incorporated SDM, local adaptation and genetic load to predict the fate of populations under future climate change.

Hill-Robertson interference may bias the inference of fitness effects of new mutations in highly selfing species

The accurate estimation of the distribution of fitness effects (DFE) of new mutations is critical for population genetic inference but remains a challenging task. While various methods have been developed for DFE inference using the site frequency spectrum of putatively neutral and selected sites, their applicability in species with diverse life history traits and complex demographic scenarios is not well understood. Selfing is common among eukaryotic species and can lead to decreased effective recombination rates, increasing the effects of selection at linked sites, including interference between selected alleles. We employ forward simulations to investigate the limitations of current DFE estimation approaches in the presence of selfing and other model violations, such as linkage, departures from semidominance, population structure, and uneven sampling. We find that distortions of the site frequency spectrum due to Hill-Robertson interference in highly selfing populations lead to mis-inference of the deleterious DFE of new mutations. Specifically, when inferring the distribution of selection coefficients, there is an overestimation of nearly neutral and strongly deleterious mutations and an underestimation of mildly deleterious mutations when interference between selected alleles is pervasive. In addition, the presence of cryptic population structure with low rates of migration and uneven sampling across subpopulations leads to the false inference of a deleterious DFE skewed towards effectively neutral/mildly deleterious mutations. Finally, the proportion of adaptive substitutions estimated at high rates of selfing is substantially overestimated. Our observations apply broadly to species and genomic regions with little/no recombination and where interference might be pervasive.

The population genetics of convergent adaptation in maize and teosinte is not locally restricted

What is the genetic architecture of local adaptation and what is the geographic scale over which it operates? We investigated patterns of local and convergent adaptation in five sympatric population pairs of traditionally cultivated maize and its wild relative teosinte (Zea mays subsp. parviglumis). We found that signatures of local adaptation based on the inference of adaptive fixations and selective sweeps are frequently exclusive to individual populations, more so in teosinte compared to maize. However, for both maize and teosinte, selective sweeps are also frequently shared by several populations, and often between subspecies. We were further able to infer that selective sweeps were shared among populations most often via migration, though sharing via standing variation was also common. Our analyses suggest that teosinte has been a continued source of beneficial alleles for maize, even after domestication, and that maize populations have facilitated adaptation in teosinte by moving beneficial alleles across the landscape. Taken together, our results suggest local adaptation in maize and teosinte has an intermediate geographic scale, one that is larger than individual populations but smaller than the species range.

Spontaneous and environment induced genomic alterations in yeast model

While genomic alterations are fundamental to biological evolution, enabling adaptation and diversity, they can also result in detrimental outcomes, such as the development of genetic diseases including cancer. The budding yeast Saccharomyces cerevisiae serves as an exemplary model for investigating the mechanisms behind various genomic alterations, including point mutations, chromosomal rearrangements, and whole-chromosome aneuploidy. In this review, we highlight the application of genetic screening systems to assess the mutagenic effects of physical and chemical agents efficiently. Additionally, we discuss the utilization of high-throughput sequencing technologies to uncover comprehensive genomic alterations and rare genetic events. We provide a detailed summary of the features of genomic alterations and discuss the genetic mechanisms driving these changes under both spontaneous and stress-induced conditions. Given the high conservation of DNA replication and repair machinery across different organisms, the insights gained from studies on yeast offer valuable perspectives for understanding the delicate balance between genome plasticity and integrity in other species.

Decoding resilience: ecology, regulation, and evolution of biosynthetic gene clusters

Secondary metabolism is crucial for plant survival and can generate chemistry with nutritional, therapeutic, and industrial value. Biosynthetic genes of selected secondary metabolites cluster within localised chromosomal regions. The arrangement of these biosynthetic gene clusters (BGCs) challenges the long-held model of random gene order in eukaryotes, raising questions about their regulation, ecological significance, and evolution. In this review, we address these questions by exploring the contribution of BGCs to ecologically relevant plant-biotic interactions, while also evaluating the molecular-(epi)genetic mechanisms controlling their coordinated stress- and tissue-specific expression. Based on evidence that BGCs have distinct chromatin signatures and are enriched with transposable elements (TEs), we integrate emerging hypotheses into an updated evolutionary model emphasising how stress-induced epigenetic processes have shaped BGC formation.

DNA methylation dynamics in gymnosperm duplicate genes: implications for genome evolution and stress adaptation

Duplicate genes are pivotal in driving evolutionary innovation, often exhibiting expression divergence that offers a system to investigate the role of DNA methylation in transcriptional regulation. However, previous studies have predominantly focused on angiosperms, leaving the methylation patterns in major lineages of land plants still unclear. This study explores DNA methylation evolution in duplicate genes across representative gymnosperm species with large genomes, spanning over 300 million years, using genomic, transcriptomic, and high-depth DNA methylomic data. We observed variations in DNA methylation levels along gene bodies, flanking regions, and methylation statuses of coding regions across different duplication types. Biased divergences in DNA methylation and gene expression frequently occurred between duplicate copies. Specifically, methylation divergences in the 2-kb downstream regions negatively correlated with gene expression. Both CG and CHG DNA methylation in gene bodies were positively correlated with gene length, suggesting these methylation types may function as an epigenomic buffer to mitigate the adverse impact of gene length on expression. Duplicate genes exhibiting both methylation and expression divergences were notably enriched in adaptation-related biological processes, suggesting that DNA methylation may aid adaptive evolution in gymnosperms by regulating stress response genes. Changes in expression levels correlated with switches in methylation status within coding regions of transposed duplicates. Specifically, depletion for CG methylation or enrichment for non-CG methylation significantly reduced the expression of translocated copies. This correlation suggests that DNA methylation may reduce genetic redundancy by silencing translocated copies. Our study highlights the significance of DNA methylation in plant genome evolution and stress adaptation.

The gap-free genome assembly and multi-omics analyses illustrate the evolutionary history and the synthesis of medicinal components of Ligustrum lucidum

Ligustrum lucidum, an important ornamental and medicinal plant in the Oleaceae family, has become a globally notorious invasive species because of its high adaptability. However, the lack of a high quality genome has hindered the understanding of the genetic basis for its broad adaptability and the mechanisms of its medicinal component synthesis. In this study, we successfully assembled a gap free telomere to telomere genome of L. lucidum. Through the reconstruction of the phylogenetic tree, we presented evidence that Jasmineae diverged from Oleaceae earlier than Forsythieae, which challenges the existing phylogenetic relationships within Oleaceae. Comparative genomics has found that two rounds of gene family expansions in L. lucidum significantly increased the number of genes related to its biotic and abiotic stress resistances, which may lay the genetic foundation for its broad adaptability. Among the L. lucidum fruits from January to March in our investigation, the results of metabolome show that March may be the optimal harvesting time. By integrating genomic, transcriptomic, and metabolomic data, we identified candidate genes involved in the synthesis of compounds, such as specnuezhenide, salidroside, and oleanolic acid. This study provides valuable genomic resources for comparative genomics studies within Oleaceae and for the genetic research of L. lucidum.

Ultra-deep sequencing of somatic mutations induced by a maize transposon

Cells accumulate mutations throughout development, contributing to cancer, aging, and evolution. Quantitative data on the abundance of de novo mutations within plants or animals are limited, as new mutations are often rare within a tissue and fall below the limits of current sequencing depths and error rates. Here, we show that mutations induced by the maize Mutator (Mu) transposon can be reliably quantified down to a detection limit of 1 part in 12,000. We measured the abundance of millions of de novo Mu insertions across four tissue types. Within a tissue, the distribution of de novo Mu allele frequencies was highly reproducible between plants, showing that, despite the stochastic nature of mutation, repeated statistical patterns of mutation abundance emerge. In contrast, there were significant differences in the allele frequency distribution between tissues. At the extremes, root was dominated by a small number of highly abundant de novo insertions, while endosperm was characterized by thousands of insertions at low allele frequencies. Finally, we used the measured pollen allele frequencies to reinterpret a classic genetic experiment, showing that evidence for late Mu activity in pollen are better explained by cell division statistics. These results provide insight into the complexity of mutation accumulation in multicellular organisms and a system to interrogate the factors that shape mutation abundance.

TIPPo: A User-Friendly Tool for De Novo Assembly of Organellar Genomes with High-Fidelity Data

Plant cells have two major organelles with their own genomes: chloroplasts and mitochondria. While chloroplast genomes tend to be structurally conserved, the mitochondrial genomes of plants, which are much larger than those of animals, are characterized by complex structural variation. We introduce TIPPo, a user-friendly, reference-free assembly tool that uses PacBio high-fidelity long-read data and that does not rely on genomes from related species or nuclear genome information for the assembly of organellar genomes. TIPPo employs a deep learning model for initial read classification and leverages k-mer counting for further refinement, significantly reducing the impact of nuclear insertions of organellar DNA on the assembly process. We used TIPPo to completely assemble a set of 54 complete chloroplast genomes. No other tool was able to completely assemble this set. TIPPo is comparable with PMAT in assembling mitochondrial genomes from most species but does achieve even higher completeness for several species. We also used the assembled organelle genomes to identify instances of nuclear plastid DNA (NUPTs) and nuclear mitochondrial DNA (NUMTs) insertions. The cumulative length of NUPTs/NUMTs positively correlates with the size of the nuclear genome, suggesting that insertions occur stochastically. NUPTs/NUMTs show predominantly C:G to T:A changes, with the mutated cytosines typically found in CG and CHG contexts, suggesting that degradation of NUPT and NUMT sequences is driven by the known elevated mutation rate of methylated cytosines. Small interfering RNA loci are enriched in NUPTs and NUMTs, consistent with the RdDM pathway mediating DNA methylation in these sequences.

Genomic implications of the repeated shift to self-fertilization across a species' geographic distribution

The ability to self-fertilize often varies among closely related hermaphroditic plant species, though, variation can also exist within species. In the North American Arabidopsis lyrata, the shift from self-incompatibility (SI) to selfing established in multiple regions independently, mostly since recent postglacial range expansion. This has made the species an ideal model for the investigation of the genomic basis of the breakdown of SI and its population genetic consequences. By comparing nearby selfing and outcrossing populations across the entire species' geographic distribution, we investigated variation at the self-incompatibility (S-)locus and across the genome. Furthermore, a diallel crossing experiment on one mixed-mating population was performed to gain insight into the inheritance of mating system variation. We confirmed that the breakdown of SI had evolved in several S-locus backgrounds. The diallel suggested the involvement of biparental contributions with dominance relations. Though, the population-level genome-wide association study did not single out clear-cut candidate genes but several regions with one near the S-locus. On the implication side, selfing as compared to outcrossing populations had less than half of the genomic diversity, while the number and length of runs of homozygosity (ROHs) scaled with the degree of inbreeding. Selfing populations with a history of long expansion had the longest ROHs. The results highlight that mating system shift to selfing, its genetic underpinning and the likely negative genomic consequences for evolutionary potential can be strongly interlinked with past range dynamics.

Spontaneous and salt stress-induced molecular instability in the progeny of MSH7 deficient Arabidopsis thaliana plants

The MSH7 protein is a binding partner of MSH2 forming the MutSγ complex. This complex contributes to the plant mismatch repair (MMR) system by recognizing DNA base-base mismatches. Here, we evaluated the impact of MSH7 on genetic diversity of the tenth generation (G10) of wild type and MSH7 deficient Arabidopsis thaliana plants before and after two days exposure to 100 mM NaCl. Genetic diversity was assessed using inter simple sequence repeats (ISSR) and high-resolution melting (HRM) analyses. ISSR analyses revealed a 6.7 % or 5.8 % average polymorphism in the G10 of wild type before and after a short-term salt stress, respectively, and a 64.4 % or 72.1 % average polymorphism in the G10 of msh7 mutant plants before and after salt treatment, respectively. Interestingly, several ISSR markers showed different polymorphism patterns after salt stress compared with the control before treatment. We next compared the percentage of the G10 of wild type and msh7 seedlings with polymorphic bands. Statistically significant differences between genotypes but not due to the salt treatment were observed. In addition, co-amplification at lower temperature-PCR followed by HRM analysis was performed. Of the five assayed HRM loci, two loci allowed the discrimination of fragment alleles between genotypes and two loci, between conditions. We conclude that MSH7 deficient A. thaliana mutants accumulated mutations over 10 generations, and that two days of salt stress caused a further increase in new mutations, thus enhancing genetic diversity that may favor new traits associated with stress tolerance, fitness, and adaptation.

Epigenetic gene regulation in plants and its potential applications in crop improvement

DNA methylation, also known as 5-methylcytosine, is an epigenetic modification that has crucial functions in plant growth, development and adaptation. The cellular DNA methylation level is tightly regulated by the combined action of DNA methyltransferases and demethylases. Protein complexes involved in the targeting and interpretation of DNA methylation have been identified, revealing intriguing roles of methyl-DNA binding proteins and molecular chaperones. Structural studies and in vitro reconstituted enzymatic systems have provided mechanistic insights into RNA-directed DNA methylation, the main pathway catalysing de novo methylation in plants. A better understanding of the regulatory mechanisms will enable locus-specific manipulation of the DNA methylation status. CRISPR-dCas9-based epigenome editing tools are being developed for this goal. Given that DNA methylation patterns can be stably transmitted through meiosis, and that large phenotypic variations can be contributed by epimutations, epigenome editing holds great promise in crop breeding by creating additional phenotypic variability on the same genetic material.

Allometric scaling of somatic mutation and epimutation rates in trees

How long-lived trees escape "mutational meltdown" despite centuries of continuous growth remains puzzling. Here we integrate recent studies to show that the yearly rate of somatic mutations and epimutations (μY) scales inversely with generation time (G), and follows the same allometric power law found in mammals (μY ∝ G-1). Deeper insights into the scaling function may permit predictions of somatic (epi)mutation rates from life-history traits without the need for genomic data.

Retrotransposon-driven environmental regulation of FLC leads to adaptive response to herbicide

The mobilization of transposable elements is a potent source of mutations. In plants, several stransposable elements respond to external cues, fuelling the hypothesis that natural transposition can create environmentally sensitive alleles for adaptation. Here we report on the detailed characterization of a retrotransposon insertion within the first intron of the Arabidopsis floral-repressor gene FLOWERING LOCUS C (FLC) and the discovery of its role for adaptation. The insertion mutation augments the environmental sensitivity of FLC by affecting the balance between coding and non-coding transcripts in response to stress, thus expediting flowering. This balance is modulated by DNA methylation and orchestrated by IBM2, a factor involved in the processing of intronic heterochromatic sequences. The stress-sensitive allele of FLC has spread across populations subjected to recurrent chemical weeding, and we show that retrotransposon-driven acceleration of the life cycle represents a rapid response to herbicide application. Our work provides a compelling example of a transposable element-driven environmentally sensitive allele that confers an adaptive response in nature.

Climate-related naturally occurring epimutation and their roles in plant adaptation in A. thaliana

DNA methylation has been proposed to be an important mechanism that allows plants to respond to their environments sometimes entirely uncoupled from genetic variation. To understand the genetic basis, biological functions and climatic relationships of DNA methylation at a population scale in Arabidopsis thaliana, we performed a genome-wide association analysis with high-quality single nucleotide polymorphisms (SNPs), and found that ~56% on average, especially in the CHH sequence context (71%), of the differentially methylated regions (DMRs) are not tagged by SNPs. Among them, a total of 3235 DMRs are significantly associated with gene expressions and potentially heritable. 655 of the 3235 DMRs are associated with climatic variables, and we experimentally verified one of them, HEI10 (HUMAN ENHANCER OF CELL INVASION NO.10). Such epigenetic loci could be subjected to natural selection thereby affecting plant adaptation, and would be expected to be an indicator of accessions at risk. We therefore incorporated these climate-related DMRs into a gradient forest model, and found that the natural A. thaliana accessions in Southern Europe that may be most at risk under future climate change. Our findings highlight the importance of integrating DNA methylation that is independent of genetic variations, and climatic data to predict plants' vulnerability to future climate change.

Global characterization of somatic mutations and DNA methylation changes during vegetative propagation in strawberries

Somatic mutations arise and accumulate during tissue culture and vegetative propagation, potentially affecting various traits in horticultural crops, but their characteristics are still unclear. Here, somatic mutations in regenerated woodland strawberry derived from tissue culture of shoot tips under different conditions and 12 cultivated strawberry individuals are analyzed by whole genome sequencing. The mutation frequency of single nucleotide variants is significantly increased with increased hormone levels or prolonged culture time in the range of 3.3 × 10-8-3.0 × 10-6 mutations per site. CG methylation shows a stable reduction (0.71%-8.03%) in regenerated plants, and hypoCG-DMRs are more heritable after sexual reproduction. A high-quality haplotype-resolved genome is assembled for the strawberry cultivar "Beni hoppe." The 12 "Beni hoppe" individuals randomly selected from different locations show 4731-6005 mutations relative to the reference genome, and the mutation frequency varies among the subgenomes. Our study has systematically characterized the genetic and epigenetic variants in regenerated woodland strawberry plants and different individuals of the same strawberry cultivar, providing an accurate assessment of somatic mutations at the genomic scale and nucleotide resolution in plants.

Somatic mutation rates scale with time not growth rate in long-lived tropical trees

The rates of appearance of new mutations play a central role in evolution. However, mutational processes in natural environments and their relationship with growth rates are largely unknown, particular in tropical ecosystems with high biodiversity. Here, we examined the somatic mutation landscapes of two tropical trees, Shorea laevis (slow-growing) and S. leprosula (fast-growing), in central Borneo, Indonesia. Using newly constructed genomes, we identified a greater number of somatic mutations in tropical trees than in temperate trees. In both species, we observed a linear increase in the number of somatic mutations with physical distance between branches. However, we found that the rate of somatic mutation accumulation per meter of growth was 3.7-fold higher in S. laevis than in S. leprosula. This difference in the somatic mutation rate was scaled with the slower growth rate of S. laevis compared to S. leprosula, resulting in a constant somatic mutation rate per year between the two species. We also found that somatic mutations are neutral within an individual, but those mutations transmitted to the next generation are subject to purifying selection. These findings suggest that somatic mutations accumulate with absolute time and older trees have a greater contribution towards generating genetic variation.

Whole-genome assembly and annotation of the firecracker penstemon (Penstemon eatonii)

The penstemons are ornamental annual flowering plants native to the Intermountain West and Rocky Mountains and commonly used for urban landscaping. Elite commercial penstemons are generally susceptible to abiotic stresses, including drought, root rot, cold, and high salinity. Firecracker penstemon (Penstemon eatonii), however, is much more tolerant to these stresses than most elite cultivars. Importantly, firecracker penstemon has been reported to hybridize with many other penstemons and therefore provides the opportunity to develop more tolerant elite cultivars through strategic crossing. To facilitate the study and utilization of firecracker penstemon, we sequenced and annotated the genome of a P. eatonii accession collected from Utah, USA. We also performed low-coverage, whole-genome sequencing of 26 additional accessions from three different varieties of P. eatonii. This chromosome-scale genome assembly is the most contiguous and complete Penstemon genome sequenced to date.

Resilience of genetic diversity in forest trees over the Quaternary

The effect of past environmental changes on the demography and genetic diversity of natural populations remains a contentious issue and has rarely been investigated across multiple, phylogenetically distant species. Here, we perform comparative population genomic analyses and demographic inferences for seven widely distributed and ecologically contrasting European forest tree species based on concerted sampling of 164 populations across their natural ranges. For all seven species, the effective population size, Ne, increased or remained stable over many glacial cycles and up to 15 million years in the most extreme cases. Surprisingly, the drastic environmental changes associated with the Pleistocene glacial cycles have had little impact on the level of genetic diversity of dominant forest tree species, despite major shifts in their geographic ranges. Based on their trajectories of Ne over time, the seven tree species can be divided into three major groups, highlighting the importance of life history and range size in determining synchronous variation in genetic diversity over time. Altogether, our results indicate that forest trees have been able to retain their evolutionary potential over very long periods of time despite strong environmental changes.

Exploring the Relationship Between Gene Expression and Low-Frequency Somatic Mutations in Arabidopsis with Duplex Sequencing

Intragenomic mutation rates can vary dramatically due to transcription-associated mutagenesis or transcription-coupled repair, which vary based on local epigenomic modifications that are nonuniformly distributed across genomes. One feature associated with decreased mutation is higher expression level, which depends on environmental cues. To understand the magnitude of expression-dependent mutation rate variation, we perturbed expression through a heat treatment in Arabidopsis thaliana. We quantified gene expression to identify differentially expressed genes, which we then targeted for mutation detection using duplex sequencing. This approach provided a highly accurate measurement of the frequency of rare somatic mutations in vegetative plant tissues, which has been a recent source of uncertainty. Somatic mutations in plants may be useful for understanding drivers of DNA damage and repair in the germline since plants experience late germline segregation and both somatic and germline cells share common repair machinery. We included mutant lines lacking mismatch repair (MMR) and base excision repair (BER) capabilities to understand how repair mechanisms may drive biased mutation accumulation. We found wild-type (WT) and BER mutant mutation frequencies to be very low (mean variant frequency 1.8 × 10-8 and 2.6 × 10-8, respectively), while MMR mutant frequencies were significantly elevated (1.13 × 10-6). Interestingly, in the MMR mutant lines, there was no difference in the somatic mutation frequencies between temperature treatments or between highly versus lowly expressed genes. The extremely low somatic variant frequencies in WT plants indicate that larger datasets will be needed to address fundamental evolutionary questions about whether environmental change leads to gene-specific changes in mutation rate.

The identification and analysis of meristematic mutations within the apple tree that developed the RubyMac sport mutation

BACKGROUND

Understanding the molecular basis of sport mutations in fruit trees has the potential to accelerate generation of improved cultivars.

RESULTS

For this, we analyzed the genome of the apple tree that developed the RubyMac phenotype through a sport mutation that led to the characteristic fruit coloring of this variety. Overall, we found 46 somatic mutations that distinguished the mutant and wild-type branches of the tree. In addition, we found 54 somatic gene conversions (i.e., loss-of-heterozygosity mutations) that also distinguished the two parts of the tree. Approximately 20% of the mutations were specific to individual cell lineages, suggesting that they originated from the corresponding meristematic layers. Interestingly, the de novo mutations were enriched for GC =  > AT transitions while the gene conversions showed the opposite bias for AT =  > GC transitions, suggesting that GC-biased gene conversions have the potential to counteract the AT-bias of de novo mutations. By comparing the gene expression patterns in fruit skins from mutant and wild-type branches, we found 56 differentially expressed genes including 18 involved in anthocyanin biosynthesis. While none of the differently expressed genes harbored a somatic mutation, we found that some of them in regions of the genome that were recently associated with natural variation in fruit coloration.

CONCLUSION

Our analysis revealed insights in the characteristics of somatic change, which not only included de novo mutations but also gene conversions. Some of these somatic changes displayed strong candidate mutations for the change in fruit coloration in RubyMac.

A review on point mutations via proton transfer in DNA base pairs in the absence and presence of electric fields

This comprehensive review focuses on spontaneous mutations that may occur during DNA replication, the fundamental process responsible for transferring genetic information. In 1963, Löwdin postulated that these mutations are primarily a result of proton transfer reactions within the hydrogen-bonded DNA base pairs. The single and double proton transfer reactions within the base pairs in DNA result in zwitterions and rare tautomers, respectively. For persistent mutations, these products must be generated at high rates and should be thermodynamically stable. This review covers the proton transfer reactions studied experimentally and computationally. The review also examines the influence of externally applied electric fields on the thermodynamics and kinetics of proton transfer reactions within DNA base pairs, and their biological implications.

Microhomology-Mediated End-Joining Chronicles: Tracing the Evolutionary Footprints of Genome Protection

The fidelity of genetic information is essential for cellular function and viability. DNA double-strand breaks (DSBs) pose a significant threat to genome integrity, necessitating efficient repair mechanisms. While the predominant repair strategies are usually accurate, paradoxically, error-prone pathways also exist. This review explores recent advances and our understanding of microhomology-mediated end joining (MMEJ), an intrinsically mutagenic DSB repair pathway conserved across organisms. Central to MMEJ is the activity of DNA polymerase theta (Polθ), a specialized polymerase that fuels MMEJ mutagenicity. We examine the molecular intricacies underlying MMEJ activity and discuss its function during mitosis, where the activity of Polθ emerges as a last-ditch effort to resolve persistent DSBs, especially when homologous recombination is compromised. We explore the promising therapeutic applications of targeting Polθ in cancer treatment and genome editing. Lastly, we discuss the evolutionary consequences of MMEJ, highlighting its delicate balance between protecting genome integrity and driving genomic diversity.

Local cryptic diversity in salinity adaptation mechanisms in the wild outcrossing Brassica fruticulosa

It is normally supposed that populations of the same species should evolve shared mechanisms of adaptation to common stressors due to evolutionary constraint. Here, we describe a system of within-species local adaptation to coastal habitats, Brassica fruticulosa, and detail surprising strategic variability in adaptive responses to high salinity. These different adaptive responses in neighboring populations are evidenced by transcriptomes, diverse physiological outputs, and distinct genomic selective landscapes. In response to high salinity Northern Catalonian populations restrict root-to-shoot Na+ transport, favoring K+ uptake. Contrastingly, Central Catalonian populations accumulate Na+ in leaves and compensate for the osmotic imbalance with compatible solutes such as proline. Despite contrasting responses, both metapopulations were salinity tolerant relative to all inland accessions. To characterize the genomic basis of these divergent adaptive strategies in an otherwise non-saline-tolerant species, we generate a long-read-based genome and population sequencing of 18 populations (nine inland, nine coastal) across the B. fruticulosa species range. Results of genomic and transcriptomic approaches support the physiological observations of distinct underlying mechanisms of adaptation to high salinity and reveal potential genetic targets of these two very recently evolved salinity adaptations. We therefore provide a model of within-species salinity adaptation and reveal cryptic variation in neighboring plant populations in the mechanisms of adaptation to an important natural stressor highly relevant to agriculture.

Improved inference of population histories by integrating genomic and epigenomic data

With the availability of high-quality full genome polymorphism (SNPs) data, it becomes feasible to study the past demographic and selective history of populations in exquisite detail. However, such inferences still suffer from a lack of statistical resolution for recent, for example bottlenecks, events, and/or for populations with small nucleotide diversity. Additional heritable (epi)genetic markers, such as indels, transposable elements, microsatellites, or cytosine methylation, may provide further, yet untapped, information on the recent past population history. We extend the Sequential Markovian Coalescent (SMC) framework to jointly use SNPs and other hyper-mutable markers. We are able to (1) improve the accuracy of demographic inference in recent times, (2) uncover past demographic events hidden to SNP-based inference methods, and (3) infer the hyper-mutable marker mutation rates under a finite site model. As a proof of principle, we focus on demographic inference in Arabidopsis thaliana using DNA methylation diversity data from 10 European natural accessions. We demonstrate that segregating single methylated polymorphisms (SMPs) satisfy the modeling assumptions of the SMC framework, while differentially methylated regions (DMRs) are not suitable as their length exceeds that of the genomic distance between two recombination events. Combining SNPs and SMPs while accounting for site- and region-level epimutation processes, we provide new estimates of the glacial age bottleneck and post-glacial population expansion of the European A. thaliana population. Our SMC framework readily accounts for a wide range of heritable genomic markers, thus paving the way for next-generation inference of evolutionary history by combining information from several genetic and epigenetic markers.

Evolutionary relationship of moso bamboo forms and a multihormone regulatory cascade involving culm shape variation

Moso bamboo (Phyllostachys edulis) known as Mao Zhu (MZ) in Chinese exhibits various forms with distinct morphological characteristics. However, the evolutionary relationship among MZ forms and the mechanisms of culm shape variation are still lacking. Here, the main differences among MZ forms were identified as culm shape variation, which were confirmed by analysing MZ forms (799 bamboo culms) and MZ (458 bamboo culms) populations. To unravel the genetic basis underlying the morphological variations, 20 MZ forms were subjected to whole-genome resequencing. Further analysis yielded 3 230 107 high-quality SNPs and uncovered low genetic diversity and high genotype heterozygosity associated with MZ forms' formation. By integrating the SNP data of 427 MZ individuals representing 15 geographic regions, the origins of eight MZ forms were successfully traced using the phylogenetic tree and the identified common heterozygous loci. Meanwhile, transcriptomic analysis was performed using shoots from MZ and its two forms with culm shape variation. The results, combined with genomic analyses, demonstrated that hormone signalling related genes played crucial roles in culm variation. Co-expression network analysis uncovered genes associated with multiple plant hormone signal transduction, especially auxin and cytokinin were involved in culm shape variation. Furthermore, the regulatory relationships of a specific transcription factor and their target genes associated with auxin and ethylene signalling were validated by yeast one-hybrid, electrophoretic mobility shift assays, and dual-luciferase reporter. Overall, this study provides important insights into the culm shape variation formation in bamboo, which facilitates to breed new varieties with novel culms.

Furfural tolerance of mutant Saccharomyces cerevisiae selected via ionizing radiation combined with adaptive laboratory evolution

BACKGROUND

Lignocellulose is a renewable and sustainable resource used to produce second-generation biofuel ethanol to cope with the resource and energy crisis. Furfural is the most toxic inhibitor of Saccharomyces cerevisiae cells produced during lignocellulose treatment, and can reduce the ability of S. cerevisiae to utilize lignocellulose, resulting in low bioethanol yield. In this study, multiple rounds of progressive ionizing radiation was combined with adaptive laboratory evolution to improve the furfural tolerance of S. cerevisiae and increase the yield of ethanol.

RESULTS

In this study, the strategy of multiple rounds of progressive X-ray radiation combined with adaptive laboratory evolution significantly improved the furfural tolerance of brewing yeast. After four rounds of experiments, four mutant strains resistant to high concentrations of furfural were obtained (SCF-R1, SCF-R2, SCF-R3, and SCF-R4), with furfural tolerance concentrations of 4.0, 4.2, 4.4, and 4.5 g/L, respectively. Among them, the mutant strain SCF-R4 obtained in the fourth round of radiation had a cellular malondialdehyde content of 49.11 nmol/mg after 3 h of furfural stress, a weakening trend in mitochondrial membrane potential collapse, a decrease in accumulated reactive oxygen species, and a cell death rate of 12.60%, showing better cell membrane integrity, stable mitochondrial function, and an improved ability to limit reactive oxygen species production compared to the other mutant strains and the wild-type strain. In a fermentation medium containing 3.5 g/L furfural, the growth lag phase of the SCF-R4 mutant strain was shortened, and its growth ability significantly improved. After 96 h of fermentation, the ethanol production of the mutant strain SCF-R4 was 1.86 times that of the wild-type, indicating that with an increase in the number of irradiation rounds, the furfural tolerance of the mutant strain SCF-R4 was effectively enhanced. In addition, through genome-transcriptome analysis, potential sites related to furfural detoxification were identified, including GAL7, MAE1, PDC6, HXT1, AUS1, and TPK3.

CONCLUSIONS

These results indicate that multiple rounds of progressive X-ray radiation combined with adaptive laboratory evolution is an effective mutagenic strategy for obtaining furfural-tolerant mutants and that it has the potential to tap genes related to the furfural detoxification mechanism.

The vast majority of somatic mutations in plants are layer-specific

BACKGROUND

Plant meristems are structured organs consisting of distinct layers of stem cells, which differentiate into new plant tissue. Mutations in meristematic layers can propagate into large sectors of the plant. However, the characteristics of meristematic mutations remain unclear, limiting our understanding of the genetic basis of somaclonal phenotypic variation.

RESULTS

Here, we analyse the frequency and distribution of somatic mutations in an apricot tree. We separately sequence the epidermis (developing from meristem layer 1) and the flesh (developing from meristem layer 2) of several fruits sampled across the entire tree. We find that most somatic mutations (> 90%) are specific to individual layers. Interestingly, layer 1 shows a higher mutation load than layer 2, implying different mutational dynamics between the layers. The distribution of somatic mutations follows the branching of the tree. This suggests that somatic mutations are propagated to developing branches through axillary meristems. In turn, this leads us to the unexpected observation that the genomes of layer 1 of distant branches are more similar to each other than to the genomes of layer 2 of the same branches. Finally, using single-cell RNA sequencing, we demonstrate that layer-specific mutations were only transcribed in the cells of the respective layers and can form the genetic basis of somaclonal phenotypic variation.

CONCLUSIONS

Here, we analyse the frequency and distribution of somatic mutations with meristematic origin. Our observations on the layer specificity of somatic mutations outline how they are distributed, how they propagate, and how they can impact clonally propagated crops.

Phylogeography and genetic structure of Papaver bracteatum populations in Iran based on genotyping-by-sequencing (GBS)

Papaver bracteatum, known for its high thebaine content and absence of morphine, has emerged as a promising alternative to opium poppy for codeine production. In this study, our objective was to create a diverse panel representing the natural variation of this species in Iran. To achieve this, we employed genotyping-by-sequencing to obtain genome-wide distributed single-nucleotide polymorphisms (SNPs) for phylogeographic analysis, population structure assessment, and evaluation of genetic diversity within P. bracteatum populations. A total of 244 P. bracteatum individuals from 13 distinct populations formed seven genetic groups, along with one highly admixed population. We observed a clear split between the populations inhabiting the Alborz Mts. in the east and Zagros Mts. in the west. In between these mountain ranges, the population of Kachal Mangan exhibited a high degree of genetic admixture between both genetic groups. At or after the end of the last glacial maximum, when climate conditions rapidly changed, all P. bracteatum populations experienced a strong demographic bottleneck reducing the already small effective population sizes further before they increased to their recent strengths. Our results suggest that the ongoing climate change together with human pressure on the species' habitats and limited seed-dispersal ability are potential factors contributing today to rising genetic isolation of P. bracteatum populations. Our results provide genetic data that can be used for conservation measures to safeguard the species' genetic diversity as a resource for future breeding approaches in this medicinally important species.

Heteroplasmy Is Rare in Plant Mitochondria Compared with Plastids despite Similar Mutation Rates

Plant cells harbor two membrane-bound organelles containing their own genetic material-plastids and mitochondria. Although the two organelles coexist and coevolve within the same plant cells, they differ in genome copy number, intracellular organization, and mode of segregation. How these attributes affect the time to fixation or, conversely, loss of neutral alleles is currently unresolved. Here, we show that mitochondria and plastids share the same mutation rate, yet plastid alleles remain in a heteroplasmic state significantly longer compared with mitochondrial alleles. By analyzing genetic variants across populations of the marine flowering plant Zostera marina and simulating organelle allele dynamics, we examine the determinants of allele segregation and allele fixation. Our results suggest that the bottlenecks on the cell population, e.g. during branching or seeding, and stratification of the meristematic tissue are important determinants of mitochondrial allele dynamics. Furthermore, we suggest that the prolonged plastid allele dynamics are due to a yet unknown active plastid partition mechanism. The dissimilarity between plastid and mitochondrial novel allele fixation at different levels of organization may manifest in differences in adaptation processes. Our study uncovers fundamental principles of organelle population genetics that are essential for further investigations of long-term evolution and molecular dating of divergence events.

Polygenic selection to a changing optimum under self-fertilisation

Many traits are polygenic, affected by multiple genetic variants throughout the genome. Selection acting on these traits involves co-ordinated allele-frequency changes at these underlying variants, and this process has been extensively studied in random-mating populations. Yet many species self-fertilise to some degree, which incurs changes to genetic diversity, recombination and genome segregation. These factors cumulatively influence how polygenic selection is realised in nature. Here, we use analytical modelling and stochastic simulations to investigate to what extent self-fertilisation affects polygenic adaptation to a new environment. Our analytical solutions show that while selfing can increase adaptation to an optimum, it incurs linkage disequilibrium that can slow down the initial spread of favoured mutations due to selection interference, and favours the fixation of alleles with opposing trait effects. Simulations show that while selection interference is present, high levels of selfing (at least 90%) aids adaptation to a new optimum, showing a higher long-term fitness. If mutations are pleiotropic then only a few major-effect variants fix along with many neutral hitchhikers, with a transient increase in linkage disequilibrium. These results show potential advantages to self-fertilisation when adapting to a new environment, and how the mating system affects the genetic composition of polygenic selection.

A somatic genetic clock for clonal species

Age and longevity are key parameters for demography and life-history evolution of organisms. In clonal species, a widespread life history among animals, plants, macroalgae and fungi, the sexually produced offspring (genet) grows indeterminately by producing iterative modules, or ramets, and so obscure their age. Here we present a novel molecular clock based on the accumulation of fixed somatic genetic variation that segregates among ramets. Using a stochastic model, we demonstrate that the accumulation of fixed somatic genetic variation will approach linearity after a lag phase, and is determined by the mitotic mutation rate, without direct dependence on asexual generation time. The lag phase decreased with lower stem cell population size, number of founder cells for the formation of new modules, and the ratio of symmetric versus asymmetric cell divisions. We calibrated the somatic genetic clock on cultivated eelgrass Zostera marina genets (4 and 17 years respectively). In a global data set of 20 eelgrass populations, genet ages were up to 1,403 years. The somatic genetic clock is applicable to any multicellular clonal species where the number of founder cells is small, opening novel research avenues to study longevity and, hence, demography and population dynamics of clonal species.

Natural variation in protein kinase D modifies alcohol sensitivity in Caenorhabditis elegans

Differences in naïve alcohol sensitivity between individuals are a strong predictor of later life alcohol use disorders (AUD). However, the genetic bases for alcohol sensitivity (beyond ethanol metabolism) and pharmacological approaches to modulate alcohol sensitivity remain poorly understood. We used a high-throughput behavioral screen to measure acute behavioral sensitivity to alcohol, a model of intoxication, in a genetically diverse set of over 150 wild strains of the nematode Caenorhabditis elegans. We performed a genome-wide association study to identify loci that underlie natural variation in alcohol sensitivity. We identified five quantitative trait loci (QTL) and further show that variants in the C. elegans ortholog of protein kinase D, dkf-2, likely underlie the chromosome V QTL. We found that resistance to intoxication was conferred by dkf-2 loss-of-function mutations as well as partly by a PKD inhibitor in a dkf-2-dependent manner. Protein kinase D might represent a conserved, druggable target to modify alcohol sensitivity with application towards AUD.

Genomics of hybrid parallel origin in Aquilegia ecalcarata

BACKGROUND

The parallel evolution of similar traits or species provides strong evidence for the role of natural selection in evolution. Traits or species that evolved repeatedly can be driven by separate de novo mutations or interspecific gene flow. Although parallel evolution has been reported in many studies, documented cases of parallel evolution caused by gene flow are scarce by comparison. Aquilegia ecalcarata and A. kansuensis belong to the genus of Aquilegia, and are the closest related sister species. Mutiple origins of A. ecalcarata have been reported in previous studies, but whether they have been driven by separate de novo mutations or gene flow remains unclear.

RESULTS

In this study, We conducted genomic analysis from 158 individuals of two repeatedly evolving pairs of A. ecalcarata and A. kansuensis. All samples were divided into two distinct clades with obvious geographical distribution based on phylogeny and population structure. Demographic modeling revealed that the origin of the A. ecalcarata in the Eastern of China was caused by gene flow, and the Eastern A. ecalcarata occurred following introgression from Western A. ecalcarata population. Analysis of Treemix and D-statistic also revealed that a strong signal of gene flow was detected from Western A. ecalcarata to Eastern A. ecalcarata. Genetic divergence and selective sweep analyses inferred parallel regions of genomic divergence and identified many candidate genes associated with ecologically adaptive divergence between species pair. Comparative analysis of parallel diverged regions and gene introgression confirms that gene flow contributed to the parallel evolution of A. ecalcarata.

CONCLUSIONS

Our results further confirmed the multiple origins of A. ecalcarata and highlighted the roles of gene flow. These findings provide new evidence for parallel origin after hybridization as well as insights into the ecological adaptation mechanisms underlying the parallel origins of species.

A historical stepping-stone path for an island-colonizing cactus across a submerged "bridge" archipelago

Here we use population genomic data (ddRAD-Seq) and ecological niche modeling to test biogeographic hypotheses for the divergence of the island-endemic cactus species Cereus insularis Hemsl. (Cereeae; Cactaceae) from its sister species C. fernambucensis Lem. The Cereus insularis grows in the Fernando de Noronha Islands (FNI), a Neotropical archipelago located 350 km off the Brazilian Atlantic Forest (BAF) coast. Phylogeographic reconstructions support a northward expansion by the common ancestor of C. insularis and C. fernambucensis along the mainland BAF coast, with C. insularis diverging from the widespread mainland taxon C. fernambucensis after colonizing FNI in the late Pleistocene. The morphologically distinct C. insularis is monophyletic and nested within C. fernambucensis, as expected from a progenitor-derivative speciation model. We tested alternative biogeographic and demographic hypotheses for the colonization of the FNI using Approximate Bayesian Computation. We found the greatest support for a stepping-stone path that emerged during periods of decreased sea level (the "bridge" hypothesis), in congruence with historical ecological niche modeling that shows highly suitable habitats on stepping-stone islands during glacial periods. The outlier analyses reveal signatures of selection in C. insularis, suggesting a putative role of adaptation driving rapid anagenic differentiation of this species in FNI.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

Chromatin dynamics and RNA metabolism are double-edged swords for the maintenance of plant genome integrity

Maintenance of genome integrity is an essential process in all organisms. Mechanisms avoiding the formation of DNA lesions or mutations are well described in animals because of their relevance to human health and cancer. In plants, they are of growing interest because DNA damage accumulation is increasingly recognized as one of the consequences of stress. Although the cellular response to DNA damage is mostly studied in response to genotoxic treatments, the main source of DNA lesions is cellular activity itself. This can occur through the production of reactive oxygen species as well as DNA processing mechanisms such as DNA replication or transcription and chromatin dynamics. In addition, how lesions are formed and repaired is greatly influenced by chromatin features and dynamics and by DNA and RNA metabolism. Notably, actively transcribed regions or replicating DNA, because they are less condensed and are sites of DNA processing, are more exposed to DNA damage. However, at the same time, a wealth of cellular mechanisms cooperate to favour DNA repair at these genomic loci. These intricate relationships that shape the distribution of mutations along the genome have been studied extensively in animals but much less in plants. In this Review, we summarize how chromatin dynamics influence lesion formation and DNA repair in plants, providing a comprehensive view of current knowledge and highlighting open questions with regard to what is known in other organisms.

Regulation of gene-edited plants in Europe: from the valley of tears into the shining sun?

Some 20 years ago, the EU introduced complex regulatory rules for the growth of transgenic crops, which resulted in a de facto ban to grow these plants in fields within most European countries. With the rise of novel genome editing technologies, it has become possible to improve crops genetically in a directed way without the need for incorporation of foreign genes. Unfortunately, in 2018, the European Court of Justice ruled that such gene-edited plants are to be regulated like transgenic plants. Since then, European scientists and breeders have challenged this decision and requested a revision of this outdated law. Finally, after 5 years, the European Commission has now published a proposal on how, in the future, to regulate crops produced by new breeding technologies. The proposal tries to find a balance between the different interest groups in Europe. On one side, genetically modified plants, which cannot be discerned from their natural counterparts, will exclusively be used for food and feed and are-besides a registration step-not to be regulated at all. On the other side, plants expressing herbicide resistance are to be excluded from this regulation, a concession to the strong environmental associations and NGOs in Europe. Moreover, edited crops are to be excluded from organic farming to protect the business interests of the strong organic sector in Europe. Nevertheless, if this law passes European parliament and council, unchanged, it will present a big step forward toward establishing a more sustainable European agricultural system. Thus, it might soon be possible to develop and grow crops that are more adapted to global warming and whose cultivation will require lower amounts of pesticides. However, there is still a long way to go until the law is passed. Too often, the storm of arguments raised by the opponents, based on irrational fears of mutations and a naive understanding of nature, has fallen on fruitful ground in Europe.

The Carthamus tinctorius L. genome sequence provides insights into synthesis of unsaturated fatty acids

Domesticated safflower (Carthamus tinctorius L.) is a widely cultivated edible oil crop. However, despite its economic importance, the genetic basis underlying key traits such as oil content, resistance to biotic and abiotic stresses, and flowering time remains poorly understood. Here, we present the genome assembly for C. tinctorius variety Jihong01, which was obtained by integrating Oxford Nanopore Technologies (ONT) and BGI-SEQ500 sequencing results. The assembled genome was 1,061.1 Mb, and consisted of 32,379 protein-coding genes, 97.71% of which were functionally annotated. Safflower had a recent whole genome duplication (WGD) event in evolution history and diverged from sunflower approximately 37.3 million years ago. Through comparative genomic analysis at five seed development stages, we unveiled the pivotal roles of fatty acid desaturase 2 (FAD2) and fatty acid desaturase 6 (FAD6) in linoleic acid (LA) biosynthesis. Similarly, the differential gene expression analysis further reinforced the significance of these genes in regulating LA accumulation. Moreover, our investigation of seed fatty acid composition at different seed developmental stages unveiled the crucial roles of FAD2 and FAD6 in LA biosynthesis. These findings offer important insights into enhancing breeding programs for the improvement of quality traits and provide reference resource for further research on the natural properties of safflower.

Recent Advances in Studies of Genomic DNA Methylation and Its Involvement in Regulating Drought Stress Response in Crops

As global arid conditions worsen and groundwater resources diminish, drought stress has emerged as a critical impediment to plant growth and development globally, notably causing declines in crop yields and even the extinction of certain cultivated species. Numerous studies on drought resistance have demonstrated that DNA methylation dynamically interacts with plant responses to drought stress by modulating gene expression and developmental processes. However, the precise mechanisms underlying these interactions remain elusive. This article consolidates the latest research on the role of DNA methylation in plant responses to drought stress across various species, focusing on methods of methylation detection, mechanisms of methylation pattern alteration (including DNA de novo methylation, DNA maintenance methylation, and DNA demethylation), and overall responses to drought conditions. While many studies have observed significant shifts in genome-wide or gene promoter methylation levels in drought-stressed plants, the identification of specific genes and pathways involved remains limited. This review aims to furnish a reference for detailed research into plant responses to drought stress through epigenetic approaches, striving to identify drought resistance genes regulated by DNA methylation, specific signaling pathways, and their molecular mechanisms of action.

Genome-Wide Profile of Mutations Induced by Carbon Ion Beam Irradiation of Dehulled Rice Seeds

As a physical mutagen, carbon ion beam (CIB) irradiation can induce high-frequency mutation, which is user-friendly and environment-friendly in plant breeding. In this study, we resequenced eight mutant lines which were screened out from the progeny of the CIB-irradiated dehulled rice seeds. Among these mutants, CIB induced 135,535 variations, which include single base substitutions (SBSs), and small insertion and deletion (InDels). SBSs are the most abundant mutation, and account for 88% of all variations. Single base conversion is the main type of SBS, and the average ratio of transition and transversion is 1.29, and more than half of the InDels are short-segmented mutation (1-2 bp). A total of 69.2% of the SBSs and InDels induced by CIBs occurred in intergenic regions on the genome. Surprisingly, the average mutation frequency in our study is 9.8 × 10-5/bp and much higher than that of the previous studies, which may result from the relatively high irradiation dosage and the dehulling of seeds for irradiation. By analyzing the mutation of every 1 Mb in the genome of each mutant strain, we found some unusual high-frequency (HF) mutation regions, where SBSs and InDels colocalized. This study revealed the mutation mechanism of dehulled rice seeds by CIB irradiation on the genome level, which will enrich our understanding of the mutation mechanism of CIB radiation and improve mutagenesis efficiency.

Complexity of responses to ionizing radiation in plants, and the impact on interacting biotic factors

In nature, plants are simultaneously exposed to different abiotic (e.g., heat, drought, and salinity) and biotic (e.g., bacteria, fungi, and insects) stresses. Climate change and anthropogenic pressure are expected to intensify the frequency of stress factors. Although plants are well equipped with unique and common defense systems protecting against stressors, they may compromise their growth and development for survival in such challenging environments. Ionizing radiation is a peculiar stress factor capable of causing clustered damage. Radionuclides are both naturally present on the planet and produced by human activities. Natural and artificial radioactivity affects plants on molecular, biochemical, cellular, physiological, populational, and transgenerational levels. Moreover, the fitness of pests, pathogens, and symbionts is concomitantly challenged in radiologically contaminated areas. Plant responses to artificial acute ionizing radiation exposure and laboratory-simulated or field chronic exposure are often discordant. Acute or chronic ionizing radiation exposure may occasionally prime the defense system of plants to better tolerate the biotic stress or could often exhaust their metabolic reserves, making plants more susceptible to pests and pathogens. Currently, these alternatives are only marginally explored. Our review summarizes the available literature on the responses of host plants, biotic factors, and their interaction to ionizing radiation exposure. Such systematic analysis contributes to improved risk assessment in radiologically contaminated areas.

Genetic Causes and Genomic Consequences of Breakdown of Distyly in Linum trigynum

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.