Genome-wide DNA mutations in Arabidopsis plants after multigenerational exposure to high temperatures

Disruption of recombination machinery alters the mutational landscape in plant organellar genomes

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made toward understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar-targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana. We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA, recA1, or recA3. In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild-type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

A Next-Generation Combinatorial Genomic Strategy Scans Genomic Loci Governing Heat Stress Tolerance in Chickpea

In the wake of rising earth temperature, chickpea crop production is haunted by the productivity crisis. Chickpea, a cool season legume manifests tolerance in several agro-physiological level, which is complex quantitative in nature, and regulated by multiple genes and genetic networks. Understanding the molecular genetic basis of this tolerance and identifying key regulators can leverage chickpea breeding against heat stress. This study employed a genomics-assisted breeding strategy utilizing multi-locus GWAS to identify 10 key genomic regions linked to traits contributing to heat stress tolerance in chickpea. These loci subsequently delineated few key candidates and hub regulatory genes, such as RAD23b, CIPK25, AAE19, CK1 and WRKY40, through integrated genomics, transcriptomics and interactive analyses. The differential transcript accumulation of these identified candidates in contrasting chickpea accessions suggests their potential role in heat stress tolerance. Differential ROS accumulation along with their scavengers' transcript abundance aligning with the expression of identified candidates in the contrasting chickpea accessions persuade their regulatory significance. Additionally, their functional significance is ascertained by heterologous expression and subsequent heat stress screening. The high confidence genomic loci and the superior genes and natural alleles delineated here has great potential for swift genomic interventions to enhance heat resilience and yield stability in chickpea.

Genetic and Epigenetic Changes in Arabidopsis thaliana Exposed to Ultraviolet-C Radiation Stress for 25 Generations

Continuous exposure to stress contributes to species diversity and drives microevolutionary processes. It is still unclear, however, whether epigenetic changes, in the form of epimutations such as, for example, differential DNA methylation, are the pre-requisite to speciation events. We hypothesized that continuous stress exposure would increase epigenetic diversity to a higher extent than genetic diversity. In this work, we have analyzed the effect of 25 consecutive generations of UV-C-stress exposure on the Arabidopsis thaliana genome and epigenome. We found no evidence of increased tolerance to UV-C in the progeny of UV-C-stressed plants (F25UV) as compared to the progeny of control plants (F25C). Genetic analysis showed an increased number of single nucleotide polymorphisms (SNPs) and deletions in F25UV plants. Most common SNPs were mutations in cytosines, C to T, C to A, and C to G. Analysis of cytosine methylation showed a significant increase in the percentage of methylated cytosines at CG context in F25UV as compared to F25C or F2C (parental control). The most significant differences between F25UV and either control group were observed in CHG and CHH contexts; the number of hypomethylated cytosines at CHH contexts was over 10 times higher in the F25UC group. F25UV plants clustered separately from other groups in both genomic and epigenomic analyses. GO term analysis of differentially methylated genes revealed enrichments in "DNA or RNA metabolism", "response to stress", "response to biotic and abiotic stimulus", and "signal transduction". Our work thus demonstrates that continuous exposure to UV-C increases genomic and epigenomic diversity in the progeny, with epigenetic changes occurring in many stress-responsive pathways.

Low Mutation Rate and Atypical Mutation Spectrum in Prasinoderma coloniale: Insights From an Early Diverging Green Lineage

Mutations are the ultimate source of genetic diversity on which natural selection and genetic drift act, playing a crucial role in evolution and long-term adaptation. At the molecular level, the spontaneous mutation rate (µ), defined as the number of mutations per base per generation, thus determines the adaptive potential of a species. Through a mutation accumulation experiment, we estimate the mutation rate and spectrum in Prasinoderma coloniale, a phytoplankton species from an early-branching lineage within the Archaeplastida, characterized by an unusually high genomic guanine-cytosine (GC) content (69.8%). We find that P. coloniale has a very low total mutation rate of µ = 2.00 × 10-10. The insertion-deletion mutation rate is almost 5 times lesser than the single nucleotide mutation rate with µID = 3.40 × 10-11 and µSNM = 1.62 × 10-10. Prasinoderma coloniale also exhibits an atypical mutational spectrum: While essentially all other eukaryotes show a bias toward GC to AT mutations, no evidence of this AT-bias is observed in P. coloniale. Since cytosine methylation is known to be mutagenic, we hypothesized that this may result from an absence of C-methylation. Surprisingly, we found high levels of C-methylation (14% in 5mC, 25% in 5mCG contexts). Methylated cytosines did not show increased mutation rates compared with unmethylated ones, not supporting the prevailing notion that C-methylation universally leads to higher mutation rates. Overall, P. coloniale combines a GC-rich genome with a low mutation rate and original mutation spectrum, suggesting the almost universal AT-bias may not have been present in the ancestor of the green lineage.

Adaptive Changes and Genetic Mechanisms in Organisms Under Controlled Conditions: A Review

Adaptive changes encompass physiological, morphological, or behavioral modifications occurring in organisms in response to specific environmental conditions. These modifications may become established within a population through natural selection. While adaptive changes can influence individuals or populations over short timeframes, evolution involves the inheritance and accumulation of these changes over extended periods under environmental pressures through natural selection. At present, addressing climate change, emerging infectious diseases, and food security are the main challenges faced by scientists. A comprehensive and profound understanding of the mechanisms of adaptive evolution is of great significance for solving these problems. The genetic basis of these adaptations can be examined through classical genetics, which includes stochastic gene mutations and chromosomal instability, as well as epigenetics, which involves DNA methylation and histone modifications. These mechanisms not only govern the rate and magnitude of adaptive changes but also affect the transmission of adaptive traits to subsequent generations. In the study of adaptive changes under controlled conditions, short-term controlled experiments are commonly utilized in microbial and animal research to investigate long-term evolutionary trends. However, the application of this approach in plant research remains limited. This review systematically compiles the findings on adaptive changes and their genetic foundations in organisms within controlled environments. It aims to provide valuable insights into fundamental evolutionary processes, offering novel theoretical frameworks and research methodologies for future experimental designs, particularly in the field of plant studies.

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.

Short tandem repeats delineate gene bodies across eukaryotes

Short tandem repeats (STRs) have emerged as important and hypermutable sites where genetic variation correlates with gene expression in plant and animal systems. Recently, it has been shown that a broad range of transcription factors (TFs) are affected by STRs near or in the DNA target binding site. Despite this, the distribution of STR motif repetitiveness in eukaryote genomes is still largely unknown. Here, we identify monomer and dimer STR motif repetitiveness in 5.1 billion 10-bp windows upstream of translation starts and downstream of translation stops in 25 million genes spanning 1270 species across the eukaryotic Tree of Life. We report that all surveyed genomes have gene-proximal shifts in motif repetitiveness. Within genomes, variation in gene-proximal repetitiveness landscapes correlated to the function of genes; genes with housekeeping functions were depleted in upstream and downstream repetitiveness. Furthermore, the repetitiveness landscapes correlated with TF binding sites, indicating that gene function has evolved in conjunction with cis-regulatory STRs and TFs that recognize repetitive sites. These results suggest that the hypermutability inherent to STRs is canalized along the genome sequence and contributes to regulatory and eco-evolutionary dynamics in all eukaryotes.

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.

Genetic and neuro-epigenetic effects of divergent artificial selection for feather pecking behaviour in chickens

Feather pecking (FP) is a repetitive behaviour in chickens, influenced by genetic, epigenetic, and environmental factors, similar to behaviours seen in human developmental disorders (e.g., hyperactivity, autism). This study examines genetic and neuro-epigenetic factors in the thalamus of chickens from lines selected for seven generations for high or low FP behaviour (HFP or LFP). We integrate data on Differentially Methylated Regions (DMRs), Single Nucleotide Polymorphisms (SNPs), and Copy Number Variations (CNVs) in this controlled artificial selection process. Significant differences in behaviour, immunology, and neurology have been reported in these lines. We identified 710 SNPs in these lines that indicate new potentially important genes for FP such as TMPRSS6 (implicated in autism), and SST and ARNT2 (somatostatin function). CNV were the omic level most affected during selection. The largest CNVs found were in RIC3 (gain in HFP) and SH3RF2 (gain in LFP) genes, linked to nicotinic acetylcholine receptor regulation and human oncogenesis, respectively. Our study also suggests that promoters and introns are hotspots for CpG depletion. The overlapping of the omic levels investigated here with data from a public FP Quantitative Trait Loci (QTL) database revealed novel candidate genes for understanding repetitive behaviours, such as RTKN2, associated with Alzheimer's disease in humans. This study suggests CNVs as a crucial initial step for genomic diversification, potentially more impactful than SNPs.

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.

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.

Morphometrics and Phylogenomics of Coca (Erythroxylum spp.) Illuminate Its Reticulate Evolution, With Implications for Taxonomy

South American coca (Erythroxylum coca and E. novogranatense) has been a keystone crop for many Andean and Amazonian communities for at least 8,000 years. However, over the last half-century, global demand for its alkaloid cocaine has driven intensive agriculture of this plant and placed it in the center of armed conflict and deforestation. To monitor the changing landscape of coca plantations, the United Nations Office on Drugs and Crime collects annual data on their areas of cultivation. However, attempts to delineate areas in which different varieties are grown have failed due to limitations around identification. In the absence of flowers, identification relies on leaf morphology, yet the extent to which this is reflected in taxonomy is uncertain. Here, we analyze the consistency of the current naming system of coca and its four closest wild relatives (the "coca clade"), using morphometrics, phylogenomics, molecular clocks, and population genomics. We include name-bearing type specimens of coca's closest wild relatives E. gracilipes and E. cataractarum. Morphometrics of 342 digitized herbarium specimens show that leaf shape and size fail to reliably discriminate between species and varieties. However, the statistical analyses illuminate that rounder and more obovate leaves of certain varieties could be associated with the subtle domestication syndrome of coca. Our phylogenomic data indicate extensive gene flow involving E. gracilipes which, combined with morphometrics, supports E. gracilipes being retained as a single species. Establishing a robust evolutionary-taxonomic framework for the coca clade will facilitate the development of cost-effective genotyping methods to support reliable identification.

Rare but diverse off-target and somatic mutations found in field and greenhouse grown trees expressing CRISPR/Cas9

Introduction: CRISPR gene editing, while highly efficient in creating desired mutations, also has the potential to cause off-target mutations. This risk is especially high in clonally propagated plants, where editing reagents may remain in the genome for long periods of time or in perpetuity. We studied a diverse population of Populus and Eucalyptus trees that had CRISPR/Cas9-containing transgenes that targeted one or two types of floral development genes, homologs of LEAFY and AGAMOUS. Methods: Using a targeted sequence approach, we studied approximately 20,000 genomic sites with degenerate sequence homology of up to five base pairs relative to guide RNA (gRNA) target sites. We analyzed those sites in 96 individual tree samples that represented 37 independent insertion events containing one or multiples of six unique gRNAs. Results: We found low rates of off-target mutations, with rates of 1.2 × 10-9 in poplar and 3.1 × 10-10 in eucalypts, respectively, comparable to that expected due to sexual reproduction. The rates of mutation were highly idiosyncratic among sites and not predicted by sequence similarity to the target sites; a subset of two gRNAs showed off-target editing of four unique genomic sites with up to five mismatches relative to the true target sites, reaching fixation in some gene insertion events and clonal ramets. The location of off-target mutations relative to the PAM site were essentially identical to that seen with on-target CRISPR mutations. Discussion: The low rates observed support many other studies in plants that suggest that the rates of off-target mutagenesis from CRISPR/Cas9 transgenes are negligible; our study extends this conclusion to trees and other long-lived plants where CRISPR/Cas9 transgenes were present in the genome for approximately four years.

Disruption of recombination machinery alters the mutational landscape in plant organellar genomes

Land plant organellar genomes have extremely low rates of point mutation yet also experience high rates of recombination and genome instability. Characterizing the molecular machinery responsible for these patterns is critical for understanding the evolution of these genomes. While much progress has been made towards understanding recombination activity in land plant organellar genomes, the relationship between recombination pathways and point mutation rates remains uncertain. The organellar targeted mutS homolog MSH1 has previously been shown to suppress point mutations as well as non-allelic recombination between short repeats in Arabidopsis thaliana. We therefore implemented high-fidelity Duplex Sequencing to test if other genes that function in recombination and maintenance of genome stability also affect point mutation rates. We found small to moderate increases in the frequency of single nucleotide variants (SNVs) and indels in mitochondrial and/or plastid genomes of A. thaliana mutant lines lacking radA, recA1, or recA3. In contrast, osb2 and why2 mutants did not exhibit an increase in point mutations compared to wild type (WT) controls. In addition, we analyzed the distribution of SNVs in previously generated Duplex Sequencing data from A. thaliana organellar genomes and found unexpected strand asymmetries and large effects of flanking nucleotides on mutation rates in WT plants and msh1 mutants. Finally, using long-read Oxford Nanopore sequencing, we characterized structural variants in organellar genomes of the mutant lines and show that different short repeat sequences become recombinationally active in different mutant backgrounds. Together, these complementary sequencing approaches shed light on how recombination may impact the extraordinarily low point mutation rates in plant organellar genomes.

H3K4me1 recruits DNA repair proteins in plants

DNA repair proteins can be recruited by their histone reader domains to specific epigenomic features, with consequences on intragenomic mutation rate variation. Here, we investigated H3K4me1-associated hypomutation in plants. We first examined 2 proteins which, in plants, contain Tudor histone reader domains: PRECOCIOUS DISSOCIATION OF SISTERS 5 (PDS5C), involved in homology-directed repair, and MUTS HOMOLOG 6 (MSH6), a mismatch repair protein. The MSH6 Tudor domain of Arabidopsis (Arabidopsis thaliana) binds to H3K4me1 as previously demonstrated for PDS5C, which localizes to H3K4me1-rich gene bodies and essential genes. Mutations revealed by ultradeep sequencing of wild-type and msh6 knockout lines in Arabidopsis show that functional MSH6 is critical for the reduced rate of single-base substitution (SBS) mutations in gene bodies and H3K4me1-rich regions. We explored the breadth of these mechanisms among plants by examining a large rice (Oryza sativa) mutation data set. H3K4me1-associated hypomutation is conserved in rice as are the H3K4me1-binding residues of MSH6 and PDS5C Tudor domains. Recruitment of DNA repair proteins by H3K4me1 in plants reveals convergent, but distinct, epigenome-recruited DNA repair mechanisms from those well described in humans. The emergent model of H3K4me1-recruited repair in plants is consistent with evolutionary theory regarding mutation modifier systems and offers mechanistic insight into intragenomic mutation rate variation in plants.

Allelic variation of PmCBF03 contributes to the altitude and temperature adaptability in Japanese apricot (Prunus mume Sieb. et Zucc.)

Japanese apricot is an important subtropical deciduous fruit tree in China, widely distributed in different altitude areas. How does it adapt to the different temperature environments in these areas? In this study, we identified a low-temperature transcription factor PmCBF03 on chromosome 7 through adaptive analysis of populations at different altitudes, which has an early termination single nucleotide polymorphism mutation. There were two different types of variation, PmCBF03A type in high-altitude areas and PmCBF03T type in low-altitude areas. PmCBF03A gene increased the survival rate, Fv/Fm values, antioxidant enzyme activity, and expression levels of antioxidant enzyme genes, and reducing electrolyte leakage and accumulation of reactive oxygen species in transgenic Arabidopsis under low temperature and freezing stress. Simultaneously, PmCBF03A gene promoted the dormancy of transgenic Arabidopsis seeds than wild-type. Biochemical analysis demonstrated that PmCBF03A directly bound to the DRE/CRT element in the promoters of the PmCOR413, PmDAM6 and PmABI5 genes, promoting their transcription and enhanced the cold resistance and dormancy of the overexpressing PmCBF03A lines. While PmCBF03T gene is unable to bind to the promoters of PmDAM6 and PmABI5 genes, leading to early release of dormancy to adapt to the problem of insufficient chilling requirement in low-altitude areas.

Genomic and Epigenomic Changes in the Progeny of Cold-Stressed Arabidopsis thaliana Plants

Plants are continuously exposed to various environmental stresses. Because they can not escape stress, they have to develop mechanisms of remembering stress exposures somatically and passing it to the progeny. We studied the Arabidopsis thaliana ecotype Columbia plants exposed to cold stress for 25 continuous generations. Our study revealed that multigenerational exposure to cold stress resulted in the changes in the genome and epigenome (DNA methylation) across generations. Main changes in the progeny were due to the high frequency of genetic mutations rather than epigenetic changes; the difference was primarily in single nucleotide substitutions and deletions. The progeny of cold-stressed plants exhibited the higher rate of missense non-synonymous mutations as compared to the progeny of control plants. At the same time, epigenetic changes were more common in the CHG (C = cytosine, H = cytosine, adenine or thymine, G = guanine) and CHH contexts and favored hypomethylation. There was an increase in the frequency of C to T (thymine) transitions at the CHH positions in the progeny of cold stressed plants; because this type of mutations is often due to the deamination of the methylated cytosines, it can be hypothesized that environment-induced changes in methylation contribute to mutagenesis and may be to microevolution processes and that RNA-dependent DNA methylation plays a crucial role. Our work supports the existence of heritable stress response in plants and demonstrates that genetic changes prevail.

Investigating low frequency somatic mutations in Arabidopsis with Duplex Sequencing

Mutations are the source of novel genetic diversity but can also lead to disease and maladaptation. The conventional view is that mutations occur randomly with respect to their environment-specific fitness consequences. However, intragenomic mutation rates can vary dramatically due to transcription coupled repair and based on local epigenomic modifications, which are non-uniformly distributed across genomes. One sequence feature associated with decreased mutation is higher expression level, which can vary depending on environmental cues. To understand whether the association between expression level and mutation rate creates a systematic relationship with environment-specific fitness effects, 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 in plant mutation research. 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). These results show that somatic variant frequencies are extremely low in WT plants, indicating that larger datasets will be needed to address the fundamental evolutionary question as to whether environmental change leads to gene-specific changes in mutation rate.

Determination of Differential Alternative Splicing Under Stress Conditions

Alternative splicing (AS) is an important mechanism contributing to stress-induced regulation of gene expression and proteome diversity. Massive sequencing technologies allow the identification of transcripts generated via stress-responsive AS, potentially important for adaptation to stress conditions. Several bioinformatics tools have been developed to identify differentially expressed alternative splicing events/transcripts from RNA-sequencing results. This chapter describes a detailed protocol for differential alternative splicing analysis using the rMATS tool. In addition, we provide guidelines for validation of the detected splice variants by qRT-PCR based on the obtained output files.

Toward Transgene-Free Transposon-Mediated Biological Mutagenesis for Plant Breeding

Genetic diversity is a key factor for plant breeding. The birth of novel genic and genomic variants is also crucial for plant adaptation in nature. Therefore, the genomes of almost all living organisms possess natural mutagenic mechanisms. Transposable elements (TEs) are a major mutagenic force driving genetic diversity in wild plants and modern crops. The relatively rare TE transposition activity during the thousand-year crop domestication process has led to the phenotypic diversity of many cultivated species. The utilization of TE mutagenesis by artificial and transient acceleration of their activity in a controlled mode is an attractive foundation for a novel type of mutagenesis called TE-mediated biological mutagenesis. Here, I focus on TEs as mutagenic sources for plant breeding and discuss existing and emerging transgene-free approaches for TE activation in plants. Furthermore, I also review the non-randomness of TE insertions in a plant genome and the molecular and epigenetic factors involved in shaping TE insertion preferences. Additionally, I discuss the molecular mechanisms that prevent TE transpositions in germline plant cells (e.g., meiocytes, pollen, egg and embryo cells, and shoot apical meristem), thereby reducing the chances of TE insertion inheritance. Knowledge of these mechanisms can expand the TE activation toolbox using novel gene targeting approaches. Finally, the challenges and future perspectives of plant populations with induced novel TE insertions (iTE plant collections) are discussed.

Chasing the mechanisms of ecologically adaptive salinity tolerance

Plants adapted to challenging environments offer fascinating models of evolutionary change. Importantly, they also give information to meet our pressing need to develop resilient, low-input crops. With mounting environmental fluctuation-including temperature, rainfall, and soil salinity and degradation-this is more urgent than ever. Happily, solutions are hiding in plain sight: the adaptive mechanisms from natural adapted populations, once understood, can then be leveraged. Much recent insight has come from the study of salinity, a widespread factor limiting productivity, with estimates of 20% of all cultivated lands affected. This is an expanding problem, given increasing climate volatility, rising sea levels, and poor irrigation practices. We therefore highlight recent benchmark studies of ecologically adaptive salt tolerance in plants, assessing macro- and microevolutionary mechanisms, and the recently recognized role of ploidy and the microbiome on salinity adaptation. We synthesize insight specifically on naturally evolved adaptive salt-tolerance mechanisms, as these works move substantially beyond traditional mutant or knockout studies, to show how evolution can nimbly "tweak" plant physiology to optimize function. We then point to future directions to advance this field that intersect evolutionary biology, abiotic-stress tolerance, breeding, and molecular plant physiology.

Linking temperature dependence of fitness effects of mutations to thermal niche adaptation

Fitness effects of mutations may generally depend on temperature that influences all rate-limiting biophysical and biochemical processes. Earlier studies suggested that high temperatures may increase the availability of beneficial mutations ('more beneficial mutations'), or allow beneficial mutations to show stronger fitness effects ('stronger beneficial mutation effects'). The 'more beneficial mutations' scenario would inevitably be associated with increased proportion of conditionally beneficial mutations at higher temperatures. This in turn predicts that populations in warm environments show faster evolutionary adaptation but suffer fitness loss when faced with cold conditions, and those evolving in cold environments become thermal-niche generalists ('hotter is narrower'). Under the 'stronger beneficial mutation effects' scenario, populations evolving in warm environments would show faster adaptation without fitness costs in cold environments, leading to a 'hotter is (universally) better' pattern in thermal niche adaptation. We tested predictions of the two competing hypotheses using an experimental evolution study in which populations of two model bacterial species, Escherichia coli and Pseudomonas fluorescens, evolved for 2400 generations at three experimental temperatures. Results of reciprocal transplant experiments with our P. fluorescens populations were largely consistent with the 'hotter is narrower' prediction. Results from the E. coli populations clearly suggested stronger beneficial mutation effects at higher assay temperatures, but failed to detect faster adaptation in populations evolving in warmer experimental environments (presumably because of limitation in the supply of genetic variation). Our results suggest that the influence of temperature on mutational effects may provide insight into the patterns of thermal niche adaptation and population diversification across thermal conditions.

Optimizing genome editing efficiency in wheat: Effects of heat treatments and different promoters for single guide RNA expression

Genome editing is a promising method for simultaneously mutagenizing homoeologs in the three subgenomes of wheat (Triticum aestivum L.). However, the mutation rate via genome editing must be improved in order to analyze gene function and to quickly modify agronomic traits in wheat. Here, we examined the Cas9-induced mutation rates in wheat plants using two promoters for single guide RNA (sgRNA) expression and applying heat treatment during Agrobacterium tumefaciens-mediated transformation. Using the TaU6 promoter instead of the OsU6 promoter from rice (Oryza sativa L.) to drive sgRNA expression greatly improved the Cas9-induced mutation rate. Moreover, a heat treatment of 30°C for 1 day during tissue culture increased the Cas9-induced mutation rate and the variety of mutations obtained compared to tissue culture at the normal temperature (25°C). The same heat treatment did not affect the regeneration rates of transgenic plants but tended to increase the number of transgene integration sites in each transgenic plant. These results lay the foundation for improving the Cas9-induced mutation rate in wheat to enhance research on gene function and crop improvement.

Visualization of synthetic retroelement integration reveals determinants of permissivity to retrotransposition

Real-time visualization of retrotransposon mobilization in Arabidopsis reveals that the DNA damage response pathways are involved in the retrotranspositional process.

Environmental gradients reveal stress hubs pre-dating plant terrestrialization

Plant terrestrialization brought forth the land plants (embryophytes). Embryophytes account for most of the biomass on land and evolved from streptophyte algae in a singular event. Recent advances have unravelled the first full genomes of the closest algal relatives of land plants; among the first such species was Mesotaenium endlicherianum. Here we used fine-combed RNA sequencing in tandem with a photophysiological assessment on Mesotaenium exposed to a continuous range of temperature and light cues. Our data establish a grid of 42 different conditions, resulting in 128 transcriptomes and ~1.5 Tbp (~9.9 billion reads) of data to study the combinatory effects of stress response using clustering along gradients. Mesotaenium shares with land plants major hubs in genetic networks underpinning stress response and acclimation. Our data suggest that lipid droplet formation and plastid and cell wall-derived signals have denominated molecular programmes since more than 600 million years of streptophyte evolution-before plants made their first steps on land.

Whole-genome sequence analysis of mutations in rice plants regenerated from zygotes, mature embryos, and immature embryos

Somaclonal variation was studied by whole-genome sequencing in rice plants (Oryza sativa L., 'Nipponbare') regenerated from the zygotes, mature embryos, and immature embryos of a single mother plant. The mother plant and its seed-propagated progeny were also sequenced. A total of 338 variants of the mother plant sequence were detected in the progeny, and mean values ranged from 9.0 of the seed-propagated plants to 37.4 of regenerants from mature embryos. The natural mutation rate of 1.2 × 10-8 calculated using the variants in the seed-propagated plants was consistent with the values reported previously. The ratio of single nucleotide variants (SNVs) among the variants in the seed-propagated plants was 91.1%, which is higher than 56.1% previously reported, and not significantly different from those in the regenerants. Overall, the ratio of transitions to transversions of SNVs was lower in the regenerants as shown previously. Plants regenerated from mature embryos had significantly more variants than different progeny types. Therefore, using zygotes and immature embryos can reduce somaclonal variation during the genetic manipulation of rice.

Persistence of parental age effect on somatic mutation rates across generations in Arabidopsis

In the model plant Arabidopsis thaliana, parental age is known to affect somatic mutation rates in their immediate progeny and here we show that this age dependent effect persists across successive generations. Using a set of detector lines carrying the mutated uidA gene, we examined if a particular parental age maintained across five consecutive generations affected the rates of base substitution (BSR), intrachromosomal recombination (ICR), frameshift mutation (FS), and transposition. The frequency of functional GUS reversions were assessed in seedlings as a function of identical/different parental ages across generations. In the context of a fixed parental age, BSR/ICR rates were unaffected in the first three generations, then dropped significantly in the 4th and increased in most instances in the 5th generation (e.g. BSR (F1 38 = 0.9, F2 38 = 1.14, F3 38 = 1.02, F4 38 = 0.5, F5 38 = 0.76)). On the other hand, with advancing parental ages, BSR/ICR rates remained high in the first two/three generations, with a striking resemblance in the pattern of mutation rates (BSR (F1 38 = 0.9, F1 43 = 0.53, F1 48 = 0.79, F1 53 = 0.83 and F2 38 = 1.14, F2 43 = 0.57, F2 48 = 0.64, F2 53 = 0.94). We adopted a novel approach of identifying and tagging flowers pollinated on a particular day, thereby avoiding biases due to potential emasculation induced stress responses. Our results suggest a time component in counting the number of generations a plant has passed through self-fertilization at a particular age in determining the somatic mutation rates.

Mutation Rate and Spectrum of the Silkworm in Normal and Temperature Stress Conditions

Mutation rate is a crucial parameter in evolutionary genetics. However, the mutation rate of most species as well as the extent to which the environment can alter the genome of multicellular organisms remain poorly understood. Here, we used parents–progeny sequencing to investigate the mutation rate and spectrum of the domestic silkworm (Bombyx mori) among normal and two temperature stress conditions (32 °C and 0 °C). The rate of single-nucleotide mutations in the normal temperature rearing condition was 0.41 × 10−8 (95% confidence interval, 0.33 × 10−8–0.49 × 10−8) per site per generation, which was up to 1.5-fold higher than in four previously studied insects. Moreover, the mutation rates of the silkworm under the stresses are significantly higher than in normal conditions. Furthermore, the mutation rate varies less in gene regions under normal and temperature stresses. Together, these findings expand the known diversity of the mutation rate among eukaryotes but also have implications for evolutionary analysis that assumes a constant mutation rate among species and environments.

Potential and limits of (mal)adaptive mutation rate plasticity in plants

Genetic mutations provide the heritable material for plant adaptation to their environments. At the same time, the environment can affect the mutation rate across plant genomes. However, the extent to which environmental plasticity in mutation rates can facilitate or hinder adaptation remains a longstanding and unresolved question. Emerging discoveries of mechanisms affecting mutation rate variability provide opportunities to consider this question in a new light. Links between chromatin states, transposable elements, and DNA repair suggest cases of adaptive mutation rate plasticity could occur. Yet, numerous evolutionary and biological forces are expected to limit the impact of any such mutation rate plasticity on adaptive evolution. Persistent uncertainty about the significance of mutation rate plasticity on adaptation motivates new experimental and theoretical research relevant to understanding plant responses in changing environments.

Spectrum and Density of Gamma and X-ray Induced Mutations in a Non-Model Rice Cultivar

Physical mutagens are a powerful tool used for genetic research and breeding for over eight decades. Yet, when compared to chemical mutagens, data sets on the effect of different mutagens and dosages on the spectrum and density of induced mutations remain lacking. To address this, we investigated the landscape of mutations induced by gamma and X-ray radiation in the most widely cultivated crop species: rice. A mutant population of a tropical upland rice, Oryza sativa L., was generated and propagated via self-fertilization for seven generations. Five dosages ranging from 75 Gy to 600 Gy in both X-ray and gamma-irradiated material were applied. In the process of a forward genetic screens, 11 unique rice mutant lines showing phenotypic variation were selected for mutation analysis via whole-genome sequencing. Thousands of candidate mutations were recovered in each mutant with single base substitutions being the most common, followed by small indels and structural variants. Higher dosages resulted in a higher accumulation of mutations in gamma-irradiated material, but not in X-ray-treated plants. The in vivo role of all annotated rice genes is yet to be directly investigated. The ability to induce a high density of single nucleotide and structural variants through mutagenesis will likely remain an important approach for functional genomics and breeding.

Transgene Was Silenced in Hybrids between Transgenic Herbicide-Resistant Crops and Their Wild Relatives Utilizing Alien Chromosomes

The commercialization of transgenic herbicide-resistant (HR) crops may cause gene flow risk. If a transgene in progenies of transgenic crops and wild relatives is silencing, these progenies should be killed by the target herbicide, thus, the gene flow risk could be decreased. We obtained the progenies of backcross generations between wild Brassca juncea (AABB, 2n = 36) and glufosinate-resistant transgenic Brassica napus (AACC, 2n = 38, PAT gene located on the C-chromosome). They carried the HR gene but did not express it normally, i.e., gene silencing occurred. Meanwhile, six to nine methylation sites were found on the promoter of PAT in transgene-silencing progenies, while no methylation sites occurred on that in transgene-expressing progenies. In addition, transgene expressing and silencing backcross progenies showed similar fitness with wild Brassica juncea. In conclusion, we elaborate on the occurrence of transgene-silencing event in backcross progenies between transgenic crop utilizing alien chromosomes and their wild relatives, and the DNA methylation of the transgene promoter was an important factor leading to gene silencing. The insertion site of the transgene could be considered a strategy to reduce the ecological risk of transgenic crops, and applied to cultivate lower gene flow HR crops in the future.

Experimentally heat-induced transposition increases drought tolerance in Arabidopsis thaliana

Eukaryotic genomes contain a vast diversity of transposable elements (TEs). Formerly often described as selfish and parasitic DNA sequences, TEs are now recognised as a source of genetic diversity and powerful drivers of evolution. However, because their mobility is tightly controlled by the host, studies experimentally assessing how fast TEs may mediate the emergence of adaptive traits are scarce. We exposed Arabidopsis thaliana high-copy TE lines (hcLines) with up to c. eight-fold increased copy numbers of the heat-responsive ONSEN TE to drought as a straightforward and ecologically highly relevant selection pressure. We provide evidence for increased drought tolerance in five out of the 23 tested hcLines and further pinpoint one of the causative mutations to an exonic insertion of ONSEN in the ribose-5-phosphate-isomerase 2 gene. The resulting loss-of-function mutation caused a decreased rate of photosynthesis, plant size and water consumption. Overall, we show that the heat-induced transposition of a low-copy TE increases phenotypic diversity and leads to the emergence of drought-tolerant individuals in A. thaliana. This is one of the rare empirical examples substantiating the adaptive potential of mobilised stress-responsive TEs in eukaryotes. Our work demonstrates the potential of TE-mediated loss-of-function mutations in stress adaptation.

Challenges and opportunities for plant viruses under a climate change scenario

There is an increasing societal awareness on the enormous threat that climate change may pose for human, animal and plant welfare. Although direct effects due to exposure to heat, drought or elevated greenhouse gasses seem to be progressively more obvious, indirect effects remain debatable. A relevant aspect to be clarified relates to the relationship between altered environmental conditions and pathogen-induced diseases. In the particular case of plant viruses, it is still unclear whether climate change will primarily represent an opportunity for the emergence of new infections in previously uncolonized areas and hosts, or if it will mostly be a strong constrain reducing the impact of plant virus diseases and challenging the pathogen's adaptive capacity. This review focuses on current knowledge on the relationship between climate change and the outcome plant-virus interactions. We summarize work done on how this relationship modulates plant virus pathogenicity, between-host transmission (which include the triple interaction plant-virus-vector), ecology, evolution and management of the epidemics they cause. Considering these studies, we propose avenues for future research on this subject.

Identification of Loci Governing Agronomic Traits and Mutation Hotspots via a GBS-Based Genome-Wide Association Study in a Soybean Mutant Diversity Pool

In this study, we performed a genotyping-by-sequencing analysis and a genome-wide association study of a soybean mutant diversity pool previously constructed by gamma irradiation. A GWAS was conducted to detect significant associations between 37,249 SNPs, 11 agronomic traits, and 6 phytochemical traits. In the merged data set, 66 SNPs on 13 chromosomes were highly associated (FDR p < 0.05) with the following 4 agronomic traits: days of flowering (33 SNPs), flower color (16 SNPs), node number (6 SNPs), and seed coat color (11 SNPs). These results are consistent with the findings of earlier studies on other genetic features (e.g., natural accessions and recombinant inbred lines). Therefore, our observations suggest that the genomic changes in the mutants generated by gamma irradiation occurred at the same loci as the mutations in the natural soybean population. These findings are indicative of the existence of mutation hotspots, or the acceleration of genome evolution in response to high doses of radiation. Moreover, this study demonstrated that the integration of GBS and GWAS to investigate a mutant population derived from gamma irradiation is suitable for dissecting the molecular basis of complex traits in soybeans.

Machine Learning for Plant Stress Modeling: A Perspective towards Hormesis Management

Plant stress is one of the most significant factors affecting plant fitness and, consequently, food production. However, plant stress may also be profitable since it behaves hormetically; at low doses, it stimulates positive traits in crops, such as the synthesis of specialized metabolites and additional stress tolerance. The controlled exposure of crops to low doses of stressors is therefore called hormesis management, and it is a promising method to increase crop productivity and quality. Nevertheless, hormesis management has severe limitations derived from the complexity of plant physiological responses to stress. Many technological advances assist plant stress science in overcoming such limitations, which results in extensive datasets originating from the multiple layers of the plant defensive response. For that reason, artificial intelligence tools, particularly Machine Learning (ML) and Deep Learning (DL), have become crucial for processing and interpreting data to accurately model plant stress responses such as genomic variation, gene and protein expression, and metabolite biosynthesis. In this review, we discuss the most recent ML and DL applications in plant stress science, focusing on their potential for improving the development of hormesis management protocols.

Emerging Strategies Mold Plasticity of Vegetable Plants in Response to High Temperature Stress

As a result of energy consumption and human activities, a large amount of carbon dioxide emissions has led to global warming, which seriously affects the growth and development of plants. Vegetables are an indispensable part of people's diet. In the plant kingdom, a variety of vegetables are highly sensitive to climate change. For them, an increase of just a few degrees above their optimum temperature threshold can result in a loss of yield and quality. Emerging strategies such as practice management and breeding varieties in response to above-optimal temperatures are critical for abiotic stress resistance of vegetable crops. In this study, the function and application of multiple strategies, including breeding improvement, epigenetic modification directed generation of alleles, gene editing techniques, and accumulation of mutations in multigenerational adaptation to abiotic stress, were discussed in vegetable crops. It is believed to be meaningful for plants to build plasticity under high temperature stress, thus generating more genetic structures for heat resistant traits in vegetable products.

Genome Sequences of Two Strains of Prototheca wickerhamii Provide Insight Into the Protothecosis Evolution

The Prototheca alga is the only chlorophyte known to be involved in a series of clinically relevant opportunistic infections in humans and animals, namely, protothecosis. Most pathogenic cases in humans are caused by Prototheca wickerhamii. In order to investigate the evolution of Prototheca and the genetic basis for its pathogenicity, the genomes of two P. wickerhamii strains S1 and S931 were sequenced using Nanopore long-read and Illumina short-read technologies. The mitochondrial, plastid, and nuclear genomes were assembled and annotated including a transcriptomic data set. The assembled nuclear genome size was 17.57 Mb with 19 contigs and 17.45 Mb with 26 contigs for strains S1 and S931, respectively. The number of predicted protein-coding genes was approximately 5,700, and more than 96% of the genes could be annotated with a gene function. A total of 2,798 gene families were shared between the five currently available Prototheca genomes. According to the phylogenetic analysis, the genus of Prototheca was classified in the same clade with A. protothecoides and diverged from Chlorella ~500 million years ago (Mya). A total of 134 expanded genes were enriched in several pathways, mostly in metabolic pathways, followed by biosynthesis of secondary metabolites and RNA transport. Comparative analysis demonstrated more than 96% consistency between the two herein sequenced strains. At present, due to the lack of sufficient understanding of the Prototheca biology and pathogenicity, the diagnosis rate of protothecosis is much lower than the actual infection rate. This study provides an in-depth insight into the genome sequences of two strains of P. wickerhamii isolated from the clinic to contribute to the basic understanding of this alga and explore future prevention and treatment strategies.

Context-Dependent Substitution Dynamics in Plastid DNA Across a Wide Range of Taxonomic Groups

The influence of neighboring base composition, or context, on substitution bias at fourfold degenerate coding sites and in intergenic regions in plastid DNA is compared across the angiosperms, gymnosperms, ferns, liverworts, chlorophytes, stramenopiles and rhodophytes. An influence of flanking base G + C content on the relative rates of transitions and transversions is observed in all lineages and extends up to four nucleotides from the site of substitution in some. Despite finding context effects in all lineages, significant differences were observed between lineages. Overall, the data suggest that context is a general factor affecting mutation bias in plastid DNA but that the dynamics of the influence have evolved over time. It is also shown that, although there are similar effects of context on substitution bias at fourfold degenerate coding sites and at sites within intergenic regions, there are also small but significant differences, suggesting that there could be some selection on some of these sites and that there could be some difference in the mutation and/or repair process between coding and noncoding DNA.

A Simple Heat Treatment Increases SpCas9-Mediated Mutation Efficiency in Arabidopsis

The CRISPR/Cas9 system is now commonly employed for genome editing in various plants such as Arabidopsis, rice and tobacco. In general, in genome editing of the Arabidopsis genome, the SpCas9 and guide RNA genes are introduced into the genome by the floral dip method. Mutations induced in the target sequence by SpCas9 are confirmed after selecting transformants by screening the T1 seed population. The advantage of this method is that genome-edited plants can be isolated easily. However, mutation efficiency in Arabidopsis using SpCas9 is not as high as that achieved in rice and tobacco, which are subjected to a tissue culture step. In this study, we compared four promoters and found that the parsley UBIQITIN promoter is highly active in Arabidopsis meristem tissue. Furthermore, we examined whether a simple heat treatment could improve mutation efficiency in Arabidopsis. Just one heat treatment at 37°C for 24 h increased the mutation efficiency at all four target sites from 3 to 42%, 43 to 62%, 54 to 75% and 89 to 91%, without detectable off-target mutations. We recommend heat treatment of plate-grown plants at 37°C for 24 h as a simple method to increase the efficiency of CRISPR/Cas9-mediated mutagenesis in Arabidopsis.

Pre-mRNA alternative splicing as a modulator for heat stress response in plants

The molecular responses of plants to the important abiotic stress, heat stress (HS), have been extensively studied at the transcriptional level. Alternative splicing (AS) is a post-transcriptional regulatory process in which an intron-containing gene can generate more than one mRNA variant. The impact of HS on the pre-mRNA splicing process has been reported in various eukaryotes but seldom discussed in-depth, especially in plants. Here, we review AS regulation in response to HS in different plant species. We discuss potential molecular mechanisms controlling heat-inducible AS regulation in plants and hypothesize that AS regulation participates in heat-priming establishment and HS memory maintenance. We propose that the pre-mRNA splicing variation is an important regulator of plant HS responses (HSRs).