Transposable elements contribute to activation of maize genes in response to abiotic stress

Expression divergence of BAG gene family in maize under heat stress

Heat stress poses a significant challenge for maize production, especially during the spring when high temperatures disrupt cellular processes, impeding plant growth and development. The B-cell lymphoma-2 (Bcl-2) associated athanogene (BAG) gene family is known to be relatively conserved across various species. It plays a crucial role as molecular chaperone cofactors that are responsible for programmed cell death and tumorigenesis. Once the plant is under heat stress, the BAG genes act as co-chaperones and modulate the molecular functions of HSP70/HSC70 saving the plant from the damage of high temperature stress. The study was planned to identify and characterize the BAG genes for heat stress responsiveness in maize. Twenty-one (21) BAG genes were identified in the latest maize genome. The evolutionary relationship of Zea mays BAGs (ZmBAGs) with Arabidopsis thaliana, Solanum lycopersicum, Theobroma cacao, Sorghum bicolor, Ananas comosus, Physcomitrium patens, Oryza sativa and Populus trichocarpa were represented by the phylogenetic analysis. Differential expressions of BAG gene family in leaf, endosperm, anther, silk, seed and developing embryo depict their contribution to the growth and development. The in-silico gene expression analysis indicated ZmBAG-8 (Zm00001eb170080), and ZmBAG-11 (Zm00001eb237960) showed higher expression under abiotic stresses (cold, heat and salinity). The RT-qPCR further confirmed the expression of ZmBAG-8 and ZmBAG-11 in plant leaf tissue across the contrasting inbred lines and their F1 hybrid (DR-139, UML-1 and DR-139 × UML-1) when exposed to heat stress. Furthermore, the protein-protein interaction networks of ZmBAG-8 and ZmBAG-11 further elucidated their role in stress tolerance related pathways. This research offers a roadmap to plan functional research and utilize ZmBAG genes to enhance heat tolerance in grasses.

Harnessing the power of genomics to develop climate-smart crop varieties: A comprehensive review

Abiotic stresses arising as consequences of climate change pose a serious threat to agricultural productivity on a global scale. Most cultivated crop varieties exhibit susceptibility to such environmental pressures as drought, salinity, and waterlogging. Addressing these abiotic stresses through agronomic means is not only financially burdensome but also often impractical, particularly in the case of abiotic stresses like heat stress. Cultivating resilient varieties that can withstand such pressures emerges as an economically feasible strategy to mitigate these challenges. Nevertheless, the development of stress-tolerant cultivars is hindered by the intricate nature of abiotic stress tolerance, often characterized by low heritability values. Compounding this complexity is the dynamic and multifaceted nature of these stresses, which impede conventional breeding efforts, rendering them painstakingly slow. The identification of molecular markers has emerged as a pivotal advancement in this arena. By pinpointing genomic regions associated with tolerance to abiotic stresses, these markers serve as effective tools for selection and trait introgression. In the post-genomic era, the proliferation of high-density SNP markers has revolutionized breeding strategies. Genomic selection, leveraging these markers, has become the method of choice for addressing polygenic traits with low heritability, such as abiotic stress tolerance. With the functional characterization of many genes being done, precise manipulation through genome editing techniques is gaining significant traction. This review delves into the application of molecular markers in breeding stress-tolerant crop varieties, alongside role of recent genomic techniques in enhancing abiotic stress tolerance. It also explores success stories and identifies potential targets for marker-assisted selection.

Detection and classification of long terminal repeat sequences in plant LTR-retrotransposons and their analysis using explainable machine learning

BACKGROUND

Long terminal repeats (LTRs) represent important parts of LTR retrotransposons and retroviruses found in high copy numbers in a majority of eukaryotic genomes. LTRs contain regulatory sequences essential for the life cycle of the retrotransposon. Previous experimental and sequence studies have provided only limited information about LTR structure and composition, mostly from model systems. To enhance our understanding of these key sequence modules, we focused on the contrasts between LTRs of various retrotransposon families and other genomic regions. Furthermore, this approach can be utilized for the classification and prediction of LTRs.

RESULTS

We used machine learning methods suitable for DNA sequence classification and applied them to a large dataset of plant LTR retrotransposon sequences. We trained three machine learning models using (i) traditional model ensembles (Gradient Boosting), (ii) hybrid convolutional/long and short memory network models, and (iii) a DNA pre-trained transformer-based model using k-mer sequence representation. All three approaches were successful in classifying and isolating LTRs in this data, as well as providing valuable insights into LTR sequence composition. The best classification (expressed as F1 score) achieved for LTR detection was 0.85 using the hybrid network model. The most accurate classification task was superfamily classification (F1=0.89) while the least accurate was family classification (F1=0.74). The trained models were subjected to explainability analysis. Positional analysis identified a mixture of interesting features, many of which had a preferred absolute position within the LTR and/or were biologically relevant, such as a centrally positioned TATA-box regulatory sequence, and TG..CA nucleotide patterns around both LTR edges.

CONCLUSIONS

Our results show that the models used here recognized biologically relevant motifs, such as core promoter elements in the LTR detection task, and a development and stress-related subclass of transcription factor binding sites in the family classification task. Explainability analysis also highlighted the importance of 5'- and 3'- edges in LTR identity and revealed need to analyze more than just dinucleotides at these ends. Our work shows the applicability of machine learning models to regulatory sequence analysis and classification, and demonstrates the important role of the identified motifs in LTR detection.

Structural variants contribute to phenotypic variation in maize

Comprehensively identifying the loci shaping trait variation has been challenging, in part because standard approaches often miss many types of genetic variants. Structural variants (SVs), especially transposable elements (TEs), are likely to affect phenotypic variation but we lack methods that can detect polymorphic structural variants and TEs using short-read sequencing data. Here, we used a whole genome alignment between two maize genotypes to identify polymorphic structural variants and then genotyped a large maize diversity panel for these variants using short-read sequencing data. After characterizing SV variation in the panel, we identified SV polymorphisms that are associated with life history traits and genotype-by-environment (GxE) interactions. While most of the SVs associated with traits contained TEs, only two of the SVs had boundaries that clearly matched TE breakpoints indicative of a TE insertion, while the other polymorphisms were likely caused by deletions. One of the SVs that appeared to be caused by a TE insertion had the most associations with gene expression compared to other trait-associated SVs. All of the SVs associated with traits were in linkage disequilibrium with nearby single nucleotide polymorphisms (SNPs), suggesting that the approach used here did not identify unique associations that would have been missed in a SNP association study. Overall, we have created a technique to genotype SV polymorphisms across a large diversity panel using support from genomic short-read sequencing alignments and connecting this presence/absence SV variation to diverse traits and GxE interactions.

Genetic and epigenetic reprogramming in response to internal and external cues by induced transposon mobilization in Moso bamboo

Long terminal repeat retroelements (LTR-REs) have profound effects on DNA methylation and gene regulation. Despite the vast abundance of LTR-REs in the genome of Moso bamboo (Phyllostachys edulis), an industrial crop in underdeveloped countries, their precise implication of the LTR-RE mobility in stress response and development remains unknown. We investigated the RNA and DNA products of LTR-REs in Moso bamboo under various developmental stages and stressful conditions. Surprisingly, our analyses identified thousands of active LTR-REs, particularly those located near genes involved in stress response and developmental regulation. These genes adjacent to active LTR-REs exhibited an increased expression under stress and are associated with reduced DNA methylation that is likely affected by the induced LTR-REs. Moreover, the analyses of simultaneous mapping of insertions and DNA methylation showed that the LTR-REs effectively alter the epigenetic status of the genomic regions where they inserted, and concomitantly their transcriptional competence which might impact the stress resilience and growth of the host. Our work unveils the unusually strong LTR-RE mobility in Moso bamboo and its close association with (epi)genetic changes, which supports the co-evolution of the parasitic DNAs and host genome in attaining stress tolerance and developmental robustness.

The Genomic SSR Millets Database (GSMDB): enhancing genetic resources for sustainable agriculture

The global population surge demands increased food production and nutrient-rich options to combat rising food insecurity. Climate-resilient crops are vital, with millets emerging as superfoods due to nutritional richness and stress tolerance. Given limited genomic information, a comprehensive genetic resource is crucial to advance millet research. Whole-genome sequencing provides an unprecedented opportunity, and molecular genetic methodologies, particularly simple sequence repeats (SSRs), play a pivotal role in DNA fingerprinting, constructing linkage maps, and conducting population genetic studies. SSRs are composed of repetitive DNA sequences where one to six nucleotides are repeated in tandem and distributed throughout the genome. Different millet species exhibit genomic variations attributed to the presence of SSRs. While SSRs have been identified in a few millet species, the existing information only covers some of the sequenced genomes. Moreover, there is an absence of complete gene annotation and visualization features for SSRs. Addressing this disparity and leveraging the de-novo millet genome assembly available from the NCBI, we have developed the Genomic SSR Millets Database (GSMDB; https://bioinfo.icgeb.res.in/gsmdb/). This open-access repository provides a web-based tool offering search functionalities and comprehensive details on 6.747645 million SSRs mined from the genomic sequences of seven millet species. The database, featuring unrestricted public access and JBrowse visualization, is a pioneering resource for the research community dedicated to advancing millet cultivars and related species. GSMDB holds immense potential to support myriad studies, including genetic diversity assessments, genetic mapping, marker-assisted selection, and comparative population investigations aiming to facilitate the millet breeding programs geared toward ensuring global food security. Database URL: https://bioinfo.icgeb.res.in/gsmdb/.

Epigenetic factors direct synergistic and antagonistic regulation of transposable elements in Arabidopsis

Arabidopsis (Arabidopsis thaliana) HISTONE DEACETYLASE 6 (HDA6) and HISTONE DEMETHYLASES LSD-LIKE 1 (LDL1) and LDL2 synergistically regulate the expression of long non-coding RNAs associated with H3Ac and H3K4me2. The underlying mechanisms of such highly coordinated interactions among genetic and epigenetic factors contributing to this collaborative regulation remain largely unclear. We analyzed all transposable elements (TEs) across the Arabidopsis genome and the individual and combined roles of HDA6 and LDL1/LDL2 by dissecting multilayered epigenomes and their association with transcription. Instead of an individual synergistic effect, we observed dual synergistic and antagonistic effects, which are positively associated with H3Ac and H3K4me2 while maintaining a negative but moderate association with DNA methylation. Specifically, 2 modes of synergistic regulation were discovered in TEs: 74% are primarily regulated by HDA6, with less dependence on LDL1/LDL2, and the remaining 26% are co-regulated by both. Between the 2 modes, we showed that HDA6 has a strong effect on TE silencing, whereas LDL1/LDL2 plays a weaker yet crucial role in co-regulation with HDA6. Our results led to a model of epigenomic regulation-the differential de-repression between the 2 modes of synergistic regulation of TEs was determined by H3Ac and H3K4me2 levels, where TEs are in accessible chromatins free of DNA methylation, and this open chromatin environment precedes transcriptional changes and epigenome patterning. Our results discovered unbalanced effects of genetic factors in synergistic regulation through delicately coordinated multilayered epigenomes and chromatin accessibility.

Genetic variation in the aquaporin TONOPLAST INTRINSIC PROTEIN 4;3 modulates maize cold tolerance

Cold stress is a major abiotic stress that threatens maize (Zea mays L.) production worldwide. Understanding the molecular mechanisms underlying cold tolerance is crucial for breeding resilient maize varieties. Tonoplast intrinsic proteins (TIPs) are a subfamily of aquaporins in plants. Here, we report that TIP family proteins are involved in maize cold tolerance. The expression of most TIP genes was responsive to cold stress. Overexpressing TIP2;1, TIP3;2 or TIP4;3 reduced the cold tolerance of maize seedlings, while loss-of-function mutants of TIP4;3 exhibited enhanced cold tolerance. Candidate gene-based association analysis revealed that a 328-bp transposon insertion in the promoter region of TIP4;3 was strongly associated with maize cold tolerance. This transposon insertion conferred cold tolerance by repressing TIP4;3 expression through increased methylation of its promoter region. Moreover, TIP4;3 was found to suppress stomatal closure and facilitate reactive oxygen species (ROS) accumulation under cold stress, thereby inhibiting the expression of cold-responsive genes, including DEHYDRATION-RESPONSIVE ELEMENT BINDING FACTOR 1 (DREB1) genes and a subset of peroxidase genes, ultimately attenuating maize cold tolerance. This study thus elucidates the mechanism underlying TIP-mediated cold tolerance and identifies a favourable TIP4;3 allele as a potential genetic resource for breeding cold-tolerant maize varieties.

Natural Diversity of Heat-Induced Transcription of Retrotransposons in Arabidopsis thaliana

Transposable elements (TEs) are major components of plant genomes, profoundly impacting the fitness of their hosts. However, technical bottlenecks have long hindered our mechanistic understanding of TEs. Using RNA-Seq and long-read sequencing with Oxford Nanopore Technologies' (ONT) direct cDNA sequencing, we analyzed the heat-induced transcription of TEs in three natural accessions of Arabidopsis thaliana (Cvi-0, Col-0, and Ler-1). In addition to the well-studied ONSEN retrotransposon family, we confirmed Copia-35 as a second heat-responsive retrotransposon family with particularly high activity in the relict accession Cvi-0. Our analysis revealed distinct expression patterns of individual TE copies and suggest different mechanisms regulating the GAG protein production in the ONSEN versus Copia-35 families. In addition, analogously to ONSEN, Copia-35 activation led to the upregulation of flanking genes such as APUM9 and potentially to the quantitative modulation of flowering time. ONT data allowed us to test the extent to which read-through formation is important in the regulation of adjacent genes. Unexpectedly, our results indicate that for both families, the upregulation of flanking genes is not predominantly directly initiated by transcription from their 3' long terminal repeats. These findings highlight the intraspecific expressional diversity linked to retrotransposon activation under stress.

Unlocking plant genetics with telomere-to-telomere genome assemblies

Contiguous genome sequence assemblies will help us to realize the full potential of crop translational genomics. Recent advances in sequencing technologies, especially long-read sequencing strategies, have made it possible to construct gapless telomere-to-telomere (T2T) assemblies, thus offering novel insights into genome organization and function. Plant genomes pose unique challenges, such as a continuum of ancient to recent polyploidy and abundant highly similar and long repetitive elements. Owing to progress in sequencing approaches, for most crop plants, chromosome-scale reference genome assemblies are available, but T2T assembly construction remains challenging. Here we describe methods for haplotype-resolved, gapless T2T assembly construction in plants, including various crop species. We outline the impact of T2T assemblies in elucidating the roles of repetitive elements in gene regulation, as well as in pangenomics, functional genomics, genome-assisted breeding and targeted genome manipulation. In conjunction with sequence-enriched germplasm repositories, T2T assemblies thus hold great promise for basic and applied plant sciences.

The regulatory potential of transposable elements in maize

Since their initial discovery in maize, transposable elements (TEs) have emerged as being integral to the evolution of maize, accounting for 80% of its genome. However, the repetitive nature of TEs has hindered our understanding of their regulatory potential. Here, we demonstrate that long-read chromatin fiber sequencing (Fiber-seq) permits the comprehensive annotation of the regulatory potential of maize TEs. We uncover that only 94 LTR retrotransposons contain the functional epigenetic architecture required for mobilization within maize leaves. This epigenetic architecture degenerates with evolutionary age, resulting in solo TE enhancers being preferentially marked by simultaneous hyper-CpG methylation and chromatin accessibility, an architecture markedly divergent from canonical enhancers. We find that TEs shape maize gene regulation by creating novel promoters within the TE itself as well as through TE-mediated gene amplification. Lastly, we uncover a pervasive epigenetic code directing TEs to specific loci, including that locus that sparked McClintock's discovery of TEs.

Harbinger transposon insertion in ethylene signaling gene leads to emergence of new sexual forms in cucurbits

In flowering plants, the predominant sexual morph is hermaphroditism, and the emergence of unisexuality is poorly understood. Using Cucumis melo (melon) as a model system, we explore the mechanisms driving sexual forms. We identify a spontaneous mutant exhibiting a transition from bisexual to unisexual male flower, and identify the causal mutation as a Harbinger transposon impairing the expression of Ethylene Insensitive 2 (CmEIN2) gene. Genetics and transcriptomic analysis reveal a dual role of CmEIN2 in both sex determination and fruit shape formation. Upon expression of CmACS11, EIN2 is recruited to repress the expression of the carpel inhibitor, CmWIP1. Subsequently, EIN2 is recruited to mediate stamina inhibition. Following the sex determination phase, EIN2 promotes fruit shape elongation. Genome-wide analysis reveals that Harbinger transposon mobilization is triggered by environmental cues, and integrates preferentially in active chromatin, particularly within promoter regions. Characterization of a large collection of melon germplasm points to active transpositions in the wild, compared to cultivated accessions. Our study underscores the association between chromatin dynamics and the temporal aspects of mobile genetic element insertions, providing valuable insights into plant adaptation and crop genome evolution.

Genome-wide identification and analyses of ZmAPY genes reveal their roles involved in maize development and abiotic stress responses

Apyrase is a class of enzyme that catalyzes the hydrolysis of nucleoside triphosphates/diphosphates (NTP/NDP), which widely involved in regulation of plant growth and stress responses. However, apyrase family genes in maize have not been identified, and their characteristics and functions are largely unknown. In this study, we identified 16 apyrases (named as ZmAPY1-ZmAPY16) in maize genome, and analyzed their phylogenetic relationships, gene structures, chromosomal distribution, upstream regulatory transcription factors and expression patterns. Analysis of the transcriptome database unveiled tissue-specific and abiotic stress-responsive expression of ZmAPY genes in maize. qPCR analysis further confirmed their responsiveness to drought, heat, and cold stresses. Association analyses indicated that variations of ZmAPY5 and ZmAPY16 may regulate maize agronomic traits and drought responses. Our findings shed light on the molecular characteristics and evolutionary history of maize apyrase genes, highlighting their roles in various biological processes and stress responses. This study forms a basis for further exploration of apyrase functions in maize.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01474-9.

Horticulture crop under pressure: Unraveling the impact of climate change on nutrition and fruit cracking

Climate change is increasingly affecting the nutritional content and structural integrity of horticultural crops, leading to challenges such as diminished fruit quality and the exacerbation of fruit cracking. This manuscript systematically explores the multifaceted impacts of these changes, with a particular focus on the nutritional quality and increased incidence of fruit cracking. An exhaustive review of current research identifies the critical role of transcription factors in mediating plant responses to climatic stressors, such as drought, temperature extremes, and saline conditions. The significance of transcription factors, including bHLH, bZIP, DOF, MDP, HD-ZIP, MYB, and ERF4, is highlighted in the development of fruit cracking, underscoring the genetic underpinnings behind stress-related phenotypic outcomes. The effectiveness of greenhouse structures in mitigating adverse climatic effects is evaluated, offering a strategic approach to sustain crop productivity amidst CO2 fluctuations and water scarcity, which are shown to influence plant physiology and lead to changes in fruit development, nutrient dynamics, and a heightened risk of cracking. Moreover, the manuscript delves into advanced breeding strategies and genetic engineering techniques, such as genome editing, to enhance crop resilience against climatic challenges. It also discusses adaptation strategies vital for sustainable horticulture, emphasizing the need to integrate novel genetic insights with controlled environment horticulture to counteract climate change's detrimental effects. The synthesis presented here underscores the urgent need for innovative breeding strategies aimed at developing resilient crop varieties that can withstand climatic uncertainty while preserving nutritional integrity.

Barley AGO4 proteins show overlapping functionality with distinct small RNA-binding properties in heterologous complementation

KEY MESSAGE

Barley AGO4 proteins complement expressional changes of epigenetically regulated genes in Arabidopsis ago4-3 mutant and show a distinct affinity for the 5' terminal nucleotide of small RNAs, demonstrating functional conservation and divergence. The function of Argonaute 4 (AGO4) in Arabidopsis thaliana has been extensively characterized; however, its role in monocots, which have large genomes abundantly supplemented with transposable elements (TEs), remains elusive. The study of barley AGO4 proteins can provide insights into the conserved aspects of RNA-directed DNA methylation (RdDM) and could also have further applications in the field of epigenetics or crop improvement. Bioinformatic analysis of RNA sequencing data identified two active AGO4 genes in barley, HvAGO4a and HvAGO4b. These genes function similar to AtAGO4 in an Arabidopsis heterologous complementation system, primarily binding to 24-nucleotide long small RNAs (sRNAs) and triggering methylation at specific target loci. Like AtAGO4, HvAGO4B exhibits a preference for binding sRNAs with 5' adenine residue, while also accepting 5' guanine, uracil, and cytosine residues. In contrast, HvAGO4A selectively binds only sRNAs with a 5' adenine residue. The diverse binding capacity of barley AGO4 proteins is reflected in TE-derived sRNAs and in their varying abundance. Both barley AGO4 proteins effectively restore the levels of extrachromosomal DNA and transcript abundancy of the heat-activated ONSEN retrotransposon to those observed in wild-type Arabidopsis plants. Our study provides insight into the distinct binding specificities and involvement in TE regulation of barley AGO4 proteins in Arabidopsis by heterologous complementation.

Transposition of HOPPLA in siRNA-deficient plants suggests a limited effect of the environment on retrotransposon mobility in Brachypodium distachyon

Long terminal repeat retrotransposons (LTR-RTs) are powerful mutagens regarded as a major source of genetic novelty and important drivers of evolution. Yet, the uncontrolled and potentially selfish proliferation of LTR-RTs can lead to deleterious mutations and genome instability, with large fitness costs for their host. While population genomics data suggest that an ongoing LTR-RT mobility is common in many species, the understanding of their dual role in evolution is limited. Here, we harness the genetic diversity of 320 sequenced natural accessions of the Mediterranean grass Brachypodium distachyon to characterize how genetic and environmental factors influence plant LTR-RT dynamics in the wild. When combining a coverage-based approach to estimate global LTR-RT copy number variations with mobilome-sequencing of nine accessions exposed to eight different stresses, we find little evidence for a major role of environmental factors in LTR-RT accumulations in B. distachyon natural accessions. Instead, we show that loss of RNA polymerase IV (Pol IV), which mediates RNA-directed DNA methylation in plants, results in high transcriptional and transpositional activities of RLC_BdisC024 (HOPPLA) LTR-RT family elements, and that these effects are not stress-specific. This work supports findings indicating an ongoing mobility in B. distachyon and reveals that host RNA-directed DNA methylation rather than environmental factors controls their mobility in this wild grass model.

Recent reactivation of a pathogenicity-associated transposable element is associated with major chromosomal rearrangements in a fungal wheat pathogen

Transposable elements (TEs) are key drivers of genomic variation contributing to recent adaptation in most species. Yet, the evolutionary origins and insertion dynamics within species remain poorly understood. We recapitulate the spread of the pathogenicity-associated Styx element across five species that last diverged ∼11 000 years ago. We show that the element likely originated in the Zymoseptoria fungal pathogen genus and underwent multiple independent reactivation events. Using a global 900-genome panel of the wheat pathogen Zymoseptoria tritici, we assess Styx copy number variation and identify renewed transposition activity in Oceania and South America. We show that the element can mobilize to create additional Styx copies in a four-generation pedigree. Importantly, we find that new copies of the element are not affected by genomic defenses suggesting minimal control against the element. Styx copies are preferentially located in recombination breakpoints and likely triggered multiple types of large chromosomal rearrangements. Taken together, we establish the origin, diversification and reactivation of a highly active TE with likely major consequences for chromosomal integrity and the expression of disease.

Preliminary characterization of biomolecular processes related to plasticity in Acyrthosiphonpisum

Global warming strongly impacts many organisms' development, distribution and population structure. This problem has attracted the attention of many scientists to understand and study its actual effects, especially on insects influenced by environmental temperatures. Aphids are a model for studies of the genetics and physiology of stress. Aphids are characterized by parthenogenetic reproduction, which limits the effects of recombination on evolutionary processes, and have shown resistance to various biotic and abiotic stresses. This study was based on the hypothesis that aphids have optimized, over time, genetic mechanisms capable to give them plasticity through genome modifications mediated by transposition. To understand and evaluate the effects of heat stress, the expression levels of transposases and methylases were analyzed in mothers and daughters. Our results show that after four days from the thermal shock, methylation decreases in both mothers and daughters, while transposition significantly increases in daughters, thus generating gene variability, essential for adaptation.

Transposable elements: multifunctional players in the plant genome

Transposable elements (TEs) are indispensable components of eukaryotic genomes that play diverse roles in gene regulation, recombination, and environmental adaptation. Their ability to mobilize within the genome leads to gene expression and DNA structure changes. TEs serve as valuable markers for genetic and evolutionary studies and facilitate genetic mapping and phylogenetic analysis. They also provide insight into how organisms adapt to a changing environment by promoting gene rearrangements that lead to new gene combinations. These repetitive sequences significantly impact genome structure, function and evolution. This review takes a comprehensive look at TEs and their applications in biotechnology, particularly in the context of plant biology, where they are now considered "genomic gold" due to their extensive functionalities. The article addresses various aspects of TEs in plant development, including their structure, epigenetic regulation, evolutionary patterns, and their use in gene editing and plant molecular markers. The goal is to systematically understand TEs and shed light on their diverse roles in plant biology.

The maize ATP-binding cassette (ABC) transporter ZmMRPA6 confers cold and salt stress tolerance in plants

KEY MESSAGE

ZmMRPA6 was cloned and characterized as the first ATP-binding cassette (ABC) transporter in maize to be proven to participate in cold and salt tolerance. Homologous genes AtABCC4 and AtABCC14 of ZmMRPA6 also responded to salt stress. ATP-binding cassette (ABC) proteins are major transmembrane transporters that play significant roles in plant development against various abiotic stresses. However, available information regarding stress-related ABC genes in maize is minimal. In this study, a maize ABC transporter gene, ZmMRPA6, was identified through genome-wide association analysis (GWAS) for cold tolerance in maize seeds germination and functionally characterized. During germination and seedling stages, the zmmrpa6 mutant exhibited enhanced resistance to cold or salt stress. Mutated of ZmMRPA6 did not affect the expression of downstream response genes related cold or salt response at the transcriptional level. Mass spectrometry analysis revealed that most of the differential proteins between zmmrpa6 and wild-type plants were involved in response to stress process including oxidative reduction, hydrolase activity, small molecule metabolism, and photosynthesis process. Meanwhile, the plants which lack the ZmMRPA6 homologous genes AtABCC4 or AtABCC14 were sensitive to salt stress in Arabidopsis. These results indicated that ZmMRPA6 and its homologous genes play a conserved role in cold and salt stress, and functional differentiation occurs in monocotyledonous and dicotyledonous plants. In summary, these findings dramatically improved our understanding of the function of ABC transporters resistance to abiotic stresses in plants.

A multi-organ maize metabolic model connects temperature stress with energy production and reducing power generation

Climate change has adversely affected maize productivity. Thereby, a holistic understanding of metabolic crosstalk among its organs is important to address this issue. Thus, we reconstructed the first multi-organ maize metabolic model, iZMA6517, and contextualized it with heat and cold stress transcriptomics data using expression distributed reaction flux measurement (EXTREAM) algorithm. Furthermore, implementing metabolic bottleneck analysis on contextualized models revealed differences between these stresses. While both stresses had reducing power bottlenecks, heat stress had additional energy generation bottlenecks. We also performed thermodynamic driving force analysis, revealing thermodynamics-reducing power-energy generation axis dictating the nature of temperature stress responses. Thus, a temperature-tolerant maize ideotype can be engineered by leveraging the proposed thermodynamics-reducing power-energy generation axis. We experimentally inoculated maize root with a beneficial mycorrhizal fungus, Rhizophagus irregularis, and as a proof-of-concept demonstrated its efficacy in alleviating temperature stress. Overall, this study will guide the engineering effort of temperature stress-tolerant maize ideotypes.

RNA pseudouridine modification in plants

Pseudouridine is one of the well-known chemical modifications in various RNA species. Current advances to detect pseudouridine show that the pseudouridine landscape is dynamic and affects multiple cellular processes. Although our understanding of this post-transcriptional modification mainly depends on yeast and human models, the recent findings provide strong evidence for the critical role of pseudouridine in plants. Here, we review the current knowledge of pseudouridine in plant RNAs, including its synthesis, degradation, regulatory mechanisms, and functions. Moreover, we propose future areas of research on pseudouridine modification in plants.

The first high-quality genome assembly and annotation of Patiria pectinifera

The blue bat star, a highly adaptive species in the East Sea of Korea, has displayed remarkable success in adapting to recent climate change. The genetic mechanisms behind this success were not well-understood, prompting our report on the first chromosome-level assembly of the Patiria genus. We assembled the genome using Nanopore and Illumina sequences, yielding a total length of 615 Mb and a scaffold N50 of 24,204,423 bp. Hi-C analysis allowed us to anchor the scaffold sequences onto 22 pseudochromosomes. K-mer based analysis revealed 5.16% heterozygosity rate of the genome, higher than any previously reported echinoderm species. Our transposable element analysis exposed a substantial number of genome-wide retrotransposons and DNA transposons. These results offer valuable resources for understanding the evolutionary mechanisms behind P. pectinifera's successful adaptation in fluctuating environments.

Regulation of plant epigenetic memory in response to cold and heat stress: towards climate resilient agriculture

Plants have evolved to adapt and grow in hot and cold climatic conditions. Some also adapt to daily and seasonal temperature changes. Epigenetic modifications play an important role in regulating plant tolerance under such conditions. DNA methylation and post-translational modifications of histone proteins influence gene expression during plant developmental stages and under stress conditions, including cold and heat stress. While short-term modifications are common, some modifications may persist and result in stress memory that can be inherited by subsequent generations. Understanding the mechanisms of epigenomes responding to stress and the factors that trigger stress memory is crucial for developing climate-resilient agriculture, but such an integrated view is currently limited. This review focuses on the plant epigenetic stress memory during cold and heat stress. It also discusses the potential of machine learning to modify stress memory through epigenetics to develop climate-resilient crops.

From parasites to partners: exploring the intricacies of host-transposon dynamics and coevolution

Transposable elements, often referred to as "jumping genes," have long been recognized as genomic parasites due to their ability to integrate and disrupt normal gene function and induce extensive genomic alterations, thereby compromising the host's fitness. To counteract this, the host has evolved a plethora of mechanisms to suppress the activity of the transposons. Recent research has unveiled the host-transposon relationships to be nuanced and complex phenomena, resulting in the coevolution of both entities. Transposition increases the mutational rate in the host genome, often triggering physiological pathways such as immune and stress responses. Current gene transfer technologies utilizing transposable elements have potential drawbacks, including off-target integration, induction of mutations, and modifications of cellular machinery, which makes an in-depth understanding of the host-transposon relationship imperative. This review highlights the dynamic interplay between the host and transposable elements, encompassing various factors and components of the cellular machinery. We provide a comprehensive discussion of the strategies employed by transposable elements for their propagation, as well as the mechanisms utilized by the host to mitigate their parasitic effects. Additionally, we present an overview of recent research identifying host proteins that act as facilitators or inhibitors of transposition. We further discuss the evolutionary outcomes resulting from the genetic interactions between the host and the transposable elements. Finally, we pose open questions in this field and suggest potential avenues for future research.

Identification and characterization of RuvBL DNA helicase genes for tolerance against abiotic stresses in bread wheat (Triticum aestivum L.) and related species

Recombination UVB (sensitivity) like (RuvBL) helicase genes represent a conserved family of genes, which are known to be involved in providing tolerance against abiotic stresses like heat and drought. We identified nine wheat RuvBL genes, one each on nine different chromosomes, belonging to homoeologous groups 2, 3, and 4. The lengths of genes ranged from 1647 to 2197 bp and exhibited synteny with corresponding genes in related species including Ae. tauschii, Z. mays, O. sativa, H. vulgare, and B. distachyon. The gene sequences were associated with regulatory cis-elements and transposable elements. Two genes, namely TaRuvBL1a-4A and TaRuvBL1a-4B, also carried targets for a widely known miRNA, tae-miR164. Gene ontology revealed that these genes were closely associated with ATP-dependent formation of histone acetyltransferase complex. Analysis of the structure and function of RuvBL proteins revealed that the proteins were localized mainly in the cytoplasm. A representative gene, namely TaRuvBL1a-4A, was also shown to be involved in protein-protein interactions with ten other proteins. On the basis of phylogeny, RuvBL proteins were placed in two sub-divisions, namely RuvBL1 and RuvBL2, which were further classified into clusters and sub-clusters. In silico studies suggested that these genes were differentially expressed under heat/drought. The qRT-PCR analysis confirmed that expression of TaRuvBL genes differed among wheat cultivars, which differed in the level of thermotolerance. The present study advances our understanding of the biological role of wheat RuvBL genes and should help in planning future studies on RuvBL genes in wheat including use of RuvBL genes in breeding thermotolerant wheat cultivars.

A dolabralexin-deficient mutant provides insight into specialized diterpenoid metabolism in maize

Two major groups of specialized metabolites in maize (Zea mays), termed kauralexins and dolabralexins, serve as known or predicted diterpenoid defenses against pathogens, herbivores, and other environmental stressors. To consider the physiological roles of the recently discovered dolabralexin pathway, we examined dolabralexin structural diversity, tissue-specificity, and stress-elicited production in a defined biosynthetic pathway mutant. Metabolomics analyses support a larger number of dolabralexin pathway products than previously known. We identified dolabradienol as a previously undetected pathway metabolite and characterized its enzymatic production. Transcript and metabolite profiling showed that dolabralexin biosynthesis and accumulation predominantly occur in primary roots and show quantitative variation across genetically diverse inbred lines. Generation and analysis of CRISPR-Cas9-derived loss-of-function Kaurene Synthase-Like 4 (Zmksl4) mutants demonstrated dolabralexin production deficiency, thus supporting ZmKSL4 as the diterpene synthase responsible for the conversion of geranylgeranyl pyrophosphate precursors into dolabradiene and downstream pathway products. Zmksl4 mutants further display altered root-to-shoot ratios and root architecture in response to water deficit. Collectively, these results demonstrate dolabralexin biosynthesis via ZmKSL4 as a committed pathway node biochemically separating kauralexin and dolabralexin metabolism, and suggest an interactive role of maize dolabralexins in plant vigor during abiotic stress.

Genomic imbalance modulates transposable element expression in maize

Genomic imbalance refers to the more severe phenotypic consequences of changing part of a chromosome compared with the whole genome set. Previous genome imbalance studies in maize have identified prevalent inverse modulation of genes on the unvaried chromosomes (trans) with both the addition or subtraction of chromosome arms. Transposable elements (TEs) comprise a substantial fraction of the genome, and their reaction to genomic imbalance is therefore of interest. Here, we analyzed TE expression using RNA-seq data of aneuploidy and ploidy series and found that most aneuploidies showed an inverse modulation of TEs, but reductions in monosomy and increases in disomy and trisomy were also common. By contrast, the ploidy series showed little TE modulation. The modulation of TEs and genes in the same experimental group were compared, and TEs showed greater modulation than genes, especially in disomy. Class I and II TEs were differentially modulated in most aneuploidies, and some superfamilies in each TE class also showed differential modulation. Finally, the significantly upregulated TEs in three disomies (TB-7Lb, TB9Lc, and TB-10L19) did not increase the proportion of adjacent gene expression when compared with non-differentially expressed TEs, indicating that modulations of TEs do not compound the effect on genes. These results suggest that the prevalent inverse TE modulation in aneuploidy results from stoichiometric upset of the regulatory machinery used by TEs, similar to the response of core genes to genomic imbalance.

R3-MYB repressor Mybr97 is a candidate gene associated with the Anthocyanin3 locus and enhanced anthocyanin accumulation in maize

Anthocyanin3 inhibits the anthocyanin and monolignol pathways in maize. Transposon-tagging, RNA-sequencing, and GST-pulldown assays determine Anthocyanin3 may be R3-MYB repressor gene Mybr97. Anthocyanins are colorful molecules receiving recent attention due to their numerous health benefits and applications as natural colorants and nutraceuticals. Purple corn is being investigated as a more economical source of anthocyanins. Anthocyanin3 (A3) is a known recessive intensifier of anthocyanin pigmentation in maize. In this study, anthocyanin content was elevated 100-fold in recessive a3 plants. Two approaches were used to discover candidates involved with the a3 intense purple plant phenotype. First, a large-scale transposon-tagging population was created with a Dissociation (Ds) insertion in the nearby Anthocyanin1 gene. A de novo a3-m1::Ds mutant was generated, and the transposon insertion was found to be located in the promoter of Mybr97, which has homology to R3-MYB repressor CAPRICE in Arabidopsis. Second, a bulked segregant RNA-sequencing population found expression differences between pools of green A3 plants and purple a3 plants. All characterized anthocyanin biosynthetic genes were upregulated in a3 plants along with several genes of the monolignol pathway. Mybr97 was highly downregulated in a3 plants, suggesting its role as a negative regulator of the anthocyanin pathway. Photosynthesis-related gene expression was reduced in a3 plants through an unknown mechanism. Numerous transcription factors and biosynthetic genes were also upregulated and need further investigation. Mybr97 may inhibit anthocyanin synthesis by associating with basic helix-loop helix transcription factors like Booster1. Overall, Mybr97 is the most likely candidate gene for the A3 locus. A3 has a profound effect on the maize plant and has many favorable implications for crop protection, human health, and natural colorant production.

Evolutionary Analysis of Respiratory Burst Oxidase Homolog (RBOH) Genes in Plants and Characterization of ZmRBOHs

The respiratory burst oxidase homolog (RBOH), as the key producer of reactive oxygen species (ROS), plays an essential role in plant development. In this study, a bioinformatic analysis was performed on 22 plant species, and 181 RBOH homologues were identified. A typical RBOH family was identified only in terrestrial plants, and the number of RBOHs increased from non-angiosperms to angiosperms. Whole genome duplication (WGD)/segmental duplication played a key role in RBOH gene family expansion. Amino acid numbers of 181 RBOHs ranged from 98 to 1461, and the encoded proteins had molecular weights from 11.1 to 163.6 kDa, respectively. All plant RBOHs contained a conserved NADPH_Ox domain, while some of them lacked the FAD_binding_8 domain. Plant RBOHs were classified into five main subgroups by phylogenetic analysis. Most RBOH members in the same subgroup showed conservation in both motif distribution and gene structure composition. Fifteen ZmRBOHs were identified in maize genome and were positioned in eight maize chromosomes. A total of three pairs of orthologous genes were found in maize, including ZmRBOH6/ZmRBOH8, ZmRBOH4/ZmRBOH10 and ZmRBOH15/ZmRBOH2. A Ka/Ks calculation confirmed that purifying selection was the main driving force in their evolution. ZmRBOHs had typical conserved domains and similar protein structures. cis-element analyses together with the expression profiles of the ZmRBOH genes in various tissues and stages of development suggested that ZmRBOH was involved in distinct biological processes and stress responses. Based on the RNA-Seq data and qRT-PCR analysis, the transcriptional response of ZmRBOH genes was examined under various abiotic stresses, and most of ZmRBOH genes were up-regulated by cold stress. These findings provide valuable information for further revealing the biological roles of ZmRBOH genes in plant development and abiotic stress responses.

Chromosomal-level genome assembly of the high-quality Xian/Indica rice (Oryza sativa L.) Xiangyaxiangzhan

The indica rice variety XYXZ carries elite traits including appearance and eating quality. Here, we report the de novo assembly of XYXZ using Illumine paired-end whole-genome shotgun sequencing and Nanopore sequencing. We annotated 39,722 protein-coding genes in the 395.04 Mb assembly. In comparison to other cultivars, XYXZ showed a larger gene size including the transcripts and introns, and more exons per gene. And hundreds of ultra-long genes were also detected. A total of 4362 complete LTRs were annotated, and among them, many were located next to or in protein-coding genes including several genes related to rice quality. We observed the different distributions of LTRs in these genes among XYXZ, Nipponbare, and R498, implying these LTRs might potentially affect expressions of the proximal genes and rice quality. Overall, This chromosome-length genome assembly of XYXZ provides a valuable resource for gene discovery, genetic variation and evolution, and the breeding of high-quality rice.

Genome-wide identification and characterization of exapted transposable elements in the large genome of sunflower (Helianthus annuus L.)

Transposable elements (TEs) are an important source of genome variability, playing many roles in the evolution of eukaryotic species. Besides well-known phenomena, TEs may undergo the exaptation process and generate the so-called exapted transposable element genes (ETEs). Here we present a genome-wide survey of ETEs in the large genome of sunflower (Helianthus annuus L.), in which the massive amount of TEs, provides a significant source for exaptation. A library of sunflower TEs was used to build TE-specific Hidden Markov Model profiles, to search for all available sunflower gene products. In doing so, 20 016 putative ETEs were identified and further investigated for the characteristics that distinguish TEs from genes, leading to the validation of 3530 ETEs. The analysis of ETEs transcription patterns under different stress conditions showed a differential regulation triggered by treatments mimicking biotic and abiotic stress; furthermore, the distribution of functional domains of differentially regulated ETEs revealed a relevant presence of domains involved in many aspects of cellular functions. A comparative genomic investigation was performed including species representative of Asterids and appropriate outgroups: the bulk of ETEs that resulted were specific to the sunflower, while few ETEs presented orthologues in the genome of all analyzed species, making the hypothesis of a conserved function. This study highlights the crucial role played by exaptation, actively contributing to species evolution.

Graph-based pan-genomes: increased opportunities in plant genomics

Due to the development of sequencing technology and the great reduction in sequencing costs, an increasing number of plant genomes have been assembled, and numerous genomes have revealed large amounts of variations. However, a single reference genome does not allow the exploration of species diversity, and therefore the concept of pan-genome was developed. A pan-genome is a collection of all sequences available for a species, including a large number of consensus sequences, large structural variations, and small variations including single nucleotide polymorphisms and insertions/deletions. A simple linear pan-genome does not allow these structural variations to be intuitively characterized, so graph-based pan-genomes have been developed. These pan-genomes store sequence and structural variation information in the form of nodes and paths to store and display species variation information in a more intuitive manner. The key role of graph-based pan-genomes is to expand the coordinate system of the linear reference genome to accommodate more regions of genetic diversity. Here, we review the origin and development of graph-based pan-genomes, explore their application in plant research, and further highlight the application of graph-based pan-genomes for future plant breeding.

PlanTEnrichment: A How-to Guide on Rapid Identification of Transposable Elements Associated with Regions of Interest in Select Plant Genomes

Transposable elements (TEs) are repeat elements that can relocate or create novel copies of themselves in the genome and contribute to genomic complexity and expansion, via events such as chromosome recombination or regulation of gene expression. However, given the large number of such repeats across the genome, identifying repeats of interest can be a challenge in even well-annotated genomes, especially in more complex, TE-rich plant genomes. Here, we describe a protocol for PlanTEnrichment, a database we created comprising information on 11 plant genomes to analyze stress-associated TEs using publicly available data. By selecting a genome and providing a list of genes or genomic regions whose TE associations the user wants to identify, the user can rapidly obtain TE subfamilies found near the provided regions, as well as their superfamily and class, and the enrichment values of the repeats. The results also provide the locations of individual repeat instances found, alongside the input regions or genes they are associated with, and a bar graph of the top ten most significant repeat subfamilies identified. PlanTEnrichment is freely available at http://tools.ibg.deu.edu.tr/plantenrichment/ and can be used by researchers with rudimentary or no proficiency in computational analysis of TE elements, allowing for expedience in the identification of TEs of interest and helping further our understanding of the potential contributions of TEs in plant genomes.

A long non-coding RNA PelncRNA1 is involved in Phyllostachys edulis response to UV-B stress

Phyllostachys edulis (moso bamboo) is China's most widespread bamboo species, with significant economic and ecological values. Long non-coding RNA (lncRNA) is a type of regulatory RNA that is longer than 200 nucleotides and incapable of encoding proteins, and is frequently involved in regulating biotic and abiotic stress and plant development. However, the biological functions of lncRNA in moso bamboo are unknown. In this study, a lncRNA (named PelncRNA1) differentially expressed following UV-B treatment was discovered in the whole transcriptome sequencing database of moso bamboo. The target genes were filtered and defined by correlation analysis of PelncRNA1 and gene expression pattern. The expression levels of PelncRNA1 and its target genes were verified using qRT-PCR. The results demonstrated that the expression levels of PelncRNA1 and its target genes increased during UV-B treatment. In Arabidopsis transgenic seedlings and moso bamboo protoplasts, PelncRNA1 was discovered to influence the expression of its target genes when overexpressed. In addition, transgenic Arabidopsis showed higher tolerance to UV-B stress. These results suggest that PelncRNA1 and its target genes are involved in the response of moso bamboo to UV-B stress. The novel findings would contribute to our understanding of how lncRNAs regulate the response to abiotic stresses in moso bamboo.

A review of strategies used to identify transposition events in plant genomes

Transposable elements (TEs) were initially considered redundant and dubbed 'junk DNA'. However, more recently they were recognized as an essential element of genome plasticity. In nature, they frequently become active upon exposition of the host to stress conditions. Even though most transposition events are neutral or even deleterious, occasionally they may happen to be beneficial, resulting in genetic novelty providing better fitness to the host. Hence, TE mobilization may promote adaptability and, in the long run, act as a significant evolutionary force. There are many examples of TE insertions resulting in increased tolerance to stresses or in novel features of crops which are appealing to the consumer. Possibly, TE-driven de novo variability could be utilized for crop improvement. However, in order to systematically study the mechanisms of TE/host interactions, it is necessary to have suitable tools to globally monitor any ongoing TE mobilization. With the development of novel potent technologies, new high-throughput strategies for studying TE dynamics are emerging. Here, we present currently available methods applied to monitor the activity of TEs in plants. We divide them on the basis of their operational principles, the position of target molecules in the process of transposition and their ability to capture real cases of actively transposing elements. Their possible theoretical and practical drawbacks are also discussed. Finally, conceivable strategies and combinations of methods resulting in an improved performance are proposed.

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.

The ecology of palm genomes: repeat-associated genome size expansion is constrained by aridity

Genome size varies 2400-fold across plants, influencing their evolution through changes in cell size and cell division rates which impact plants' environmental stress tolerance. Repetitive element expansion explains much genome size diversity, and the processes structuring repeat 'communities' are analogous to those structuring ecological communities. However, which environmental stressors influence repeat community dynamics has not yet been examined from an ecological perspective. We measured genome size and leveraged climatic data for 91% of genera within the ecologically diverse palm family (Arecaceae). We then generated genomic repeat profiles for 141 palm species, and analysed repeats using phylogenetically informed linear models to explore relationships between repeat dynamics and environmental factors. We show that palm genome size and repeat 'community' composition are best explained by aridity. Specifically, Ty3-gypsy and TIR elements were more abundant in palm species from wetter environments, which generally had larger genomes, suggesting amplification. By contrast, Ty1-copia and LINE elements were more abundant in drier environments. Our results suggest that water stress inhibits repeat expansion through selection on upper genome size limits. However, elements that may associate with stress-response genes (e.g. Ty1-copia) have amplified in arid-adapted palm species. Overall, we provide novel evidence of climate influencing the assembly of repeat 'communities'.

Transcription factor ZmWRKY20 interacts with ZmWRKY115 to repress expression of ZmbZIP111 for salt tolerance in maize

Maize (Zea mays) is an important cereal crop worldwide. However, its yield and quality are adversely affected by salt stress resulting from soil hypersalinity. Exploring the regulatory mechanisms of stress responses is of vital importance to increase maize seed production. In the present study, we screened ethyl methanesulfonate-induced maize mutants and identified a salt-tolerant mutant. A single base was mutated in ZmWRKY20, leading to the formation of a truncated protein variant. A detailed phenotypic analysis revealed that this mutant had significantly higher resistance to wilting and lower reactive oxygen species levels than the inbred line B73. ZmWRKY20 showed transcriptional activity in yeast and specifically bound W-boxes according to the results of our yeast one-hybrid, electrophoretic mobility shift, and dual-luciferase assays. Overexpression of ZmWRKY20 decreased salt tolerance in maize. Transcriptome profiling revealed that ZmWRKY20 overexpression extensively reprogrammed genes involved in regulating defense and oxidation-reduction responses. The results substantiate that ZmWRKY20 is directly targeted to the basic leucine zipper (bZIP) motif in the transcription factor ZmbZIP111. It was also verified that ZmWRKY20 interacts with ZmWRKY115 and both proteins act jointly to enhance ZmbZIP111 repression. The results indicate that the ZmWRKY20 and ZmWRKY115 transcription factors interact in the nucleus, leading to repression of ZmbZIP111 expression by directly binding its promoter, and increase the sensitivity of maize seedlings to salt stress. The current study improves our understanding of the complicated responses of maize to salt stress.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant-pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.

Hybrid Assembly and Annotation of the Genome of the Indian Punica granatum, a Superfood

The wonder fruit pomegranate (Punica granatum, family Lythraceae) is one of India’s economically important fruit crops that can grow in different agro-climatic conditions ranging from tropical to temperate regions. This study reports high-quality de novo draft hybrid genome assembly of diploid Punica cultivar “Bhagwa” and identifies its genomic features. This cultivar is most common among the farmers due to its high sustainability, glossy red color, soft seed, and nutraceutical properties with high market value. The draft genome assembly is about 361.76 Mb (N50 = 40 Mb), ∼9.0 Mb more than the genome size estimated by flow cytometry. The genome is 90.9% complete, and only 26.68% of the genome is occupied by transposable elements and has a relative abundance of 369.93 SSRs/Mb of the genome. A total of 30,803 proteins and their putative functions were predicted. Comparative whole-genome analysis revealed Eucalyptus grandis as the nearest neighbor. KEGG-KASS annotations indicated an abundance of genes involved in the biosynthesis of flavonoids, phenylpropanoids, and secondary metabolites, which are responsible for various medicinal properties of pomegranate, including anticancer, antihyperglycemic, antioxidant, and anti-inflammatory activities. The genome and gene annotations provide new insights into the pharmacological properties of the secondary metabolites synthesized in pomegranate. They will also serve as a valuable resource in mining biosynthetic pathways for key metabolites, novel genes, and variations associated with disease resistance, which can facilitate the breeding of new varieties with high yield and superior quality.

Spatial Features and Functional Implications of Plant 3D Genome Organization

The advent of high-throughput sequencing-based methods for chromatin conformation, accessibility, and immunoprecipitation assays has been a turning point in 3D genomics. Altogether, these new tools have been pushing upward the interpretation of pioneer cytogenetic evidence for a higher order in chromatin packing. Here, we review the latest development in our understanding of plant spatial genome structures and different levels of organization and discuss their functional implications. Then, we spotlight the complexity of organellar (i.e., mitochondria and plastids) genomes and discuss their 3D packing into nucleoids. Finally, we propose unaddressed research axes to investigate functional links between chromatin-like dynamics and transcriptional regulation within organellar nucleoids.

Comprehensive analysis of both long and short read transcriptomes of a clonal and a seed-propagated model species reveal the prerequisites for transcriptional activation of autonomous and non-autonomous transposons in plants

Background

Transposable element (TE) transcription is a precursor to its mobilisation in host genomes. However, the characteristics of expressed TE loci, the identification of self-competent transposon loci contributing to new insertions, and the genomic conditions permitting their mobilisation remain largely unknown.

Results

Using Vitis vinifera embryogenic callus, we explored the impact of biotic stressors on transposon transcription through the exposure of the callus to live cultures of an endemic grapevine yeast, Hanseniaspora uvarum. We found that only 1.7–2.5% of total annotated TE loci were transcribed, of which 5–10% of these were full-length, and the expressed TE loci exhibited a strong location bias towards expressed genes. These trends in transposon transcription were also observed in RNA-seq data from Arabidopsis thaliana wild-type plants but not in epigenetically compromised Arabidopsis ddm1 mutants. Moreover, differentially expressed TE loci in the grapevine tended to share expression patterns with co-localised differentially expressed genes. Utilising nanopore cDNA sequencing, we found a strong correlation between the inclusion of intronic TEs in gene transcripts and the presence of premature termination codons in these transcripts. Finally, we identified low levels of full-length transcripts deriving from structurally intact TE loci in the grapevine model.

Conclusion

Our observations in two disparate plant models representing clonally and seed propagated plant species reveal a closely connected transcriptional relationship between TEs and co-localised genes, particularly when epigenetic silencing is not compromised. We found that the stress treatment alone was insufficient to induce large-scale full-length transcription from structurally intact TE loci, a necessity for non-autonomous and autonomous mobilisation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00271-5.

Patterns of selection in the evolution of a transposable element

Transposable elements are a major component of most eukaryotic genomes. Here, we present a new approach which allows us to study patterns of natural selection in the evolution of transposable elements over short time scales. The method uses the alignment of all elements with intact gag/pol genes of a transposable element family from a single genome. We predict that the ratio of nonsynonymous to synonymous variants in the alignment should decrease as a function of the frequency of the variants, because elements with nonsynonymous variants that reduce transposition will have fewer progeny. We apply our method to Sirevirus long-terminal repeat retrotransposons that are abundant in maize and other plant species and show that nonsynonymous to synonymous variants declines as variant frequency increases, indicating that negative selection is acting strongly on the Sirevirus genome. The asymptotic value of nonsynonymous to synonymous variants suggests that at least 85% of all nonsynonymous mutations in the transposable element reduce transposition. Crucially, these patterns in nonsynonymous to synonymous variants are only predicted to occur if the gene products from a particular transposable element insertion preferentially promote the transposition of the same insertion. Overall, by using large numbers of intact elements, this study sheds new light on the selective processes that act on transposable elements.

Transcriptional Contribution of Transposable Elements in Relation to Salinity Conditions in Teleosts and Silencing Mechanisms Involved

Fish are an interesting taxon comprising species adapted to a wide range of environments. In this work, we analyzed the transcriptional contribution of transposable elements (TEs) in the gill transcriptomes of three fish species exposed to different salinity conditions. We considered the giant marbled eel Anguilla marmorata and the chum salmon Oncorhynchus keta, both diadromous, and the marine medaka Oryzias melastigma, an euryhaline organism sensu stricto. Our analyses revealed an interesting activity of TEs in the case of juvenile eels, commonly adapted to salty water, when exposed to brackish and freshwater conditions. Moreover, the expression assessment of genes involved in TE silencing mechanisms (six in heterochromatin formation, fourteen known to be part of the nucleosome remodeling deacetylase (NuRD) complex, and four of the Argonaute subfamily) unveiled that they are active. Finally, our results evidenced for the first time a krüppel-associated box (KRAB)-like domain specific to actinopterygians that, together with TRIM33, might allow the functioning of NuRD complex also in fish species. The possible interaction between these two proteins was supported by structural prediction analyses.

Genome-wide analyses of introgression between two sympatric Asian oak species

Introgression can be an important source of new alleles for adaption under rapidly changing environments, perhaps even more important than standing variation. Though introgression has been extensively studied in many plants and animals, key questions on the underlying mechanisms of introgression still remain unanswered. In particular, we are yet to determine the genomic distribution of introgressed regions along the genome; whether the extent and patterns of introgression are influenced by ecological factors; and when and how introgression contributes to adaptation. Here, we generated high-quality genomic resources for two sympatric widespread Asian oak species, Quercus acutissima and Q. variabilis, sampled in multiple forests to study introgression between them. We show that introgressed regions are broadly distributed across the genome. Introgression was affected by genetic divergence between pairs of populations and by the similarity of the environments in which they live-populations occupying similar ecological sites tended to share the same introgressed regions. Introgressed genomic footprints of adaptation were preferentially located in regions with suppressed recombination rate. Introgression probably confers adaptation in these oak populations by introducing allelic variation in cis-regulatory elements, in particular through transposable element insertions, thereby altering the regulation of genes related to stress. Our results provide new avenues of research for uncovering mechanisms of adaptation due to hybridization in sympatric species.

Meta-QTL analysis and candidate genes identification for various abiotic stresses in maize (Zea mays L.) and their implications in breeding programs

Global climate change leads to the concurrence of a number of abiotic stresses including moisture stress (drought, waterlogging), temperature stress (heat, cold), and salinity stress, which are the major factors affecting maize production. To develop abiotic stress tolerance in maize, many quantitative trait loci (QTL) have been identified, but very few of them have been utilized successfully in breeding programs. In this context, the meta-QTL analysis of the reported QTL will enable the identification of stable/real QTL which will pave a reliable way to introgress these QTL into elite cultivars through marker-assisted selection. In this study, a total of 542 QTL were summarized from 33 published studies for tolerance to different abiotic stresses in maize to conduct meta-QTL analysis using BiomercatorV4.2.3. Among those, only 244 major QTL with more than 10% phenotypic variance were preferably utilised to carry out meta-QTL analysis. In total, 32 meta-QTL possessing 1907 candidate genes were detected for different abiotic stresses over diverse genetic and environmental backgrounds. The MQTL2.1, 5.1, 5.2, 5.6, 7.1, 9.1, and 9.2 control different stress-related traits for combined abiotic stress tolerance. The candidate genes for important transcription factor families such as ERF, MYB, bZIP, bHLH, NAC, LRR, ZF, MAPK, HSP, peroxidase, and WRKY have been detected for different stress tolerances. The identified meta-QTL are valuable for future climate-resilient maize breeding programs and functional validation of candidate genes studies, which will help to deepen our understanding of the complexity of these abiotic stresses.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01294-9.

Genomes shed light on the evolution of Begonia, a mega-diverse genus

Clarifying the evolutionary processes underlying species diversification and adaptation is a key focus of evolutionary biology. Begonia (Begoniaceae) is one of the most species-rich angiosperm genera with c. 2000 species, most of which are shade-adapted. Here, we present chromosome-scale genome assemblies for four species of Begonia (B. loranthoides, B. masoniana, B. darthvaderiana and B. peltatifolia), and whole genome shotgun data for an additional 74 Begonia representatives to investigate lineage evolution and shade adaptation of the genus. The four genome assemblies range in size from 331.75 Mb (B. peltatifolia) to 799.83 Mb (B. masoniana), and harbor 22 059-23 444 protein-coding genes. Synteny analysis revealed a lineage-specific whole-genome duplication (WGD) that occurred just before the diversification of Begonia. Functional enrichment of gene families retained after WGD highlights the significance of modified carbohydrate metabolism and photosynthesis possibly linked to shade adaptation in the genus, which is further supported by expansions of gene families involved in light perception and harvesting. Phylogenomic reconstructions and genomics studies indicate that genomic introgression has also played a role in the evolution of Begonia. Overall, this study provides valuable genomic resources for Begonia and suggests potential drivers underlying the diversity and adaptive evolution of this mega-diverse clade.

A B73×Palomero Toluqueño mapping population reveals local adaptation in Mexican highland maize

Generations of farmer selection in the central Mexican highlands have produced unique maize varieties adapted to the challenges of the local environment. In addition to possessing great agronomic and cultural value, Mexican highland maize represents a good system for the study of local adaptation and acquisition of adaptive phenotypes under cultivation. In this study, we characterize a recombinant inbred line population derived from the B73 reference line and the Mexican highland maize variety Palomero Toluqueño. B73 and Palomero Toluqueño showed classic rank-changing differences in performance between lowland and highland field sites, indicative of local adaptation. Quantitative trait mapping identified genomic regions linked to effects on yield components that were conditionally expressed depending on the environment. For the principal genomic regions associated with ear weight and total kernel number, the Palomero Toluqueño allele conferred an advantage specifically in the highland site, consistent with local adaptation. We identified Palomero Toluqueño alleles associated with expression of characteristic highland traits, including reduced tassel branching, increased sheath pigmentation and the presence of sheath macrohairs. The oligogenic architecture of these three morphological traits supports their role in adaptation, suggesting they have arisen from consistent directional selection acting at distinct points across the genome. We discuss these results in the context of the origin of phenotypic novelty during selection, commenting on the role of de novo mutation and the acquisition of adaptive variation by gene flow from endemic wild relatives.

Genome-Wide Identification and Expression Analysis of Pseudouridine Synthase Family in Arabidopsis and Maize

Pseudouridine (Ψ), the isomer of uridine (U), is the most abundant type of RNA modification, which is crucial for gene regulation in various cellular processes. Pseudouridine synthases (PUSs) are the key enzymes for the U-to-Ψ conversion. However, little is known about the genome-wide features and biological function of plant PUSs. In this study, we identified 20 AtPUSs and 22 ZmPUSs from Arabidopsis and maize (Zea mays), respectively. Our phylogenetic analysis indicated that both AtPUSs and ZmPUSs could be clustered into six known subfamilies: RluA, RsuA, TruA, TruB, PUS10, and TruD. RluA subfamily is the largest subfamily in both Arabidopsis and maize. It's noteworthy that except the canonical XXHRLD-type RluAs, another three conserved RluA variants, including XXNRLD-, XXHQID-, and XXHRLG-type were also identified in those key nodes of vascular plants. Subcellular localization analysis of representative AtPUSs and ZmPUSs in each subfamily revealed that PUS proteins were localized in different organelles including nucleus, cytoplasm and chloroplasts. Transcriptional expression analysis indicated that AtPUSs and ZmPUSs were differentially expressed in various tissues and diversely responsive to abiotic stresses, especially suggesting their potential roles in response to heat and salt stresses. All these results would facilitate the functional identification of these pseudouridylation in the future.